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ABSTRACT 

This paper provides a detailed analysis of the main component of the slow neu¬ 
tron capture process (the s-process), which accounts for the solar abundances of half 
of the nuclei with 90 < A < 208. We examine the impact of the uncertainties of the 
two neutron sources operating in low-mass asymptotic giant branch (AGB) stars: the 
13 C(a, n) 16 0 reaction, which releases neutrons radiatively during interpulse periods 
(,kT ~ 8 keV), and the 22 Ne(a, n) 25 Mg reaction, partially activated during the con¬ 
vective thermal pulses (TPs). We focus our attention on the branching points that 
mainly influence the abundance of s-only isotopes. 

In our AGB models, the 13 C is fully consumed radiatively during interpulse. In this 
case, we find that the present uncertainty associated to the 13 C(a, n) 16 0 reaction has 
marginal effects on s-only nuclei. On the other hand, a reduction of this rate may 
increase the amount of residual (or unburned) 13 C at the end of the interpulse: in this 
condition, the residual 13 C is burned at higher temperature in the convective zone 
powered by the following TP. 

The neutron burst produced by the 22 Ne(a, n) 25 Mg reaction has major effects on the 
branches along the s path. The contributions of s-only isotopes with 90 < A ^ 204 
are reproduced within solar and nuclear uncertainties, even if the 22 Ne(a, n) 25 Mg rate 
is varied by a factor of two. Improved /3-decay and neutron capture rates of a few key 
radioactive nuclides would help to attain a comprehensive understanding of the solar 
main component. 

Key words: Stars: AGB - Stars 


1 INTRODUCTION 

The classical analysis of the s process (the slow neutron 
capture process) has provided a first phenomenological ap¬ 
proach to inte rpret the solar s distribution by means of an - 
alytical tools llClavton fc Wardlll974l : iKappeler et al.lll989l) . 
but it was soon clear that the solar s process requires a 
multiplicity of neutron exposures, intuitively associated to 
different astrophysical sites. 

* E-mail: bisterzo@ph.unito.it; sarabisterzo@gmail.com 


The existence of the main component of the s-process was 
advanced by the classical analysis to reproduce the solar 
abundances of s isotopes between 90 < A ^ 204. Low-mass 
Asymptotic Giant Branch (AGB) stars (M < 3 Mq) were 
recog nised to be a most promising site fo r the main compo¬ 
nent (lUlrichl Tl973l : ITberi fc Ren zi n il 1 1 9831 ) . The contribution 
of two additional s-process components was required to re¬ 
produce the solar s-process distribution, the so-called weak 
and strong components. The weak component is partly re¬ 
sponsible for the nucleosynthesis of s nuclei with A < 90 
during the hydrostatic evolutionary phases of massive stars 
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jAmett_&JThielemann 19991 : iLimongi. Straniero fc Chieffil 
120001 : Rausche r et al.ll 2002l). The stron g component was pos- 
tulated bv IClavton & Rassbachl (| 19fi7h to explain about half 
of solar 208 Pb, despite the astrophysical site was not identi¬ 
fied. Successively, the strong component found a natural ex¬ 
planation in AGB star s with low metallicity ([Fe/Hl < — 1) 
and low initial mass (Gall ino et al.l 1 19981 : iTravaglio et al.l 

l200lll . 


This paper is focused on the study of the main compo¬ 
nent. 

Compared to the classical analysis, the development of the 
first AGB stellar models has provided a more adequate de¬ 
scription of the dynamic environment in which the main 
component takes place. The s-process nucleosynthesis in 
AGB stars occurs during the late stages of the stellar evo¬ 
lution, when the star has a degenerate C-O core, a thin 
radiative layer (He-intershcll) and an expanded convective 
envelope. During the AGB phase, the star experiences a se¬ 
ries of He-shell flashes called Thermal Pulses (TPs) triggered 
by the sudden activation of the 3a process at the base of the 
He-intershell. In such a region, free neutrons are released by 
two key reactions, 13 C(a, n) 16 0 and 22 Ne(a, n) 25 Mg. 

The 13 C(a, n) 16 0 reaction is the main neutron source in 
low mass AGB stars. The 13 C forms in a thin zone (the 13 C 
pocket) via proton captures on the abundant 12 C ( 12 C(p, 
7) 13 N(/3 + j/) 13 C). The most favorable conditions for the for¬ 
mation of the 13 C pocket occur during the period imme¬ 
diately following a Third Dredge-Up (TDU) episode. Dur¬ 
ing a TDU, the H-burning shell is switched off and, thus, 
the convective envelope can penetrate inward and carry to 
the surface the heavy elements previously synthesised in the 
He-intershell. The mechanism triggering the formation of 
the 13 C pocket is far from being understood: extant models 
assume that a partial amount of protons may be diffused 
from the convective envelope into the He- and C-rich radia - 
tive He-intershell (see discussion bv IStraniero et all [20061. 
otherwise 13 C would be further converted to 14 N via pro¬ 
ton captures, mainly acting as a neutron poison of the s 
process via the 14 N(n, p) 14 C reaction. Different mechanisms 
proposed to explain the formation of the 13 C pocket are ob¬ 
jects of studjfl When the temperature becomes larger than 
0.8 x 10 8 K (kT ~ 8 keV, which corresponds to Tg ~ 0.9), 
13 C(a, n) le O burns radiatively for an extended time scale 
(some 10 4 yr), releasing a neutr on density N n ~ 10 7 cm” 3 
with a large neutron exposur^EI ( S traniero et al .111 9951 ). 

The second neutron source, the 22 Ne(a, n) 25 Mg reaction, 
is partially activated at the bottom of the convective shells 
generated by the TPs. Starting from the large amount of 
4 N left in the H-burning shell ashes, 1 2 Ne is produced via 


1 See, e.g., the opacity-induced oversh ooting at the b ase o f 

the convective envelope obtained by ICristallo et all ( 2009! ) 
by introducing in the model an exponentially decaying 
profile of the convective velocity. Other models investi¬ 
gate the effects of diffusive ov ershooting, r o tation , mag¬ 
netic fields or gravity waves jHerwi g et' al.l 1 19971, 12003 ; 
Lajjffi3£_et__ayj^^ 2003; S iess etal] 2004 ; 

Piersanti, Crist alio. Sz Stranier j|2013l ; IBusso et al .11 20121 , see also 

Majorca et alj|2012l for results in young open clusters). 

2 The neutron exposure r is the time-integrated neutron flux, r 
= f N n v t h (It. where N„ is the neutron density and Vth the 
thermal velocity. 


the 14 N(a, 7) 18 F(,0 + i') 18 O(a, 7 ) 22 Ne nuclear chain. At a 
temperature of Tg = 2.5 - 3 (kT ~ 23 keV corresponds to 
Tg ~ 2.7), the 22 Ne(a, n) 25 Mg reaction starts releasing neu¬ 
trons, giving rise to a small neutron exposure with a high 
peak neutron density (A n (peak) ~ 10 10 cm” 3 ). Although 
this second burst accounts only for a few percent of the to¬ 
tal neutron exposure, it regulates the final abundances of the 
s-only isotopes nearby to important branch points of the s 
process. 


lArlandini et ali (| 19991 ) provided a first interpretation 
of the solar main component by adopting an average be¬ 
tween two AGB stellar models with initial masses of 1.5 
and 3 Mg, half-solar metallicity, and a specific 13 C-pocket 
c hoice (called case ST). The case ST was calibrated by 
iGallino et al.l (il998l) in order to reproduce the solar abun¬ 
dances of s-only isotopes between 96 ^ A ^ 204 with half¬ 
solar metallicity models. The two 1.5 and 3 Mq models were 
chosen as they represent the stellar mass range that bet- 
ter reprod uces the observations of peculiar s-rich disk stars 
dBusso et al.lll995h . 

The complex dependence of the s process on the ini¬ 
tial ch emical composition of the star was investigated by 
IGallino et al.1 (ll998l) . They found that the strong compo¬ 
nent derives from AGB stars of low metallicity. Indeed, for 
any 13 C-pocket strength the number of free neutrons per 
iron seed increases with the 13 C/ 56 Fe ratio, and the neutron 
fluence progressively overcomes the first two s peaks at neu¬ 
tron magic numbers N = 50 and 82, directly feeding 208 Pb 
(explaining the previously introduced strong component). 
The heterogeneity of the s process is also evidenced by spec¬ 
troscopic observations in different stellar populations: the 
discovery of the first three lead-rich low-metallicity stars 
confirms that 208 Pb may be s trongly produc ed in peculiar 
objects of the Galactic halo (IVan Eck et al.l[200ll) . More¬ 
over, the s elements observed in peculiar s-rich stars (e.g., 
MS, S, C, Ba, and Post-AGB stars, planetary nebulae, and 
lead-rich stars, later called CEMP-s stars) show a scat¬ 
ter at a given metallicity. Thi s scatter has been recog ¬ 
nised since the fi rst st udies by Smith & L ambertl (1990), 
IBusso et al.l (1200 if) . and lAbia et al.l (120021) , (f or rece nt_anal- 
ysis s ee the review by Sneden_et_aL 20081 : iKappeler et al.l 
l201ll : iKarakas fe Lattanzioll2014l . and references therein). 

At present, stellar models are not able to reproduce the 
observed scatter without employing a free pa rametrisati on 
in modelling the formation of the 13 C-POcket (|Herwig et al.l 


19971. 120031: Karakas &; Lattanzid l2007l : ICristallo et al.l 20091 : 

ierol 12013 ). Most uncertainties 
of stellar models are indeed rela ted to the treatment 
of convective /ra d iative interfaces (llben fe Renzinil Il983l : 
iFrost fe Lattanz~iolll996l ). which influence the extension and 
the 13 C profile of pocket, as well as the occurrence and deep¬ 
ness of the TDU. 

Another key uncertainty of AGB stellar models results from 
the unknown efficiency of the mass loss rate, which regu¬ 
lates the number of thermal pulses and the AGB lifetime. 


This is particularly challenging in intermediate-mass AGB 
llVentura fc D’Antonal 12003 : IVentura fc Marigol l201ph . De¬ 
spite the solar s distribution between 90 < A 204 re¬ 
ceives a dominant contribution by low-mass AGB models of 
disc metallicity, intermediate-mass or low-metallicity AGB 
models are crucial to study the chemical evolution of dwarf 
galaxies and globular clusters showing a clear s process 
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signature llTolstov. Hill fe Tosil 120091 ; IStraniero et al.l 12014 
and references therein), or the spectroscopic observation o f 
pec uli ar G ala ctic stars (e.g., Rb-rich stars , van Raai et ali 
120121 ; iKarakas. Garcfa-Hernandez. fc Luearol 20121 : low- 
metallicity s-en hanced stars as e . g., CH, CEMP-s , 
post-AGB stars, IBisterzo et al.l 1201 ll : iLugaro et al.l 120121 : 
iDe Smedt et al.ll20l4 . 

Besides the aforementioned uncertainties, nuclear and solar 
abundance uncertainties affect the s process. Their impact 
may be substantially reduced by using the unbranched s- 
only 150 Sm as a reference isotope of the s distribution (be¬ 
cause it has well known solar abundance and accurately de- 
termined neutron cap ture cross sections of nearby nuclei; see 
lArlandini et al.ll 1999h . 


The s-process abundances observed in the Sun are 
the result of the complex Galactic chemical evolution, 
which accounts of the contribution of different stellar 
generations, with various masses, metallicities and s 
process strengths. A Galactic chemical evolution (GCE) 
model is therefore needed to interpret the dynamics of 
the s process o ver the Galactic history up to the present 
epoch (see, e.g., Tnwaglio^t_aL[ 2004 ;_ Romano et al.ll20ld : 


iNomoto. Kobavashi, fc Tominaeal 20131 . and references 
therein ). 

iTravaglio et al.l (i2004h have shown that in the GCE context 
low mass AGB stars in the range of 1.5 to 3 Mq provide 
the dominant contribution to the solar system main and 
strong components. 

The impact of intermediate AGB stars (IMS; p B ~ 
4 to 8 Mq) on solar abundances is marginal, with the 
exception of a few neutron-rich isotopes ( 86 Kr and 87 Rb 
due to the branches at 85 Kr and 86 Rb; 96 Zr affected by 
the branch at 95 Zr; see also IBisterzo et aI1l2014l f. Indeed, 
the He-shell of IMS stars reaches higher temperatures 
than low mass AGB stars, and the 22 Ne(a, n) 25 Mg 
reaction is efficiently activated, producing higher peak 
neutron densities (Ts ~ 3.6-3.7; N n ~ 10 1 


Under these conditions, the 22 Ne 
becomes the major neutron source 

a, n) 25 Mg reaction 

Truran & Iben 19771; 

Karakas et al. 20061; Longland et al. 

2013: van Raai et al. 

20121; Karakas. Garcfa-Hernandez. & Lugaro; 20121: 

D’Orazi et al. 20131; Straniero et al. 

2014 Dohertvet al. 

2014)1. Otherwise, in the more massive AGB models the 


formation of the 13 C pocket and the occurrence of efficient 
TDUs may be inhibited by hot bottom bur ning and hot 
third dredge-up event s ( Karakas fc Lattanzid 120031 : iHerwid 
|2004 iGorielv fc Siessl l204 l. This significantly reduces the 
overall IMS contribution to the solar s distribution. 

GCE s predictions yield a plausible agreement with the 
solar abundances of the s-only isotopes between 134,136 Ba 
and 208 Pb. The solar distribution between 90 ^ A < 
140 is instead underestimated by our GCE model: an 
additional (unknown) 20-30% contribution is required 
dTravaglio et al.l 120041 '). The light element primary pro¬ 
cess (LEPP) postulated for this contribution is different 
from the s process in AGB stars and from the weak 
s-process component occurring in massi ve stars. Its origin 
is lar gely debated in literature (s ee, e.g.. lFrischknecht et al.l 
120121 and IPignatari et al.l 120131 f or a primary s com po- 
nent in massive stars, Arcones fc Th ielemann 120131 and 


lHansen. Andersen, fc Christliebll2014i for up process or weak 
r process induced by explosive stellar nucleosynthesis). 

For an exhaustive discussion about the s-process nu ¬ 
cleosy n thesis we ref e r to t he reviews by W allcrstcin et_ah 
119971. iBusso et ail (Il999l4 . IHerwid (120051). Istraniero et al. 
12006). iKappeler et al. (1201 llf and Karakas fc Lattanzic 
1 201411 . 


It is evident from the above considerations that the 
main component approach does not provide a realistic inter¬ 
pretation of the whole solar s abundances, because it does 
not account of all complex aspects of the s process over the 
Galactic evolution. 

Nevertheless, it is noteworthy that the solar s contributions 
predicted by main component and GCE model are compa- 
rable in the atomic mass re gion between 140 ^ A ^ 204 (see 
IBisterzo et al.ll201ll . 12014 ). Moreover, despite the remark¬ 
able discrepancy estimated by GCE and main component 
between 90 < A < 140, fairly similar s contributions are 
derived by looking at the relative isotopic predictions of a 
given element. The only exceptions are the rarest neutron- 
rich isotopes ( 86 Kr, s, Rb, 96 Zr), which may be significantly 
produced in intermediate-mass AGB stars. 

Therefore, the study of the main component can be consid¬ 
ered a useful tool to investigate the nuclear aspects and the 
sensitivity of the solar s predictions between 90 ^ A ^ 204, 
focusing on s isotopes close to the major branches of the s 
path. 


The main goal of this paper is to analyse how the un¬ 
certainties of the a-induced reactions on 13 C and 22 Ne may 
affect the main component. 

Both 13 C(a, n) 16 0 and 22 Ne(a, n) 25 Mg reaction rates may 
be affected by contributions of sub-threshold states and res¬ 
onances, making the extrapolation of the laboratory cross 
section measurements do wn to the energy range of stellar 
burning a complex task dWiescher, Kappeler fc Langankj 

l2012h . 


The 13 C(a, n) ls O rate may be influenced by the unknown 
contribution of the sub-threshold resonance at 6.356 MeV 
(owing to the M = l/2 + state in 11 0). U p to a fact or of three 
uncertainty was evaluated bv lAnguIo et al.l (|l999l f . Over the 
years, s everal analyses have been dedicated to this neutron 
source jDroyeff_etalJJl99^_ Denker^talJJl995|^ Haldll99^; 


Kubono et al.l 2003[~ Keelev et al.ll2003l : I Johnson et al 2006 : 


Pellegriti et al. 21 


Recent investigations ha ve signifi¬ 
cantly i mproved the accu racy (down to ~20-30%: lHeil et al.l 


l2008bl : I Guo et al 


.2012). Lately, the abo ve resonance has 


La Cognata et all (|2013i f with the Trojan 


been detected by 
horse method. 

The 22 Ne(a, n) 25 Mg rate is dominated by the well studied 
resonance at 832 keV. The 22 Ne(a, n) 25 Mg uncertainty at Ts 
= 3 is mainly due to the unknown influence of a resonance at 
633 keV. The large upper limit estimated by lAngulo et al.l 
d 1999tl (up to a factor of fifty) has been reduced by more 
than a factor o f ten with the experimental measurement 
bv I Jaeger et alJ d200lll . Remarkabl e experimental and the ¬ 
oretical works were carried ou t by Kappeler et al.l d 19941) . 
iKoehlerl (|2002lL Karakas et al.l d2006l ~). lUgalde et al.1 d2007Tl 
and Longland et aid ( 2012| j . 


The study of the main component is useful to investi¬ 
gate the effects of the major nuclear network uncertainties 
on AGB s predictions, especially close to the branch points. 
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A branch in the s path occurs when neutron captures com¬ 
pete with /3 decays at unstable isotopes with half-lives of a 
few weeks to a few years, so that the s path branches partly 
by /3 decays following the usual path in the stability val¬ 
ley and partly by neutron captures feeding the neutron-rich 
isotopes. The strength of the branching is described by the 
branching factor fp=\p/(\p+\ n ), where A/ 3 =l/r=ln 2 /t 1 / 2 
is the /3-decay rate corresponding to the half-life t \/2 (or to 
the mean lifetime r). The decay pattern may include the /? , 
/3 + and electron capture channels. The neutron capture rate 
X n =N n < av >=N„vt< a > is given by the neutron den¬ 
sity N n , the mean thermal velocity vt, the relative velocity 
between neutron and target v, and the Maxwellian-averaged 
cross section (MACS) defined as < a >=< av >/vt- The 
neutron capture strength of the branching point is defined 

aS fn —An /(\/3 T \n ) • 

The uncertainties in the experimental (n, 7 ) rates have 
been highly reduced in recent years reaching in some cases 
a precision smaller than one or two percents (see, e.g., 
iKappeler et al.l l201lh . However, the high s-process tem¬ 
peratures allow the low-lying excited states to be popu¬ 
lated by the intense and energetic thermal photon bath. 
Because only the ground state is accessible by experi¬ 
ment, the effect of neutron captures in excited states has 
to be eval uated theoretically and suffers from large un¬ 
certainties (iRauscher fe Thielemann! 120001 : lRauscherll2012l : 
IReifarth, Lederer, fe Kappelen [2014). Even major uncer¬ 
tainties are associated to unstable nuclei for which no ex¬ 
perimental measurement is available at the present time. 

In addition, the /3-decay rates of some radioactive isotopes 
may be largely affected by variations of temperature and 
electron density. Altho ugh t he labo ratory /3-decay rates 
are accurately known ICvburt et al.l 120100 . the contribu¬ 
tion of thermally populated excited levels and the effects 
of unknown transitions in a strongly ionised plasma can 
largely modify the /3-dec ay rates at stellar temperatures. 
iTakahashi fe Yokoil (ll987l l investigated the /3-decay rates of 
unstable heavy isotopes at temperatures and electron densi¬ 
ties typical of stellar interiors (5x10' ^ T ^ 5x10 s K; 10 26 
^ Tie ^ 3xl0 27 cm -3 ), finding large deviations from the ter¬ 
restrial values. The temperature dependence of branchings 
is even more complex, if they have isomeric states that are 
thermalised at high temperatures through transitions via 
mediating states at higher excitation energy. Consequently, 
the abundances of the affected s-only isotopes carry direct 
information on the physical conditions occurring during the 
s process, i.e., neutr on density, temperature, and density 
(IKappeler et al.lll989l ). 

In this context, the branch point isotopes 151 Sm, 163 Dy, 
179 Hf, 176 Lu are remarkable examples. The /3 - -decay rate 
of 151 Sm increases strongly at He-shell flash tempera- 
ture and regulates the abundances of t he pair 152,154 Gd 
dMarrone et a1.ll200t3 : IWisshak et al .1120061) . 163 Dy and 179 Hf 
are terrestrially stable, but start to decay at temperatures 
typical of TPs producing a non-negligible s-process contri¬ 
bution to 164 Er and 1 S 0 Ta m , both usually bypassed by th e 
s path llJaag fe KappeleJll996l : IWisshak et al.ll 200 ll . |2004) . 
Because of its long half-life, 176 Lu was long considered as a 
cosmo-chronometer until the decay rate was found to_exhibit 
a very strong temperature depen dence dKlav et al.l Il99ll : 
iHeil et ai1l2008d : iMohr et al.ll2009l '). 

The /3-decay rates calculated in stellar environments are sub¬ 


j ect t o nuclear uncertainties difficult t o estimate. [G oricly 
1999) recomputed the calculations of Takahas hi fe Yokoi 


198711 by assuming a typical uncertainty of Alog ft = 


±0.5 for the decay rates of unknown transitions and showed 
that the final stellar rates may vary by up to a factor 
of three at typical s-process conditions (Tg = 3 and n e 
= 10 27 cm -3 ). Recently, the nuclear NETwork GENera- 
tor fNE TGEN; IXu et al. 2013fl ex trapolated the /3-decay 
rates by Thkahashi fe Yokoil 1 1987ll and their uncertainties 
dGorielvl 1999l l to an extended range of temperature and 
electron density. In a few cases the /3-decay properties of ex¬ 
cited states have been measured in laboratory experiments, 
e.g. for long-lived isomeric states ( 180 Ta m determined by 
iBelic et alJ I 2 OO 2 I; or the key branch for the weak s pro¬ 
cess in massive stars 79 Se m studied by iKlav fe Kappelen 
Il988l) or for isotopes with small Q-values, which are sen¬ 
sitive to the bound-state /3-dec ay mechanis m (examples 
are 163 Dy a nd 187 Re studied by l Jung et al.l 1 19921 and by 
iBosch et al .11 1 99fil . respectively). In general, however, mea- 
surements o f suc h d ata rema in a big experimental challenge 
(jReifarth, Lederer. fe Kappeleill2014h . 


The paper is organised as follows. We start with a de¬ 
scription of the solar main component, updated with a neu¬ 
tron capture network that includes the most recent cross 
se ction measurements a s well as the solar abundance data 
bv lLodders et all (120091 1 . (Section [3]). The AGB models em¬ 
ployed ffi_this work have been outlined by iBisterzo et al.l 
1 I 2 OI 0 L120141') . 

In Sections [3] and 5J we discuss the impact of the 13 C(a, 
n) le O and 22 Ne(a, n) 25 Mg uncertainties on the branches of 
the main component. 

Relevant branches of the s path and their influence on the s 
nuclides are analysed in Section [5] We distinguish the s-only 
isotopes in different classes, according to the characteristics 
of the related branchings. First, unbranched s-only isotopes 
(e.g., 100 Ru, 104 Pd, 110 Cd, 116 Sn, 124 Te, 150 Sm, 160 Dy, and 
198 Hg) are useful to constrain the main component (Sec¬ 
tion ED- Second, we address branchings where the decay 
of the unstable branch point isotopes exhibits no effect of 
the s-process temperature (e.g., 96 Mo, which is sensitive to 
the branch at 95 Zr) or where the half-lives are essentially 
independent of temperature. These branchings may provide 
information on the s-process neutron density (e.g., 170 Yb, 
partially regulated by 169 Er; 142 Nd, marginally affected by 
142 Pr; 186 Os, which mainly depends on the branch at 185 W; 
192 Pt, influenced by 192 Ir; see Section E2J. Third, branch¬ 
ings initiated by decays with pronounced dependencies on 
temperature and/or electron density may be interpreted as 
s-process thermometers (e.g., 134 Ba, owing to the branch 
at 134 Cs; 152 ’ 154 Gd, strongly sensitive to 151 Sm and 154 Eu; 
176 Lu, which has a short lived isomer; 204 Pb, because of the 
branch at 204 Tl; Section 15.311 or indicators of the turnover 
timescale of convective mixing during TP (e.g., 128 , 130 Xe, 
164 Er, 180 Ta m ; Section l5.4ll . Each of the three classes is illus¬ 
trated in detail by a few selected examples: 150 Sm, adopted 
to normalise the s distribution (Section 15.1.11) . as well as 
96 Mo and 1,0 Yb (Sections 15.2.11 and 15.2.21) . 134 Ba (Sec- 
tion l5.3.ril . and 180 Ta m ('Section l5.4.1ll . which are most sensi¬ 
tive to uncertainties related to the cross sections and /3-decay 
rates of the involved branch point nuclei. Other branchings 
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Figure 1 . Top panel : stellar (n, 7 ) cross sections at 30 keV for nuclides with A > 20. The compilation is taken from the KADoNiS 
database (v0.3): small dots and big circles refer to experimental and theoretical MACS, respectively. Recent MACS not included in 
KADoNiS, but considered in this work are indicated by squares (experimental measurements) and triangles (theoretical evaluations); 
see text for details. Neutron-magic numbers at N = 50, 82 and 126 are highlighted. Bottom panel: MACS uncertainties at 30 keV. Colors 
of error bars are the same as in the top panel. 


are briefly summarised, referring to the complete discussion 
in Section B (Supporting Information). 

In Section C (Supporting Information), we describe the 
effect of the 13 C(a, n) 16 0 , 22 Ne(a, n) 25 Mg and 22 Ne(a, 
7 ) 26 Mg uncertainties on two AGB models (a 3 Mq model at 
[Fe/H] = — 1 and a lialf-solar metallicity 5 Mq model), which 
are selected as representative of a more extended range of 
stars than usually adopted for the main component. 

An overview of the results is provided in Section [ 6 ] 


2 UPDATED SOLAR MAIN COMPONENT 

Despite the main component does not provide a realistic 
description of the solar s distribution between A = 90 and 
204, it still represents a useful approximation to investigate 
the effect of nuclear uncertainties of s isotopes in this atomic 
mass range. 

As shown bv lArlandini et ahl (Il999l ). the ’’stellar” main 
component is obtained as the averaged yields between two 
AGB models with initial masses of 1.5 and 3 A/q at half-solar 
metallicity. The s-process nucleosynth esis is com puted with 
the post-process method described bv lGallino et al.l d 1998T) . 
which is b ased on input da t a of f u ll evo lutionary FRANEC 
models by IStraniero et al. | (j 19971 . 120031 ). like the temporal 


history of the temperature and density gradients during the 
convective TPs, the number of TPs, the mass of the H shell 
and He intershell, the overlapping factor, and the residual 
mass of the envelope. 

In both AGB models, the 22 Ne(a, n) 25 Mg reaction oper¬ 
ates in a similar range of temperatures (Tg ~ 2.5-3). How¬ 
ever, the M = 3 Mq model achieves Tg ~ 3 for the last 
16 TPs, while the M = 1.5 Mq model for the last 8 TPs. 
The marginal activation of the 22 Ne(a, n) 25 Mg neutron 
burst mainly occurs in the advanced TPs for a rather short 
timescale (~6 yr). 

Most of the s-process nucleosynthesis takes place radiatively 
during the interpulse phases via the 13 C(a, n) le O neutron 
source, which lasts for 3-6 xlO 4 yr (this time decreases with 
increasing pulse number). We assum e a single 13 C - pocke t 
profile close to case ST described by iGallino et alj (Il99a) . 
The case ST 13 C-pocket contains 5xl0 -6 Mq of 13 C and 
2xl0 -7 Mq of 14 N, and extends in mass for ^lxlCP 3 Mq 
(~1/10 of the typical mass involved in a TP). 

Th e mass los s is estimated with the Reimers formula 
(R eiiner 3E9 z 3) with the parameter 77 = 0.3 for 1.5 Mq and 
r/ = 1 for 3 Mq. This allows the occurrence of 19 and 25 TPs 
with TDU, respectively. Updated opacities and a revised lu¬ 
minosity function of Galactic car bon stars (iLederer et al.l 
120091 : iGuandalini fe Cristalloll2013h suggest a more efficient 
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mass loss than the one adopted in our models. This signifi¬ 
cantly reduces the number of TDUs (e.g., ~15 TDU for the 
above models: ICristallo et al.ll 201 1 ). To this purpose, note 
that the overall s yields of our AGB models are marginally 


influenced by the contribution of the s abundances predicted 
in the envelope for TDU numbers higher than 15. Moreover, 
the normalisation of the average s yields to that of 150 Sm al¬ 
lows us to overcome some of the major uncertainties of AGB 
stellar models (e.g., the efficiency of the TDU, the mass loss, 
the number of TPs), which would otherwise dominate the 
s-process predictions (see Section [T]). 

For major details on the A GB models employed we refer to 
iBisterzo et al.l J2010l . 12014 ). 


The neutron capture cross section network is based on 
KADoNiS f|. Different MACS ar e adopted for ^ °' 21 ’ 22 Ne by 
iHeil et al](l20l4. 24,25,2 6 Mg bvlMassimi et al.l (|2()I2|). 41 K 
and 4 5 Sc by Heil et all J2009I), 58 ' 6 2 ^ 3 ]\Ji by l^iugec et al.l 
d20l4 and jLederer et ahl (120131 . I20l4. 64,70 Zn b y 
iReifarth et all 1 I 2012 IU soTs2? 3.84,86 Kr bv lMutti et all ll2005l) 
85 Kr by Raut et al.l J2013I) . 92,94,96 Zr by Thgliente et al l 
j20Rl_ 20111 14. "^"Zr by iTaeliente et all (1201311 an d 
iLuearo et al.l ll20l4. 136,13 7 Cs b y Patronis et al.1 (l2004. 
i48,i49p m b v |R e jf ar th et all (120031') . 152,154 Eu by Bestetall 
d200ll') . the p- only 180 W by iMarganiec et ahl 1 201(f) . 
!86 ,i87,i88q s |Mosc oni et al. Ml . 192,194,195,196 p t fay 

iKoehler fc Gubei ( 20131 ') . and their theoretical evaluation of 


the 192 Ir MACS. 

For a few heavy isotopes, the stellar temperatures are high 
enough for a significant population of low-lying nuclear ex¬ 
cited states, which modify the neutron capture cross sec¬ 
tions measured in laboratory. These effects are considered 
in the stellar MACS with the so-called stellar enhancement 
factor (SEF), which accounts of the stellar average over the 
thermally populated st ates and the laboratory cross section. 
The SEF estimated bv lBao et alj (t2000l l are updated with 
those recommended by KAD oNiS. An excep tion is 187 Os, 
for which we have adopted iFuiii et all (l2010| j, who studied 
the inelastic channel of the neutron capture cross section 
of 18 'Os and employed a Hauser-Feshbach statistical model 
theory (HFSM) to obtain a reliable estimate of the SEF. 

The state of the art of all MACS at kT = 30 keV ver¬ 
sus atomic mass is shown in Fig. [T] (top panel). We dis¬ 
tinguish between experimental measurements (small dots) 
and theoretical evaluations (big circles). Updated MACS 
with respect to KADoNiS are indicated by big squares (ex¬ 
perimental measurements) or triangles (theoretical evalu¬ 
ations). The lower panel of Fig. [l] represents the related 
percentage M ACS uncertaintie s. _ 

As outlined bv lKappeler et ahl (1201 il l, most of the MACS in¬ 
volved in the nucleosynthesis of stable isotopes heavier than 
A ~ 90 are known with an accuracy <5%. Uncertainties of 
~10% are evaluated for 104 Pd (10%), 139 La (9.6%), 159 Tb 
(9.5%; near the s-only isotope 160 Dy), 166 , 167 , 168 p] r (10 to 
13%), 200 Hg (10%), 209 Bi (11.7%). 

Among the unstable isotopes, the largest MACS uncertain¬ 
ties are quoted for 85 Kr (^50%; recently evaluated from the 
inverse 86 Kr( 7 , n) 85 Kr reaction a nd from the 86 Kr( 7 , 7 ’) 
measurement bv lRaut et alj|2013h . 86 Rb (~80%), 95 Zr (up 


3 Karlsruhe Astrophysical Data Base of Nucleosynthesis in Stars, 
web site http://www.kadonis.org/, version KADoNiSv0.3. 


Mass average; case ST; [Fe/H] = -0.3 



Figure 2. The solar main component versus atomic mass is repro¬ 
duced by assuming a ST 13 C-pocket, and by averaging between 
AGB models with initial masses of AT = 1.5 and 3 Mq at half solar 
metallicity. In this Figure, the s-production factors are normalised 
to that of the unbranched s-only nucleus 150 Sm. The s-production 
factors in the He-intershell (given in mass fraction ’Xi over the 
solar-scaled initial values) correspond to the material cumula¬ 
tively mixed with the envelope by the recurrent TDU episodes 
and eventually dispersed in the interstellar medium by efficient 
stellar winds. Different symbols are used for 86,87 Sr (which receive 
additional s contributions, e.g., the weak component in massive 
stars), and for 2(13 Pb, ~50% of which is produced by low-mass, 
low-metallicity AGB stars (strong s process). The error bars dis¬ 
played for the s-only isotopes accou n t for t he uncertainties of the 
solar abundances by iLodders et ahl (i200 9 ). Note that the solar 
uncertainties are not shown for 128 Xe, 152 Gd, 164 Er, and 180 Ta m 
because they receive additional non-negligibl e contr ib utions from 
the (^)jo-process by SNIa and SNII (e.g., iTravaglio et al.l [20141 ; 
iRauscher et alj|2002t IPignatari et al.ll20 1 3a[) . The solid line cor¬ 
responds to a 100% contribution from main component, dashed 
lines represent 10% of uncertainty. 


to a factor of two uncertainty. iLugaro et al . 203 ), 141 Ce and 
142 Pr (~40%). Other unstable isotopes with uncertainties 
from 15% to 25% are 79 Se (important for the weak s-process 
component), 81 Kr, "Mo, 103 Ru, 110 Ag, 115 Cd, 147 Nd, 160 Tb, 
170 Tm, 179 , 182 Ta, 181 Hf, 186 Re, 191 Os, 192 Ir, 193 Pt, 198 Au, 
203 Hg, 204 T1, 20 S Pb. 

The solar main component is shown in Fig. [2] We exam¬ 
ine the s-production factors in the He-intershell (correspond¬ 
ing to the overabundances with respect to the solar-scaled 
initial values given in mass fractions ’AV) of isotopes with 
A > 80 normalised to the unbranched s-only nucleus 150 Sm. 
A plausible reproduction of s-only isotopes heavier than A 
= 90 (full circles) is obtained within th e sola r abundance 
uncertainties quoted by ILodders et_alj _(|2009lb Varia tions 
with respect to p revio us results by Bisterzo et, al.l ll201ll ) and 
IKappele r et al.| (|2011 ) partly d erive from new solar abun- 
dances bv lLodders et alj (120091 ) and partly from updates in 
the nuclear reaction network. Among the nuclides with A > 
90, a noteworthy variatio n is t he incr eased solar abundance 
of Hg estimated by ILodders et al.l (120091 ') (+35%), which 
reduces the 198 Hg solar s-contribution by iKappeler et ali 
d201ll) from ~100% to 8 2%. The slight l y larg er solar Sn 
un certainty estimated by Lodders et al.l 1 20091) (from 10% 
bv I Anders fc Grevessel[T989l to 15%) agrees better with the 
s prediction for 116 Sn, although the solar abundance is still 
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13 % higher. Note that t he so lar Nb and Xe values estimated 
by iKashiv et al.l d2006h and Reifarth et al] (l2002fl were al¬ 
ready considered in lKapneler et al.l I 201ll '). 

The updated network produces small variations (<5%) 
among the s-only isotopes. Differences larger than ~10% 
are obtained for 164 Er and 180 Ta m , where t he already ascer¬ 
tained dominant s-p rocess contributions dBest et al.ll200ll : 
IWisshak et al.ll200ll ') have been increased: 


• The s contribution of solar 164 Er is enhanced by 8% 
(from 74% to 82%), because a smaller SEF is evaluated for 
the MACS of 164 Er (the old SEF = 1.24 at 30 keV estimated 
bv lBao et al.|[2000l is reduced to 1.08 by KADoNiS). 

• 1S0 Ta m increases from 75% to 82% (+7%), mostly due 
to a decreased SE F of the MACS for 180 Ta m (from SEF = 
0.96 at 30 keV bv lBao et alJ[200ol to 0.87 by KADoNiS). 


Starting from these updated results, we analyse in the 
following sections the effect of the two AGB neutron source 
uncertainties on the s distribution and discuss the conse¬ 
quences of the MACS and /3-decay rate uncertainties of the 
most relevant branching points of the main com ponent. 

Recently, t heore tical analyses bv lRauscher et all (1201 lh 
and iRauschen ll2012h showed that the SEF correction can 
cause larger MACS uncertainties than assumed so far. This 
holds especially for cases, where the (experimentally mea¬ 
sured) ground state cross section constitutes only a minor 
fraction of the MACS. If the ground state contribution, ex¬ 
pressed by a new factor A', is close to unity, the uncertainty 
of the stellar rate is directly connected to the experimental 
one; otherwise, the uncertainty of stellar rate may be larger 

than that of the experimental measurement._ 

The t heoretical uncertainties discussed by iR.auscher et al.1 
ll201lll are not included in the present work, because they 
should be considered as upper limits of future more realistic 
analyses, which need to be carried out individually for each 
nucleus, by accounting of all available nuclear details. The 
situation is even more challenging when the experimental 
rate is inferred from the in verse neutron c apture reaction, 
(y,n) (e.g., the case of 185 W: lRauscheijf2014l l. On the base of 
the branches analysed in Section [5] and in Section B (Sup¬ 
porting Information), we highlight outstanding isotopes that 
would need specific theoretical investigations in Section [6] 


3 UNCERTAINTY OF THE 13 C(a, n) le O 
NEUTRON SOURCE 

As anticipated in Section [T] the 13 C(a, n) le O rate is the 
dominant neutron source in low mass AGBs and shapes 
the overall s-abundance component. It burns radiatively 
during the interpulse periods in the top layers of the He- 
intershell when the temperature reaches Tg ~ 0.8-1.0. Un¬ 
der these conditions, the 13 C(a, n) le O rate is extremely 
small (10~ 14 1(T 13 cm 3 /mol s) and, therefore, not yet ac¬ 
cessible to direct measurements in terrestrial laboratories, 
owing to the high background signals induced by cosmic 
7-rays. Current reaction rates have been determined by ex¬ 
trapolation of direct measurements performed at energies 
higher than 270 keV, far from the Gamow window at 140- 
230 keV whe re t he 13 C(a, n) 16 0 reaction is effective (see, 
e.g.. lDrotleff et ahlll993l'l . This entails an unknown influence 
of broad sub-threshold states. A major impact is expected 


Table 1. The 13 C(ct, n) le O rate at T = lx 10 s K (i n units 


Denker et al. (Il995l), Kubono 

et alj (i2003h,lJohnson et al 

(2006 

L 

Heil et alJ (]2008bl), Guo et al. (20121), NACRE 

II (X 

u et al.l 

2013al) and La Cognata et al 

(20131). Additional studies 

are not 

listed explicitly, but found 

13 C(a, n) le O rates 

within these 

ranges. 





Reference 

Rate 

Lower 

Upper 


(IQ- 14 cm 3 /mol s) 

Limit 

Limit 

Caughlan & Fowler 11988) 

2.58 

- 

- 


Kubono et al. 12003) 

2.02 

1.49 

2.54 


Johnson et al. 12006) 

2.64 

2.11 

3.30 


NACRE II 

4.86 

3.62 

6.48 


Denker et al. 11995) 

4.32 

- 

- 


Heil et al. 12008b) 

4.6 

3.6 

5.6 


Guo et al. 120121 

4.19 

3.46 

4.90 


La Cognata et al. 12013) 

6.20 

5.32 

7.13 



from the contribution of the sub-threshold resonance due to 
the U = l/2 + state in 17 0 at 6.356 MeV (E)f b = —3 keV), 
while minor effects are expected due to sub-threshold reso¬ 
nances at 5.939 MeV (U = 1/2“; = —547 keV) and 

5.869 MeV (U = 3/2+; E l £ b = -641 keV). 

Past analyses showed large uncertainties at 8 keV (up to 
a factor of three) owing to t he unexplored contributio n of 
the 6.356 MeV state (see e.g.. ICaughlan fe Fowledll988l and 
lAngulo et al.l 1 19991). M o re rec ent evaluations by Heil et aD 
( 2008bh and IGuq et al.l (l2012l l are in reasonable agreement 
and provide a significantly reduced uncertainty (down to 
~20-30%, in the AGB energy range of interest). Later, 
lLa Cognata et ahl (l2013l l attained an accuracy of 15% for 
the 13 C(a, n) 16 0 rate by means of the Trojan horse method, 
which allows one to study reactions of astrophysical inter¬ 
est free of Coulomb suppression and electron screening at 
astrophysical energies with no need for extrapolation. This 
method clearly displays the presence of the sub-threshold 
resonance corresponding to the 6.356 MeV 17 Q sta t e. Thei r 
rate is about 35% higher than that by iHeil et al.l (l2008b|j . 
still in agreement within uncertainties with most of the re¬ 
sults in literature at Tg ~ 1. 

In Table [T] we provide a summary of the most relevant 
13 C(a, n) le O rates found in literature at Tg = 1. 

In our AGB models we adopt the 13 C(a, n) 16 0 rate by 
Denker et al. (Il995l l. w hich is pre t ty clo se to the results by 


Heil et akl ~i 2008bf Tand ICuo et all (l2012h . 


We carried out two extreme tests in order to analyse the 
effect of the 13 C(a, n) ls O uncertainty on the main compo¬ 
nent: 


• Test 1 - Adopted 13 C(a, n) le O ra te multiplied by a 
factor of two, close to the upper limit bv lLa Cognata et al.l 

Uni); 

• Test 2 - Adopted 13 C(a, n) ls O rate divided by a factor 
of three, close to the lower limit bv lKubono et all (I2003T) . 

The results of the two tests are shown in Fig. Q5] for 
isotopes with A > 80. We focus on s-only isotopes for clarity 
(circles). The ratio of the results obtained in the tests and 
with the adopted 13 C(a, n) le O rate are plotted in the upper 
and lower panels together with the 10% limits indicated by 
dashed lines. 

We see that the main component is marginally affected 
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s-only isotopes with 90 < A < 205: 
main s-process component 



Atomic mass 


Figure 3. Ratios between the main component obtained with a 
two times higher and a three times lower 13 C(a, n) le O rate than 
our adopted rate shown in Fig[2](Test 1 and Test 2 correspond¬ 
ing to the top and bottom panels ). We focus on s-only isotopes 
for clarity. The shaded area between 80 < A < 90 indicates the 
atomic mass region affected by additional s contributions (e.g., 
the weak s process in massive stars; IMS AGB stars, see Section 
B and C, Supporting Information). The two neutron-magic nuclei 
86 Kr and 87 Rb, which exhibit the largest variations, are marked 
by crosses (see text). A complete version of this Figure (which 
includes all isotopes from 90 ^ A 7 210) is given in Section A 
(Supporting Information). 

by both tests: variations are smaller than 1% for the s-only 
isotopes with A > 90 and smaller than 3% between A = 80 
and 90. 

Significant changes were obtained, however, for the two 
neutron-magic nuclei 86 Kr and 87 Rb (represented by 
crosses ), because the branching at 85 I<r is partially acti¬ 
vated during the 13 C(a, n) 16 0 phase. For a comprehensive 
discussion see Section B (Supporting Information). We note, 
however, that both isotopes receive a small contribution 
from AGB stars: only ~16% and ~18% of solar 86 Kr and 
87 Rb are produced by the main component, respectively. 

In our AGB models with initial masses 1.5 and 3 Mg 
and half solar metallicity the amount of 13 C in the pocket 
assumed for case ST burns radiatively during the interpulse 
even by reducing the 13 C(a, n) le O rate by a factor of three. 
Only a 13 C mass fraction negligible for the s process is in¬ 


gested in the next convective TP. 

Otherwise, when a substantial amount of 13 C is ingested in 
the subsequent TP, it burns at the bottom of the convective 
shell (at Ts ~ 1 . 6 ) in a very short time scale (of the order of 
few years), producing enhanced neutron densities of up to 
~10 10-11 cm -3 . A new gene ration of AGB model s avail able 
at the FRUITY database bv lCristallo et ahl (l2009l . [201 ll) ex¬ 
periences a partial convective burning of 13 C during the first 
(and second) TP(s). This mainly occurs in low-mass metal- 
rich AGB models and influences the abundances of a few 
neutron-rich isotopes during the first TPs (e. g., 86 Kr, 87 R b, 
96 Zr: and radioactive nuclides as 60 Fe, see ICristallo et al.l 

l2006U201ll ). 

Actua lly, the recent measurement by lLa Coenata et all 
HH) suggests that the 13 C(a, n) le O rate adopted in our 
models must be increased by 40%, rather than decreased 
(see Table |T|). In this case, the ingestion of 13 C into TPs is 
strongly reduced. 

Despite the main component presented in this work is 
marginally influenced by the uncertainty of the 13 C(a, n) ls O 
rate, AGB predictions in general may show larger sensitiv¬ 
ity, when e.g., a different amount of 13 C is assumed in the 
pocket, or a lower initial metallicity is adopted. Moreover, 
the physical characteristics of different stellar models em¬ 
ployed in the nucleosynthesis calculations may influence the 
effect of the 13 C(a, n) 16 0 rate. For instance, recent low- 
mass AGB models are character ized by more efficient TDU 
episodes (see the discussion in IStraniero et al.1 [200(0) ■ In 
particular, the first TDU occurs at the very beginning of 
the TP-AGB phase, when the core mass is rather small and, 
in turn, the temperature developed in the He -rich i nt ersh cll 
during the interpulse period is lower (ICristallo et alJfeoogh . 
In this case, a not negligible amount of 13 C that survives 
after the interpulse is engulfed in the convective zone gener¬ 
ated by the subsequent TP and burns at relati vely high tem¬ 
peratu re. For example , in 1 .8 Mg models bv lLuearo et al.l 
( 120121 ) and I Guo et al.l (l2012h . the convective burning of the 
partial 13 C ingested in the TPs affects s predictions up to 
Pb. This effect may even become dominant for lower mass 
AGB models (e.g., in M^p B = 1.25 AI@ models of half-solar 
metallicity the 13 C neutron source burns convectively rather 
than radiatively; iR.aut et al.|j2013h . 


4 UNCERTAINTY OF THE 22 Ne(a, n) 25 Mg 
NEUTRON SOURCE AND THE 22 Ne(a, 

7) 26 Mg RATE 

While the 13 C(a, n) 16 0 neutron source determines the over¬ 
all shape of the main component, the contribution of the 
22 Ne(a, n) 25 Mg reaction affects the s-only isotopes close to 
branching points. 

For the 22 Ne(a, n) 25 Mg reaction, the center-of-mass energy 
region of interest to the s process is E c . m . = 300 to 900 
keV. At these low energies, direct 22 Ne(a, n) 25 Mg measure¬ 
ments are limited by cosmic 7 -ray background because the 


4 A more efficient TDU is required to reproduce the observed lu¬ 
minosity function of t he Carbon stars in th e Milky Way and i n the 
Magellanic Clouds llCristallo et al.l 1201 ll : [Guandalini fe Cristallol 
120131) . 
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cross section rate is extremely small. The lowest well stud¬ 
ied r = 2+ resonanc e at ~832 keV , whic h dominates the 
22 Ne(a, n) 25 Mg rate (iDrotleff et al.1 fl993h . lies above the 
relevant astrophysical energies. Theoretical extrapolations 
of the 22 Ne(a, n) 25 Mg reaction measured at higher energies 
may be affected by the unknown influence of low-energy res¬ 
onances just below the neutron threshold. 

The resonances affecting the (a, n) and (a, 7) rates of 22 Ne 
correspond to different levels in the compound nucleus. First 
investigations derived a large 22 Ne(a, n) 25 Mg uncertainty 
because of the influence of a possible resonance at ^635 keV 


Angulo et al. 

19991: KaDDeler et al. 1994: Jaeger et al. 200ll: 

and Koehler 

2002). Successively, Longland et al. 

200§) de- 


termined the energy and quantum numbers of excited states 
in 26 Mg through a nuclear resonance fluorescence experi¬ 
ment and demonstrated that the corresponding level at an 
excitation energy of E x = 11 154 keV has unnatural parity 
= 1 + , contrary to the previous spin-parity assignment. 
Because both 22 Ne and 4 He have = 0 + , by angular mo¬ 
mentum selection rules only natural-parity (0 + , 1“, 2 + , etc.) 
states in 26 Mg can participate in the 22 Ne(a, n) 25 Mg reac¬ 
tion, thus excluding the = 1 + resonance. Nevertheless, 
the 22 Ne(a, n) 25 Mg rate may be affected by the possible con- 
tributi on of other unknown low-threshold sta tes of natural 
parity (lUealde et al.ll2007l ; IdeBoer et al.ll2010l f . Moreover, if 
yet unknown low-energy resonances have a strong influence 
on the 22 Ne(a, y) 26 Mg rate, the 22 Ne(a, 7) 26 Mg and 22 Ne(a, 
n) 25 Mg reactions may compete during TPs, further affecting 
th e neutro n production in AGB stars (see, e.g., IMS stars, 
iKarakas et al.ll2006l ). 

The recent theoretical estimate by iLongland et al.l (120121 ) 
based on a Monte Carlo evaluation suggests a much smaller 
22 Ne(a, n) 25 Mg uncertainty (~20%). 

_ In previous AGB models (e.g., iBisterzo et al.l l2010l. 

Hon]), we have adopted the lower limit of the 22 Ne(a , 
n) 25 Mg and 22 Ne(a, 7) 26 Mg rates bv lKappeler et ahl ll 19941) . 
These rates were evaluated by neglecting the resonance con¬ 
tribution at 633 keV. At Ts = 3, the 22 Ne(a, n) 2 3 Mg rate was 
about 50% higher than that reco mmended by Jaeger et al.l 
(l200lll . close to the upper limit bv lLongland et alj ( 2012h . 

In Table [5] we summarise some of the most significant 
results achieved in the estimate of the 22 Ne(a, n) 25 Mg and 
22 Ne(a, 7) 26 Mg rates for the maximum temperature at the 
bottom of the advanced TPs in the AGB models adopted 
here (Ts = 3). 


We have calculated the rates of the 22 Ne(a, n) 25 Mg 
and 22 Ne(a, 7) 26 Mg reactions (’’This work” in Table0 see 
also Fig. 0 by accounting for the most recent studies of all 
known (directly and/or indirectly detected) resonances as 
well as nuclear data available in literature. 

In particular, we adopted the rates calculated on the basis of 
experimentally detecte d levels, i.e. mainly on the basis of the 
values of I Jaeger et al.l d200ll '). For the 22 Ne(a, 7) 26 Mg rate 
we added the states which have an indirect d eterminati on 
of the a width using the 22 Ne( 6 Li, d) reaction dGiesen et alj 
11993I J. In the temperature range of low mass AGB stars (Ts 
= 2.5-3), the 22 Ne(a, n) 25 Mg recommended value is about 
a factor of two lower than that adopted in our models so far, 
while the recommended 22 Ne(a, 7) 26 Mg rate is essentially 
unchanged. 


Table 2. The 22 Ne(a, n) 25 Mg and 22 Ne(ct, 7 ) 26 Mg rates 
at 3xl0 8 K (in un i t of 10 —11 cm 3 /mol s) estimate d 
b\^JCau ghlan_&: Fowled j l988|), NACR E lAngulo_et_al] [l999|), 

I Kag peleret*aLl (ll994l V. [j aeger et al. j200ll) . iLongland et al.l 
(2012). The newly evaluated values are labelled with ’’This work” 
(see text for explanations). We distinguish between experimental 
measurements (exp) and theoretical evaluations ( th ). 


Reference 

Recommended 

Rate 

(10 -11 cm 3 /mol s) 

Lower 

Limit 

Upper 

Limit 

22 Ne(o!, n) 25 Mg rate: 

Caughlan &; Fowler (1988): (th) 

1.86 

- 

- 

NACRE; (th) 

4.06 

3.37 

192.0 

KaDDeler et al. ( exv ") 

9.09 

4.14 

14.4 

Jaeger et al. (2001); (exp) 

2.69 

2.63 

3.20 

Longland et al. (2012); (th) 

3.36 

2.74 

4.15 

This work; (th) 

2.24 

1.99 

2.92 

22 Ne(o:, 7 ) 26 Mg rate: 

Cauehlan & Fowler (1988)} (th) 

0.46 

- 

- 

NACRE; (th) 

2.56 

0.59 

20.30 

KaDDeler et al^^)94); (exv) 

1.22 

0.81 

1.63 

Longland et al. (2012); (th) 

1.13 

0.93 

1.38 

This work; (th) 

0.80 

0.72 

2.62 


We have also determined lower and upper values of both 
rates considering additional resonances in order to investi¬ 
gate their potential impact (see Fig. 0 bottom panel). The 
lower limit is determined by the lower experimental limit 
of the 703 keV resonance (11%). This reduces the recom¬ 
mended rates by about 10% at all temperatures. For the 
upper rate we introduced all the known (but not directly 
measured) states below the 703 keV resonance, adopting 
a spectroscopical factor of 0.01, which represents a conser¬ 
vative upper limit. Moreover we have estimated a possible 
larger effect due to one of the low energy states. We then 
concentrated on the lower energy states, which give substan¬ 
tial contributions to the reaction rate for T < 0.3 GK. This 
was done to investigate the influence of possible a cluster 
states in 26 Mg. These states could have spectroscopic fac¬ 
tors as large as 0.1. Larger values are unlikely because they 
should have been observed in a transfer reactions. For the 
( a , 7) reaction we identified two states, one above (568 keV) 
and one below (433 keV) the neutron threshold. Also for the 
(a, n) reaction there are two states: one at 497 and another 
at 548 keV. This corresponds to an increase of the recom¬ 
mended 22 Ne(a, n) 2S Mg recommended rate by a factor of 
two at Ts = 2.5 and by +30% at Ts = 3. The recommended 
22 Ne(a, 7) 26 Mg rate increases up to a factor of seven at Ts 
= 2.5 and by a factor of three at Ts = 3 (see dashed lines 
in Fig. 0 bottom panel). Note that in network calculation 
one should not use the upper rate for both 22 Ne(a, n) 25 Mg 
and 22 Ne(a, 7) 26 Mg reactions at the same time because the 
change of spectroscopic factors has influence on both rates. 

Based on the above results, we have updated the so¬ 
lar main component discussed in Section 0 with the rec¬ 
ommended 22 Ne(a, n) 25 Mg and 22 Ne(a, 7) 26 Mg rates (see 
Fig- 0 filled circles). By comparing these revised results 
with the previous s-distribution displayed in Fig. 0 major 
variations are shown close to the branches of the s path. 
The impact of the uncertainties evaluated for the recom¬ 
mended 22 Ne(a, n) 25 Mg and 22 Ne(a, 7) 26 Mg rates (Fig. 0 
has a rather small effect on the s-distribution (see thick bars 
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T (x 10 9 K) 


Figure 4. Top panel : new evaluated values of the 22 Ne(a;, n) 25 Mg 
and 22 Ne(a, 7 ) 26 Mg rates with black and red lines respectively. 
The difference between lower and recommended rates are almost 
unnoticeable (about 10% at all temperatures; dotted lines). Our 
upper rates are shown by dashed lines. See text for details on the 
computed uncertainties. Bottom panel : The uncertainty bands for 
the 22 Ne(a, n) 25 Mg and 22 Ne(a, 7 ) 26 Mg reactions. The uncer¬ 
tainties are the result of upper limit resonance contributions and 
of resonance strength uncertainties. The shaded areas at Tg = 
2.5-3 indicates the temperature range reached during TPs by low 
mass AGB stars. Explicit values of both rates at Tg = 3 are given 
in Table [2] 


in Fig. 0. The abundances of s-only isotopes are reproduced 
within the solar uncertainties. 

However, it should be noted that the several ambiguities re¬ 
maining in the nuclear data for the 22 Ne + a reaction rates 
need to be resolved. Indeed, the theoretical 22 Ne(a, n) 25 Mg 
and 22 Ne(a, 7) 26 Mg rates presented in this work (and all 
estimations given in literature as well) ar e b ased o n old ex¬ 
perimental measurements (e.g., Ijaeger et al.ll200ll ) and the 
uncertainties may be underestimated. The presence of low 
energy unknown states, which have been identified in several 
indirect experiments, makes the evaluation at relevant en¬ 
ergy still uncertain. New direct experimental investigations 
are needed to measure the resonances with higher accuracy. 
Future investigations are in progress to shed light on this 
issue (e.g., via ind ir ect measurements a t n_TOF-EAR2, 
IChiaveri et al. ] 120121 . iMassimi et al.l 120141 : or at accelera¬ 
tor facilities deep underground where the cosmic-ray back¬ 



Figure 5. The solar main component versus atomic mass (as 
shown in Fig. 0 obtained with the recommended 22 Ne(a, n) 25 Mg 
and 22 Ne(a, 7 ) 26 Mg rates (see Fig. [Tj. We focus on s-only iso¬ 
topes for clarity (filled circles). Also 86,87 Sr (which receive ad¬ 
ditional s contributions) and 208 Pb (not an s-only nuclide) are 
displayed for completeness. Thin err or-bars account for th e uncer¬ 
tainties of the solar abundances bvILo dders et al . (120091) . Thick 
bars are obtained by including the uncertainties of the recom¬ 
mended rates discussed in the text (lower and upper limits) and 
shown in Fig. [4] We indicate the results of an improved treat¬ 
ment of the half-life of a few key isotopes strongly sensitive to 
temperature (and electron density), which influence the nearby 
isotopes ( 134 ’ 136 Ba, 152,154 Gd, 164 Er, 180 Ta m and 204 Pb; filled 
diamonds , see text). 


ground into detectors is reduced by several orders of magni- 
tude, LUNA0 and DIANA0). 

Our post-process AGB models generally adopt constant 
/3-decay rates, which are based on a geome tric average of the 
rates given by iTakahashi fe Yokoil rtl987l l at different tem¬ 
peratures (and electron densities) over the convective pulse. 
Note that this approximation does not affect the prediction 
of most of the s-only branched isotopes, because the /3-decay 
rates of close unstable nuclei are almost constant in the tem¬ 
perature range of the TP, or they do not compete with neu¬ 
tron capture rates. Relevant exceptions are 134 Ba, 152 ’ 154 Gd, 
164 Er, 176 Lu/ 176 Hf, 1S0 Ta m , 204 Pb: in these cases, the ,0- 
decay rates of close unstable isotopes (e.g., 134 Cs, 151 Sm and 
154 Eu, 164 Ho, 176 Lu, 180 Ta m , 204 T1, respectively) vary by or- 
der(s) of magnitude over the large temperature and density 
gradients that characterise the convective zone (0.2 < Tg < 
3 and 10 < p < 10 4 g/cm 3 ), competing with neutron cap¬ 
tures. The above branches require an improved treatment of 
the /3-decay rates over the TP. 

The treatment of the branches close to 176 Lu/ 176 Hf and 
1S0 Ta m was already refined in recen t AGB models (see 
iHeil et aPl2008 d : IWisshak et al.ll2006al l. We extend the im¬ 
provement to the branches close to 134 Ba, 152,154 Gd, 164 Er 
and 204 Pb. We further implement the s prediction of 180 Ta m , 
by including the same treatment to nearby branched iso¬ 
topes ( 179 Ta and 179 Hf). 

We provide detailed calculations by dividing each TP in 30 


5 See ICostantini et all ll2009l f and references therein; web 
http: //l una.lngs.inf n.it/ 

D SeelL emut et aid ll201ll ). 
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convective meshes of constant mass. In each mesh temper¬ 
ature (and density) can be considered constant during each 
time step. We follow the production and destruction of the 
unstable isotopes close to the above s-only nuclides in each 
mesh with the /3-decay rates computed at each mesh tem¬ 
perature. The neutron density resulting from the 22 Ne(a, 
n) 25 Mg reaction is computed in the various meshes at each 
time step. Neutron densities for an efficient s process are 
only reached in the bottom region of the convective zone 
(see Fig. Al, Section A, Supporting Information). The abun¬ 
dances resulting in each zone after neutron irradiations are 
periodically mixed to account for the turnover time of the 
convective zone (of the order of a few hours). 

As shown in Fig. H 134 Ba, 152 ’ 154 Gd, 164 Er and 180 Ta m are 
mainly affected by the improved treatment of the /3-decay 
rates (see filled diamonds); a variation smaller than 5% is 
displayed by 204 Pb. We will provide a more detailed discus¬ 
sion on these branches in Section [5] and Section B (Support¬ 
ing Information). 


4.1 Tests of the 22 Ne(a, n) 25 Mg and "Ne(a, 

7 ) 26 Mg rates 

As discussed in Section |4j future direct experimental mea¬ 
surements of the 22 Ne(a, n) 25 Mg and 22 Ne(a, 7 ) 26 Mg reac¬ 
tion rates may evidence the presence of low energy unknown 
states. Also the contribution of known resonances needs to 
be determined with high precision. 

These considerations encourage us to work in a more con¬ 
servative range of uncertainty than that estimated in Sec- 
tion[4] starting from our recommended 22 Ne(a, n) 25 Mg and 
22 Ne(a, 7 ) 26 Mg rates (hereafter adopted to compute the ref¬ 
erence s distribution; see Fig. El), we investigate the effects 
of the following tests on the solar main component: 

• Test A - Recommended 22 Ne(a, n) 25 Mg rate multi- 
pli ed by a factor o f four , close to the recommended value 
bv lKappeler et ahl (1 1994T) . while the 22 Ne(a, 7 ) 26 Mg rate is 
almost unchanged. 

• Test B - Recommended lower limits of the 22 Ne(a, 
n) 25 Mg and 22 Ne(a, 7 ) 26 Mg rates. 

• Test C Recommended 22 Ne(a, n) 25 Mg rate and up¬ 
per limit for the 22 Ne(a, 7 ) 26 Mg rate. This test evaluates 
the competition between (a, 7 ) and (a, n) rates in AGB 
models of low initial mass. 

I 11 Fig. [5] we show the effect of Test A ( upper panel) 
and Test B ( lower panel) on the solar main component 
of s-only isotopes with A > 80. In both panels, a 10% un¬ 
certainty is indicated by dashed lines. A detailed list of the 
results of both tests that includes all isotopes with A ^ 70 
is given in Table Al, (Supporting Information). 

Values are obtained by normalising the s-production fac¬ 
tors to that of lso Sm in both tests. Note that for 150 Sm we 
obtain an s-production factor of X( 150 Sm)/Xi rl i( 150 Sm) = 
1133.9 with the recommended 22 Ne(a, n) 25 Mg and 22 Ne(a, 
7 ) 26 Mg rates. Test A provides X( 150 Sm)/X;„i( 150 Sm) = 
1275.5 (+12.5%), and Test B yields X( 150 Sm)/X ini ( 150 Sm) 
= 1128.4 (—1.0%). The s production of 150 Sm increases by 
+ 12.5% with increasing the 22 Ne(a, n) 25 Mg by a factor of 
four. This difference is independent of any branch (see Sec¬ 
tion 15.1.11) and affects the entire s-process distribution in 
the same way. 
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Figure 6. Ratios between the main component obtained with 
Tests A and B and our results obtained with the recommended 
22 Ne(ct, n) 25 Mg and 22 Ne(a, 7 ) 26 Mg rates (Fig[TJ. We focus on s- 
only isotopes for clarity (filled circles). An improved treatment 
of the branches close to 134 Ba, 152,154 Gd, 164 Er, 180 Ta m and 
204 Pb is included (filled diamonds). A complete version of this 
Figure (which includes all isotopes from 90 Si A ^ 210) is given in 
Section A (Supporting Information). As in Fig. [3] the shaded area 
between 80 < A < 90 indicates the atomic mass region affected 
by additional s contributions. 


As displayed by Test A, the 22 Ne(a, n) 25 Mg rate mainly 
affects 86,87 Sr, which decrease by ~30% being regulated by 
the branches at 85 Kr and 86 Rb. We remind that the main 
component synthesises about half of solar 86,87 Sr (see dis¬ 
cussion in Section B and C, Supporting Information). 
Noteworthy variations (^10-20%) are obtained for 96 Mo, 
134 Ba, 142 Nd, 152 Gd, 170 Yb, and 176 Lu. The s prediction to 
96 Mo is affected by the branch at 95 Zr (Sections 15. 2. Ill : the 
large sensitivity of the branch at 134 Cs to neutron density 
modifies the 134 Ba/ 136 Ba ratio (up to 10%; Section 15.3.11) : 
142 Nd is regulated by the branch at 141 Ce (Section B, Sup¬ 
porting Information); the branches at 151 Sm and lo4 Eu in¬ 
fluence the s production of lo2,154 Gd (e.g., variations of 20% 
of 152 Gd; Section B, Supporting Information); 170 Yb is af¬ 
fected by the branch at 170 Tm (10%; Section T5. 2.2 II : finally, 
the branch at 176 Lu modifies the 176 Lu/ 176 Hf ratio (see dis¬ 
cussion in Section B, Supporting Information). 

Note that the s contribution to 100 Ru and 104 Pd also shows 
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Figure 7. Ratios between the main component obtained with the 
recommended 22 Ne(a, n) 25 Mg and 22 Ne(a, 7 ) 26 Mg rates com¬ 
pared to Test C in which the upper limit for the 22 Ne(a, 7 ) 26 Mg 
rate is adopted, while the 22 Ne(«, n) 25 Mg reaction is unchanged. 
Symbols are the same as Fig. [6] A complete version of this Fig¬ 
ure (which includes all isotopes from 90 ^ A ^ 210) is given in 
Section A (Supporting Information). 

important variation, although the branches at "Mo and 
103 Ru (with strongly reduced half-lives at stellar temper¬ 
atures) are only marginally open. In both cases, the ~16% 
variation is mainly explained by the cumulative effect of a 
different s contribution to 150 Sm, adopted to normalise the 
s distribution. 

Other s-only isotopes show differences of less than 5%. 
A comprehensive description of the most important branch¬ 
ings of the main component is provided in Section [5] 

Note that the variations shown by Test B are marginal be¬ 
cause our recommended 22 Ne(a, n) 25 Mg rate lies close to 
a plausible lower limit (~ —10%; Section 0. Indeed, the 
dominant effect of the resonance at 800 keV indicates that 
smaller values are unlikely. 

In our AGB models with initial masses 1.5 and 3 
Mq and half-solar metallicity, the competition between the 
22 Ne(a, 7 ) 26 Mg and 22 Ne(a, n) 25 Mg reactions is marginal: 
in spite of the rather generous margins considered for the 
22 Ne(a, 7 ) 26 Mg rate, Test C confirms that the effect of the 
22 Ne(a, 7 ) 26 Mg rate for the s-abundance distribution is neg¬ 
ligible (see Fig. 0). 

Appreciable variations are only seen for 26 Mg, directly in¬ 
volved in the (a, 7 ) reaction. Due to its very small MACS 
(0.126 ± 0.009 mb at 30 keV. lBao et ahlfePOOl ). 26 Mg is ac¬ 
cumulated in the s process. 

In stellar models with higher initial mass (e.g., IMS 
stars), where the temperature at the bottom of the TP 
increases sufficiently to efficiently activate both 22 Ne(a, 
n) 25 Mg and 22 Ne(a, 7 ) 26 Mg reactions, their competition be- 
com es important for the nucleosynthesis of Mg isotopes (see 
e.g.. iKarakas et ahll2006l . iLongland et al1l2012l) . Above all, 
the 22 Ne(a, n) 25 Mg rate plays a key role in the production of 
s isotopes, as the 13 C(a, n) 16 0 neutron source is expected to 


have small or negligible effects in IMS stars (see Section 0. 
In Section C (Supporting Information), we discuss the im¬ 
pact of the 22 Ne(a, n) 25 Mg and 22 Ne(a, 7 ) 26 Mg rates on 
two AGB models: a 5 Mq model at half-solar metallicity 
taken as example of IMS stars and a 3 Mq model at [Fe/H] 
= — 1 chosen as representative of low-metallicity models. 

In a 5 Mq model at half-solar metallicity, the maximum 
temperature at the bottom of the convective TPs is Tg ~ 
3.6. In this condition, the 22 Ne(a, n) 25 Mg neutron burst ef¬ 
ficiently operates, reaching a neutron density of N n ~ 10 11 
cnF 3 . The whole s distribution increases by up to a fac¬ 
tor of ~2 by including the upper limit of the recommended 
22 Ne(a, n) 25 Mg rate, and by one order of magnitude with 
Test A (recommended [ 22 Ne(a, n) 25 Mg ]x4). Key neutron- 
rich isotopes ( 86 Kr, s 'Rb, 96 Zr) are largely produced being 
the branches at 85 Kr and 95 Zr easily open. The impact of 
the 22 Ne(a, 7 ) 26 Mg rate on s isotopes remains rather small 
(< 6 %). 

In low-mass AGB stars with [Fe/H] < —0.3 the 22 Ne(a, 
n) 25 Mg uncertainty also produces a larger impact on s pre¬ 
dictions than that observed on the solar main component. 
For instance, in a 3 Mq model at [Fe/H] = — 1 the maximum 
temperature at the bottom of the advanced TPs reaches Tg 
= 3.5. In this case, both 22 Ne(a, n) 25 Mg and 13 C(a, n) 16 0 
neutron sources operate efficiently. In addition, the s dis¬ 
tribution is largely modified by the lower initial metallicity: 
as discussed in Section 0 by decreasing the metallicity the 
abundances of isotopes with neutron magic numbers N = 50 
and 82 are progressively overcome (thus reducing the whole 
distribution between 90 < A < 130 and 140 < A < 204, 
respectively), while 208 Pb is mainly produced. 

We refer to Section C (Supporting Information) for major 
details. 


5 UNCERTAINTIES OF MAJOR BRANCHES 
OF THE MAIN COMPONENT 

We distinguish the s-only isotopes in three classes, according 
to their related unstable branch point isotopes: 

• unbranched s-only isotopes (Section l5.Hl . with unstable 
isobars having half-lives shorter than a couple of days (thus 
forbidding neutron captures during TPs); 

• s-only isotopes sensitive to neutron density only (Sec¬ 
tion [50, with unstable isobars having half-lives (almost) 
constant in stellar environments; 

• s-only isotopes affected by branch point with unstable 
isobars having /3-decay rates quickly changing under stel¬ 
lar conditions: we distinguish between nuclides sensitive to 
both neutron density and stellar temperature (and/or elec¬ 
tron density; Section [50 and nuclides less affected by neu¬ 
tron density, but dominated by stellar temperature and/or 
electron density gradients during TP (Section 15.411 . 

5.1 Unbranched s-only isotopes 

As mentioned in Section 0 a few s-only isotopes of the 
s path ( 100 Ru , 110 Cd, 116 Sn, 124 Te, 150 Sm, 160 Dy, see 
lArlandini et al.l 1 19991 : as well as 104 Pd and 198 Hg) are 
marginally affected by nearby branch points with short half- 
lives. Thus, only a few percent of the s path (<3%) bypasses 
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Figure 8. Schematic representation of the branches close to the 
s-only isotopes 148 ’ 150 Sm (red squares). Thick lines represent the 
s-process nucleosynthesis during the 13 C(ct, n) le O neutron irradi¬ 
ation, while thin lines correspond to the neutron capture channels 
open by the marginal activation of the 22 Ne(a, n) 25 Mg reaction 
during TPs. Major branches at 147 Nd and 148 Pm, which regu¬ 
late the s contribution to 148 Sm, are highlighted by circles. All 
branches of the s path join at 150 Sm, which is adopted to nor¬ 
malise the overall s distribution. ( This and the following Figures 
are adapted from http://www.kadonis.org /. See the electronic edi¬ 
tion of the Journal for a colour version of the Figures .) 

the above s-only nuclides, which become useful constraints 
for the main component. 

5.1.1 The s-only isotope 150 Sm 

150 Sm is one of the few s-only isotopes that is exposed to 
the full s flow, because all branches in the Nd-Pm-Sm region 
join at 149,150 Sm (see Fig. [SJ. Indeed, the short half-lives of 
149 Nd (ti /2 = 2.21 d) and 150 Pm (t 1/2 = 2.68 h) leave 150 Sm 
virtually unbranched. 

As anticipated in Section O 150 Sm is particularly suited to 
normalise the overall s-process predictions. We select 150 Sm 
because, besides being an unbranched s-only isotope, 150 Sm 
has very well known MACS (~1% uncertainty, <r[ 150 Sm(n, 
7 )] = 422 ± 4 mbarn; KADoNiS) and solar abundance 
(5% unce rtainty, 0.265 ± 0.013 number of atoms per 10® 
Si atoms: iLodderset al.ll2009l) . 

Unlike 150 Sm, the 148 Sm abundance and, thus, the 
148 Sm/ 150 Sm ratio are otherwise regulated by the branches 
at 147 Nd, 147 Pm, and 148 Pm. 

During the 13 C-pocket phase, the neutron density is not suf¬ 
ficient to bypass 148 Sm at the branching points 147 Nd (fi /2 
= 11 d) and 148 Pm (ti / 2 = 5.37 d) efficiently, whereas the 
branch at 14 'Pm (ti /2 = 2.62 yr) is partially activated (at 
N n = lxlO 7 cm -3 , ~20% of the s path feeds 148 Pm, by¬ 
passing the long-lived 147 Sm). 

On the other hand, the 147 Nd and 148 Pm neutron capture 
channels are activated during the 22 Ne(a, n) 25 Mg irradi¬ 
ation. At the beginning of the TP, the neutron density 
strongly increases and 148 Sm is largely bypassed by the s 
path. Under these conditions, the amount of 148 Sm pro¬ 
duced during the previous 13 C-pocket phase is progressively 
depleted but starts to recover as the neutron density falls 
during the TP and is almost completely restored at the end 
of the TP (~94% of the value produced at the end of the 
previous interpulse is reestablished). In Fig. [9j we plot the 
temporal evolution of the neutron density and of the isotopic 


148 Sm • 147 Sm 147 Pm * l49 Pm * 148 Nd -w 



Figure 9. Temporal evolution of the neutron density and of 
the isotopic abundances of 147 d48,i49,i50g m ^ 147,148^ an( j 
147,148,149 p m during the 15 th He shell flash in an AGB star with 
1.5 Mq and half-solar metallicity. The abundance values are given 
as number fractions normalised to 150 Sm at the end of the He shell 
flash. The time-scale starts when the temperature at the bottom 
of the convective TP reaches Ts = 2.5, which corresponds to the 
onset of the 22 Ne(o:, n) 25 Mg reaction. The initial abundances are 
those left from the previous 13 C-pocket phase. 


abundances of l 47 ,148749,150^ 147,148 Nd and 147,148,149p m 
during the 15 th He shell flash for an AGB model with 1.5 
Mq and half-solar metallicity. The hnal 148 Sm abundance 
is predominantly determined by the freeze-out of the neu¬ 
tron supply, intuitively defined by the time whe n an isotopic 
abundance reaches 90% of its Hnal value llArlandini et al.l 
Il999l : ICosner et~aDll980i ). The abundance of a branched nu¬ 
cleus is frozen when the probability of further neutron cap¬ 
tures (for the nearby unstable isotope) is marginal. Thus, the 
starting moment of the freeze-out depends on the MACS of 
the unstable isotope: the larger the MACS is, the later the 
freeze-out occurs. 

The branch at 147 Pm (U / 2 = 2.62 yr, which following 
iTakahashi fc Yokoilll987l is reduced to 1.2 yr at Tg = 3) 
mainly affects the production of 14 'Sm. The remaining 
14 'Pm abundance decays eventually into 147 Sm at the end 
of the TP. 

The main component reproduces the solar abundance 
of 148 Sm (+5% with respect to solar), in agreement within 
the solar uncertainty. 

As for 150 Sm, also the MACS of 148 Sm is accurately 
known (241 ±2 mbarn at 30 keV, KADoNiS; <1%). 

While the short half-lives of 147 Nd and 148 Pm are prohibitive 
for measuring their (n, 7 ) cross sections with present tech¬ 
niques, the MACS of 14 'Pm could be successfully deter¬ 
mined via th e activa tion method (709 ± 100 mbarn, 14%; 
iReifarth et al.l I 2 OO 3 I : KADoNiS). Uncertainties of ~17% 
have been estimated for the calculated MACS values of 
147 Nd and 148 Pm (at 30 keV, cr[ 147 Nd(n, 7 )] = 544 ± 
90 mb, KADoNiS; a[ 148 Pm(n, 7 )] = 1014 ± 175 mbarn, 
IReifarth et al.l[2003l ). The MACS uncertainties of the branch 
point isotopes 147 Nd, 147 Pm and 148 Pm have small effects 
on 148 Sm (<4% in total). 
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Table 3. The s-production factors of 148 Sm and 150 Sm and their relative ratio (which corresponds to the s contribution to solar 148 Sm) 
for different choices of the 22 Ne(a, n) 25 Mg rate. 


s-production factors Recomrn. 22 Ne(a, n) 25 Mg rate ( a ) TEST A TEST B [ 22 Ne(a, n) 25 Mg ]=0 ( b ) 

Recomm. 22 Ne(a, n)x4 Recomm. 22 Ne(a, n)x0.9 


148 Sm 

1194.1 

1286.9 (+7.8%) 

1195.4 (<1%) 

1258.5 (+5.4%) 

150 Sm 

1133.9 

1275.5 (+12.5%) 

1128.4 (<1%) 

1088.4 (-4%) 

148 Sm/ 150 Sm (in %) 

1.05 

1.01 (-4%) 

1.06 (+1%) 

1.156 (+10.1%) 


( a ) Results obtained with our adopted 22 Ne(a, n) 25 Mg rate (the same as Fig. 0. 

( b ) Results with the contribution of the 13 C(a, n) le O neutron source alone. 


Although 13 C(a, n) le O is the major source of neutrons 
in low-mass AGB stars, the 148,150 Sm s abundances are 
slightly influenced by the efficiency of the 22 Ne(a, n) 25 Mg 
reaction: the s production of 150 Sm increases by 4% by in¬ 
cluding the recommended 22 Ne(a, n) 25 Mg rate, and rises 
up to +12.5% with Test A, while the relative contribution 
148 Snr/ 150 Sm shows up to ~4% variations (see Table[3]|. This 
variation reflects a more (or less) efficient s-process contri¬ 
bution during the TPs neutron burst independently of any 
branch, and affects the entire s-process distribution. On the 
other hand, the 22 Ne(a, n) 25 Mg neutron source is crucial for 
regulating the 148,150 Sm abundances: the solar 148 Sm/ 150 Sm 
ratio would be overestimated by ~16% by excluding its ad¬ 
ditional neutron irradiation (see column 5 of Table [3J. 

5.1.2 Additional unbranched isotopes 

Besides 150 Sm (Section [5.1.11) . 100 Ru, 104 Pd, 110 Cd, 116 Sn, 
124 Te, 160 Dy, and 198 Hg are noteworthy s-only isotopes es¬ 
sentially unaffected by branches, and therefore also impor¬ 
tant for characterising the entire s-process distribution. 

At TP temperatures, the /3-decay half-lives of their poten¬ 
tial branch point isotopes are of the order of a couple of 
days, so that their decay rates clearly dominate over the 
respective neutron capture rates (compare e.g., the stellar 
half-liv es of "Mo, 103 Ru 110 Ag , 115 Cd, 122 Sb, 160 Tb, and 
198 Au; iTakahashi fe Yokoilll987l) . 

Accordingly, the variations of the s contributions to 100 Ru, 
104 Pd, and 110 Cd in Fig. [S] include the cumulative effect of 
the normalisation to 150 Sm (Section 15.1. 111 . The most well 
known isotopes are 110 Cd, 124 Te, and 150 Sm, which have 5- 
7% solar uncertainty and MACS values are known at ~1%. 
Less accurate MACS data with uncertainties of 6 and 10% 
are affecting the s prediction of 100 Ru and 104 Pd, although 
their solar abunda n ces ar e well determined. 
lutsunomiva et all (120131 ) have indirectly determined the 
"Mo MACS via the inverse ( 7 , n) reaction. Although their 
estimated value at 30 keV is significantly larger than the 
theoretical MACS recommended by KADoNiS (410 instead 
of 240 mb), the effect on the s prediction for 100 Ru is only 
~ 1 %. 

The s abundances of 116 Sn and 124 Te are reproduced in 
agreement within the solar un certainties (15 and 7%, re¬ 
spectively; iLodders et alif2009l ) . 

The s abundances of 160 Dy and 198 Hg are difficult to 
determine as discussed in detail in Section B (Supporting 
Information). 

About 10% of solar 160 Dy is missing. This dearth may 
be reconciled within the uncertainties: although the (n, 7 ) 


cross section of 160 Dy is very well determined in laboratory 
(u[ 160 Dy(n, 7 )] = 890 ± 12 mbarn at 30 keV, 1.4%; KADo¬ 
NiS), the MACS remains rather uncertain in this case due to 
a significant stellar enhancement factor (SEF = 1.12 at 30 
keV; KADoNiS). This correction may provide the 10% miss¬ 
ing s-process contribution to solar 160 Dy. T he solar abun¬ 
danc e of Dy is instead well determined (5% bv lLodders et all 
l2009h . 

The s contribution to solar 198 Hg is very uncertain, mainly 
because Hg is a volatile elem ent, and its solar a bundance is 
affected by 20% uncertainty (iLodders et al.ll2009i h An addi¬ 
tional 8.7% uncertainty derives from the 198 H g MA CS (173 
± 15 mbarn; dated back to lBeer and Macklinll 19851 : KADo¬ 
NiS). 

Note that the large uncertainties associated to the 116 Sn and 
198 Hg s-contributions derive from the poorly known Sn and 
Hg solar abundances. 

5.2 The s-only isotopes sensitive to the neutron 
density 

In this Section, we analyse the s contribution to 96 Mo, which 
is influenced by the branch at 95 Zr. This branching depends 
only on the neutron density because the /3-decay rate of 95 Zr 
remains constant at the relevant s-process temperatures. 
The s-only isotopes 170 Yb, 142 Nd, 186 Os and 192 Pt belong 
to this class of branchings as well. Although the half-lives 
of their related branch point nuclides exhibit a marginal 
sensitivity to stellar temperature, the respective MACS un¬ 
certainties still dominate over the temperature effects. 

5.2.1 The s-only isotope 96 Mo (the branch at 95 Zr) 

In this atomic mass region, the long-lived 93 Zr behaves as a 
stable isotope during the main s process (t\/i = 1.5x10® yr; 
Fig. [Toll. The s path directly feeds 94 Zr, bypassing 92,94 Mo 
(two p-only isotopes, mainly destroyed during AGB nucle¬ 
osynthesis). The radioactive decay of 93 Zr into 93 Nb occurs 
after the end of the TP- AGB phase, w he n the s-proce ss nu¬ 
cleosynthesis stops fe.g.. IWallerstein fe Dominyill988l ). This 
decay produces a decrease of the relative [Zr/Nb] ratio. Sim¬ 
ilarly to the discovery of Tc in the spectrum of an AGB 
star of spectral type S (lMerrilllll952l : [Lambert et al.lll995h . 
the [Zr/Nb] ratio supplies spectroscopic information about 
the synthesis of heavy elements in e volved stars (see, e.g., 
Ilvans et al. I l2005l : iKashiv et al.ll2Qloh . 

The abundance of 96 Mo is dominated by the branch 
point at 95 Zr (£ 1/2 = 64.03 d). At the low neutron density 
reached during the 13 C-pocket phase (~10 7 cm -3 ) most of 
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Figure 10. The same as Fig. [8] but for the s-path region close to 
the s-only isotope 96 Mo (red rectangle). While 93 Zr is practically 
stable on the time scale of the s process, 95 Zr acts as the main 
branching point. 
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Figure 11. The same as Fig. [9] but for the isotopic abundances 
of 95 ’ 96 Mo, 95 ’ 96 Zr and 95 Nb. 


the s flow proceeds towards 95,96 Mo, but moderate neutron 
captures on 9j Zr are allowed at the peak neutron density 
achieved during TPs (N n > lxlO 9 cm -3 ; f n > 0.1). Un¬ 
der these conditions the s path feeds the neutron-rich 96 Zr 
(about 3% of solar Zr), and 95,96 Mo are partially bypassed. 
Fig. |TT| shows the temporal evolution of the neutron den¬ 
sity and of the isotopic abundances of 95,96 Mo, 95,96 Zr and 
95 Nb during the 15 th He shell flash in an AGB star with 1.5 
Mq and half-solar metallicity. At the peak neutron density 
96 Mo is depleted, but as soon as the neutron density de¬ 
creases, 96 Mo starts rising again to exceed the values at the 
end of the previous interpulse by 18%. 

While 96,97 Mo and 96 Zr are only influenced by the branch 
at 95 Zr, 95 Mo is mostly bypassed by the s path owing to the 
branch at 95 Nb (ti /2 = 34.99 d; N n > 4xl0 9 cm -3 , f n > 
0.5). 


The 95 Zr MACS is largely uncertain, and discrepant 
values are found in l iterature: K ADoN iS recommends the 
value estimated by iBao et ahl (|2000l ) (79 ± 12 mb at 
30 keV), with a rather s mall uncertainty (15%), whereas 
iToukan &; Kappeleij ( 1990ll estimat e d a va lue of 50 mb at kT 
— 30 keV. Recently, Lugaro e t al.l (|2 014T> provided a value 
that was 50% lower than that of Toukan fe Kappeled (ll99Cll 'l 
and about three times lower than in KADoNiS. Accordingly, 


the MACS of 95 Zr is still affected by an uncertainty of about 
a factor of two. 

In our AGB model s , a v alue close to that suggested by 
IToukan fc Kappeleil ( 19901 ) was adopted so far for the 95 Zr 
MACS. 


In Table IT] we list the s contributions to isotopes from 
94 Zr to 98 Mo. By including the old 22 Ne(a, n) 25 Mg rate, 
the s contribution to solar 96 Mo was overestimated by +4% 
(column 2; which include our guess on the 95 Zr MACS); 
this value slightly increases to +8% by adopting t he new 
theoretical 95 Zr MACS evaluated by iLugaro et al.l d2014l) . 
(column 3). 

Updated s predictions (which include the recommended 
22 Ne(a, n) 25 Mg rate and the 95 Zr MACS bv lLugaro et al.l 
12014T) reproduce 99% of solar 96 Mo. Note that the solar 96 Mo 
is largely overestimated (+20%) with increasing the recom¬ 
mended 22 Ne(a, n) 25 Mg rate by factor of four (Test A) . This 
exceeds the solar Mo uncertainty (T0%: lLodders et al.ll2009l l 
and the 96 Mo MACS uncertainty (7%; 112 ± 8 mbarn at 30 
keV; K ADoNiS), which dated bac k to an earlier measure¬ 
ment bv I Winters fc Macklinl (119871 ). 

Major changes are shown by 96 Zr. Although 96 Zr is a 
neutron-rich isotope, it receives an important s contribu¬ 
tion during the 22 Ne(a, n) 25 Mg neutron burst in AGB stars: 
once built-up by the activation of the branch at 95 Zr, 96 Zr 
is only marginally destroyed by neutron ca ptures owing to 
its small MACS (10.3 ± 0.5 mb at 30 keV; iTagliente et al.l 
l2011al ). Starting from lArlandini et al.l (1 19991 ). the main com¬ 
ponent was known to produce about half of the solar 
96 Zr. Present calculations obtained with the recommended 
22 Ne(a, n) 25 Mg rate reduce this contribution from 50% to 
14%. This value is extremely sensible to the 22 Ne(a, n) 25 Mg 
rate: Test A yields an increase by a factor 6 (87%). More¬ 
over, a factor of two uncertainty is still associated to the 
theoretical 95 Zr MACS. 

We remind that 96 Zr is a well known indicator of the 
initial AGB mass. An addi tional s contribution to 96 Zr 
comes from IMS AGB stars dTravaglio et al.l|2004l ). Indeed, 
although IMS stars marginally produce heavy s isotopes, 
the strong neutron density reached during TPs (N n ~ 10 11 
cm” 3 ) may significantly increase the abundance of 96 Zr (see 
Section C, Supporting Information). 


5.2.2 The s-only isotope 170 Yb (the branches at 169 Er and 
170 Tm) 

The abundance of 1,0 Yb is regulated by the branches at 
170 Tm and to a minor extent at 169 Er (whose terrestrial 
h alf-life U /2 = 9.4 d decr eases to ti/ 2 ~ 7.5 d at Tg = 
3: lTakahashi &c Yokoilll987l ). Neutron captures on 169 Er are 
only relevant at high neutron densities (N n > 3xl0 9 cm” 3 ; 
fn > 0.3). 

The s contribution to 170 Yb is essentially determined by 
the branch at 170 Tm due to its comparably long half-life 
(U/2(/3 - ) = 128.59 d, almost independent of temperature). 
170 Yb is efficiently bypassed (/„ > 0.5) by the s-process flow 
for neutron densities larger than about 1.3x10 s cm” 3 when 
neutron captures prevail over (3 decays. Fig. IB13I shows the 
temporal evolution of the neutron density and of the iso¬ 
topic abundances of 170 ’ 171 Yb, 169 ’ 170 Tm and i6M69,i70 Er 
during the 15 th He shell flash in a 1.5 Mq model and half¬ 
solar metallicity. 170 Yb is strongly depleted during the peak 
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Table 4. Old solar main component isoto pic percentage contri b ution s from 94 Zr to 98 Mo (column 2; computed with our estimated 
95 Zr MACS, close to that recom mended bvlToukan &; Kappeled d 199(f) ; ’’OLD”), compared with our previous calculations (column 3; 
which include the 95 Zr MACS bv lLugaro et alj|2014l : ” L14”) and with s predictions obtained by assuming the 95 Zr MACS by KADoNiS 
(column 4; ”KAD”). In column 5 we list the updated main component obtained with the recommended 22 Ne(a, n) 25 Mg rate (Fig.JSjl- In 
column 6 we provide the results of Test A in Section 14.II in which we multiplied by a factor of four the recommended 22 Ne(ct, n) 25 Mg 
rate. Only variations larger than 5% are given in brackets. 
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Figure 12. The same as Fig. [8] but for the .s path close to the s- 
only isotope 179 Yb (red rectangle). The branch at 170 Tm (marked 
by the circle) is activated mainly during TPs. 
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Figure 13. The same as Fig. [9] but for the isotopic abundances 
of 170.171 Yb , 169,170 Tm and 168,169,170 Er 


neutron density, when both 169 Er and 1,0 Tm branches are 
activated and the s path proceeds mainly via 171 Tm (ii /2 
= 1.92 yr), thus bypassing m Yb as well. Most of 170 Yb 
is restored as the neutron density decreases (~74% of the 
170 Yb abundance produced during the previous 13 C pocket 
is reestablished). N ote that the 170 Tm /3 + -d ecay channel to 
170 Er is negligible llTakahashi fe Yokoilll987l ). 

The main component reproduces 99% of solar 170 Yb. 
About 10% of solar 170 Yb is missing with increasing the 
recommended 22 Ne(a, n) 25 Mg rate by a factor of four. The 
MACS of 170 Yb is very well determined with less than 1% 
uncertainty (768 ± 7 at 30 keV; KADoNiS). The major un¬ 
certainty that affects the abundance of 170 Yb is related to 
its SEF (~1.08 at kT = 23 keV; KADoNiS). 

Other uncertainties derive from the theoretical MACS of 
170 Tm (1870 ± 330 at 30 keV; 17.6%), which produces up 
to 6% variation on the 170 Yb s- predict io ns, an d from the 
170 Tm f }~-decay rate, for which iGorielvl (119991) estimated 
up to a factor of 1.7 uncertainty at Tg = 3, affecting the s 
contribution to 1,0 Yb by ~4%. 


5.2.3 Additional branches sensitive to the neutron density 

Besides 96 Mo and 170 Yb, we underline 142 Nd, 186 Os and 
192 Pt among the isotopes mainly sensitive to the 22 Ne(a, 
n) 25 Mg neutron burst. 

Similar to the branch at 95 Zr, the /3-decay rates of 85 I<r 
and 86 Rb are constant at stellar temperature, and the s 
contributions to 86,8 ' Sr are mainly sensitive to the neutron 
density reached during the 22 Ne(a, n) 25 Mg irradiation. The 
effect of the above branchings on the abundances of 86,87 Sr 
as well as on the neutron-magic nuclei 86 Kr and 87 Rb during 
the 13 C-pocket and TP phases are discussed in Section B 
(Supporting Information). 

The solar abundance of the neutron-magic isotope 
142 Nd is affected by the branches at 141 Ce and 142 Pr (par¬ 
tially open during TP), and decreases by ~15% with in¬ 
creasing the recommended 22 Ne(a, n) 25 Mg rate by a factor 
of four. 

While the solar abundances of 96 Mo, 170 Yb and 142 Nd are 
plausibly represented by the main component, the interpre¬ 
tation of solar 186 Os and 192 Pt is more prob lem atic. 

The studi e s car ried out by iMosconi et al.l {2010) and 
iFuiii et ahl (l20ld ) have substantially increased the accuracy 




































Main s-process Branchings and Neutron Source Uncertainties 17 


of the 186 Os MACS (4.1%). The s abundance of 186 Os is 
mainly influenced by the branch at 185 W, while smaller 
variations are produced by the branch at 186 Re (/3 decay 
rather constant at stellar temperatures). Although the half- 
life of 185 W is slightly sensitive to temperature and elec- 
tr on density (C = 75.10 d decreases to 50 d at T$ = 
3; ITakahashi fe Yokoil 1 19871 ) the s abundance of 186 0s is 
dominated by the uncertainty of the theoretical MACS of 
185 W. The present 185 W(n, 7 ) cross section overestimates 
the solar 186 Os abundance by 20-30%, hardly to be rec- 
onciled with the 8 % uncertainty of the solar abundance 
llLodders et al.l l2009ll . The accuracy of the 185 W MACS 
recommended by KADo NiS (9%) could be largely under¬ 
estimated fsee iRauscheij 120141 1. be ing based on an aver¬ 
age among inverse ( 7 , n) reactions (ISonnabend et ahl l2003l : 

I Mohr et al.l 120041 : IShizuma et alJ 120051 ). In order to repro¬ 
duce the solar abundance with our AGB models, we adopt 
a MACS for 185 W which is about 80%) higher than recom¬ 
mended in KADoNiS. 

Present calculations underestimate the solar 192 Pt abun¬ 
dance by ~20%. The abundance of 192 Pt is mainly deter¬ 
mined by the branch at 192 Ir. Similar to 186 0s, the s contri¬ 
bution of 192 Pt is also dominated by the uncertainty of the 
neutron capture channel rather than by the /3-decay branch 
(slight reduction of the half-life of 192 Ir from 77.54 d to 55.33 
d at Tg = 3; ITakahashi fc Yokoil 1987 ). The MAC S mea- 


surements on the Pt isotopes by iKoehler fc Gubei] (|2013i ) 
included also a new calculation of the values for 192 Ir. The 
discrepancy is reduced by accounting of a 20 % uncertainty 
of the theoretical 192 Ir MACS, which increases the 192 Pt s- 
prediction to ~85%. However, a 15% missing contribution 
scarcely agrees with the much improved accuracy of the 
192 Pt MACS measured by IKoehler fc Gub ei] d2013|) (4%), 
and with an 8 % solar uncertainty ( Lodder s et al. 2009h . A 
more detailed theoretical analysis of the branch at 192 Ir 
would help to improve the 192 Pt s-prediction. 

In Section B (Supporting Information), the uncertainties af¬ 
fecting the s predictions of 142 Nd, 186 Os and 192 Pt are anal¬ 
ysed in detail. 

While the s contribution of 192 Pt is marginally influ¬ 
enced by the 22 Ne(a, n) 25 Mg uncertainty, the present s 
contributions to solar 96 Mo and 170 Yb better agree with 
our newly evaluated 22 Ne(a, n) 25 Mg rate. However, the 
MACS of 170 Yb is affected by a non-negligible SEF cor¬ 
rection (1.15 at 30 keV; KADoNiS) although its neutron 
capt ure cross section is w ell determined experimentally (see 
also iRauscher et alJ 1201 il l . 


5.3 The s-only isotopes strongly sensitive to 

stellar temperature (and/or electron density) 
and neutron density 

The /3-decay rates of several unstable isotopes exhibit a 
strong dependence on te mperature and electron density 
(Takahashi & Yokoi Il987l l. As anticipated in Section 0 
the uncertainty affecting the stellar /3-decay rate of a few 
branches may produce wide variations of the s predictions. 

The branches at issue are 134 Cs (affecting the s contri¬ 
bution to solar 134 Ba), 151 Sm and 154 Eu (which influence 
the 152,154 Gd abundances), 176 Lu (making it an s-process 
thermometer rather than a cosmic clock), and 204 T1 (with 
consequences for the s contribution of 204 Pb). 
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Figure 14. The same as Fig. [8] but for the s path close to the 
s-only isotopes 134 Ba and 136 Ba (red squares). The main branch 
point occurs at 134 Cs. The electron-capture (EC) channel of 134 Cs 
is negligible. The long-lived 135 Cs can be considered almost sta¬ 
ble during the s process, and its radiogenic contribution to 135 Ba 
occurs after the TP-AGB phase (thick dashed arrow); the abun¬ 
dance of 137 Cs stored during TPs /3~ decays to 137 Ba during 
interpulse periods (thin dashed arrow). The neutron-magic nuclei 
at N = 82 are indicated by a blue box. 


5.3.1 The s-only pair 134,138 Ba (the branch at 134 Cs) 

The production of the two s-only isotopes 134,136 Ba occurs 
via both neutron irradiations, in the 13 C pocket and during 
TPs. 

Most of solar 134 Ba is produced in the 13 C pocket, when the 
s path proceeds via 132 Xe and 133 Cs directly to 134 Ba, be¬ 
cause the low neutron density does not allow neutron cap¬ 
tures on 133 Xe and 134 Cs (Fig. fl4t . In this way the solar 
134 Ba/ 136 Ba ratio would be overestimated by about a fac¬ 
tor of two. The additional 22 Ne(a, n) 2o Mg neutron burst 
partially activated during TP is essential to regulate the 
I34,i36g a p rec [ictions in order to reproduce the solar abun¬ 
dances. The main component is expected to reproduce the 
s contribution of the s-only pair 134 Ba and 136 Ba so that 
134 Ba: 136 Ba = 1:1 within solar and nuclear uncertainties. 

The dominant uncertainty affecting the 134 Ba/ 136 Ba ra¬ 
tio derives from the /3“-decay rate of 134 Cs: the terrestrial 
half-life of 134 Cs (ti /2 = 2.07 yr) is strongl y reduced under 
stellar conditions (ITakahashi fc Yokoi 19871) : it decreases by 
a factor of three at Ts = 1 (0.67 yr) and by two orders of 
magnitude at Ts = 3 (3.8 d). Note that A _ ( 134 Cs) does not 
change with electron (or mass) density. 

In previous AGB models (see filled circles in Fig. 0, we 
have adopted a constant /3 ~-decay rate for 134 Cs, which 
was based on a geometr ic average of the rates given by 
ITakahashi fc Yo koi (1987) at different temperatures over the 
con vective pulse (A~( 13 4 Cs) = 1.6xl0 ~ 7 s _1 , ti / 2 = 50 d; 
see iBisterzo et ahll201ll : iLiu et ahl 1 20141) . This value repro¬ 
duces the solar 134 Ba/ 136 Ba ratio, and the predicted s con¬ 
tributions exceed the solar abundances of 134 Ba and 136 Ba 
by only 6 and 9%, respectively, well within the uncertain- 
ties quote d for the solar abundances ( 6 % for Ba, 5% for Sm; 
iLodders et ali l2009l l and for the MACS values of barium 
(3.2% at 30 keV; KADoNiS). 

The strong temperature dependence of the half-life of 
134 Cs is now considered by an improved treatment of the /3- 
decay rate as a function of the temperature gradient within 
the TPs. The maximum temperature reached by a given 
AGB model at the bottom of the convective intershell rises 
from pulse to pulse. The convective TP itself is characterised 
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by a very large gradient in temperature and density, which 
decreases from Tg ~ 3 and p = 2xl0 4 g/cm 3 at the bottom 
of the advanced TP to Ts ~ 0.1 and p = 10 g/cm 3 in the 
top layers. Under these conditions, the 134 Cs half-life varies 
from 2 yr in the top layers of the convective TP (where no 
neutrons are available) down to 3.8 d in the region of highest 
neutron density at the bottom. Accordingly, the production 
of 134 Cs depends strongly on the locus inside the TP: for Tg 
~ 3 one finds ti/% = 3.8 d, and at N n ~ 4xl0 9 cm -3 only 
~25% of the s flow bypasses 134 Ba. Slightly further outward, 
at Ts ~ 2.7, neutrons are still produced via 22 Ne(a, n) 25 Mg, 
but the half-life of 134 Cs has increased to 50 d: under these 
conditions, neutron captures dominate over /3 decays start¬ 
ing from N n ~ 10 9 cm -3 (/„ ~ 0.5), and 134 Ba is more 
efficiently bypassed. 

We have interpolate d (with a cubic-spline me thod) the (3~- 
decay rates given bv lTakahashi fc Yokoil ll 19871 ) as a function 
of the stellar temperature in the range from Tg = 0.5 to 5. 
The adopted rate results very close to that recommended 
by NETGElsQ. 134 Cs is freshly produced in the thin bottom 
layers of the TP, but even in this hot environment 134 Cs has 
half-life of a few days, much larger than the turnover time, 
so that it is quickly brought to the outer and cooler layers 
of the TPs. Here, the 134 Cs half-life increases up to 2 yr, but 
the lack of neutrons does not allow further neutron captures 
and it decays into 134 Ba. Because a small amount of 134 Cs 
is converted to 135 Cs during the peak neutron density at the 
bottom of the TPs, 134 Ba is temporarily reduced (by ~30%). 
The initial amount of 134 Ba is almost fully reestablished as 
soon as the neutron density decreases. Consequently, almost 
the entire s flow is directed towards 134 Ba, resulting in a 
large overestimation of solar 134 Ba (+30%). 

This problem with the branching at 134 Cs can be somewhat 
relaxed if one consi ders the uncert ainty in the temperature- 
dependent half-life. iGorie iCiil) estimated a possible de¬ 
crease of the decay rate by a factor of ~3 at Tg = 2 - 3. As¬ 
suming twice that change as 2cr-uncertainty would reduce 
the overestimation of the solar 134 Ba (see filled diamond 
in Fig. 0. In this case, 134 Ba is largely depleted at the be¬ 
ginning of the advaced TPs (see, e.g., Fig. 1151 which shows 
the temporal evolution of the isotopic abundances during 
the 15 th TP in a 1.5 Mq model at half-solar metallicity). 
As the neutron density decreases, the /3 ~-decay channel is 
favoured and the 134 Ba abundance is increasing again ( N n 
< 2x10 s cm~ 3 ; f n ^ 0.1). About 76% of 134 Ba produced 
during the previous 13 C pocket is restored at the tail of the 
neutron density. The resulting s contributions to 134 Ba and 
136 Ba exceed the solar abundances by +12% and +9%, re- 
spectiv ely. This may sugges t that the stellar 134 Cs /3 _ decay 
rate by iTakahashi fe Yokoil (1 19871 4 has been overestimated. 
Further investigations on this topic are advised. 

An additional important effect concerning the s predictions 
for the Ba isotopes derives from the 22 Ne(a, n) 25 Mg rate. 
The above s contributions to 134 Ba and 136 Ba may be re¬ 
duced to +4% and +2% with increasing the recommended 
22 Ne(a, n) 25 Mg rate by a factor of four (Test A). 

A marginal impact on the 134:136 Ba s-predictions de¬ 
rives from the branches at 133 Xe and 135,136 Cs. The short 
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Figure 15. The same as Fig. [9] but for 134 , 135,136,137 
133 , 134 X e and 133,134,135,136 Cg isotopic abundances. 


half-life of 133 Xe (ti /2 = 5.24 d) does not allow efficient neu¬ 
tron captures (less than 2% of the s flow feeds 134 Xe). 
Although the terrestrial half-life of the long-lived 135 Cs 
(ti /2 = 2.3 Myr) is reduced by four orders of magni¬ 
tude at T s = 3 (e.g., U /i = 267 yr at p = 3000 g/cm 3 ; 
ITakahashi fc Yokoil 119871 1. it is practically stable compared 
to the 6 year timescale of the 22 Ne(a, n) 25 Mg irradiation 
with a TP. Therefore, once the s flow feeds 135 Cs, it pro¬ 
ceeds towards 136 Cs, and 135 Ba may be easily bypassed (see 
Fig. 11511 . The abundance of the long-lived 135 Cs stored dur¬ 
ing TPs /3~ decays into 135 Ba at the end of the TP-AGB 
phase. 

The s flow at 136 Cs (U /2 = 13.04 d) continues mainly via /3~ 
decay to 136 Ba, where the two s branches formed at 134 Cs 
join again. Accordingly, the s contribution to 136 Ba keeps in¬ 
creasing during the entire 22 Ne(a, n) 25 Mg irradiation, thus 
enhancing the abundance from the previous 13 C pocket by 
+ 10%. 

The 134 Cs MACS h as been estimated semi-empirically 
by IPatronis et al.l (l2004h , who provide a rather small un¬ 
certainty of 9% (KADoNiS), resulting in a ~2.5% varia¬ 
tion of the 134 Ba abundance. The uncertainties of the the¬ 
oretical MACS value s of 133 Xe (KADoNiS) and 136:137 Cs 
llPatronis et al.l 120041 4 are presently estimated between 30 
and 50%, but have almost no impact (<2%) on 134,136 Ba. 


5.3.2 Additional branches strongly sensitive to stellar 

temperature (and/or electron density) and neutron 
density 

A few branches behave as thermometers of the s process 
in AGB stars because their /1-decay rates deeply rise (or 
drop) with the large temperature and ele ctro n de nsity g radi¬ 
ents that characterise a convective TP llTakahashi fe Yokoil 
Il987l4 . In a few cases, this /3-decay behaviour is responsible 
for crucial uncertainties. 

Similar to the way 134 Cs affects the s contribution to 
134 Ba, the decay rates of 151 Sm and 154 Eu depend strongly 
on temperature and electron density during TPs with cor¬ 
responding consequences for the s abundances of 152 Gd and 
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154 Gd, respectively. 

The s predictions to 152,154 Gd have been improved by in¬ 
cluding in AGB models an appropriate treatment of the 
nearby /3-decay rates over the full convective TPs (see filled 
diamonds in Fig. [5j: about 85% and 92% of solar 15 2 Gd and 
154 Gd are produced by AGB stars. According to iGorielvl 
(119991) both rates may vary up to a factor of three at T$ 
= 3. This would produce an extreme impact on the 152 Gd 
s-prediction (up to a factor of 2), and up to 10-15% vari¬ 
ations of the 154 Gd s-contribution. The dominant effect 
of the 151 Sm half-life prevents the s contribution to so¬ 
lar 152 Gd from being accurately assessed. Additional uncer- 
tainties may derive f rom the 151 Sm and 154 Eu MACS (see 
IRauscher et al .11201 ill . 

Certainly, the s predictions of the above isotopes are largely 
influenced by the uncertainties of the 22 Ne(a, n) 25 Mg rate. 
A detailed analysis of these aspects is given in Section B 
(Supporting Information). 

Two other isotopes are worth mentioning, where the 
AGB contributions depend strongly on temperature: the 
long-lived 176 Lu (more properly classified as an s-process 
thermometer rather than a cosmic clock), and 204 T1 (which 
regulates the s contribution of solar 204 Pb). 

At present, the discrepancy between experimental data and 
the astrophysical tre atment of the branch at 176 Lu re- 
mains to be solved dHeil et al.l l2008d : iMohr et al.l 120091 ; 
ICristallo et al.l l2010h . The solution may be found in the 
nuclear coupli ng scheme between thermally populated lev- 


els in 176 Lu 

Gintautas et al. 20091: Dracoulis et al. 2010i; 

Gosselin et al. 

2010|). Besides the branch at 176 Lu. the 


22 Ne(a, n) 25 Mg rate and the SEF estimate of the 176 Hf 
MACS are directly influencing the predicted 176 Lu/ 176 Hf 


ratio. 

Although the MACS of 204 Pb is accurately known (3%), the 
present 204 Pb s prediction is jeopardised by the uncertainty 
of the theoretical 204 T1 neutron capture cross section and by 
the evaluation of its stellar /3 _ -decay rate (~10% variations). 
An ad ditional 10% uncertainty from the solar Pb abundance 
(iLodders et al.ll2009h has to be considered as well. 

The impact of the present uncertainties is discussed in 
Section B (Supporting Information). 
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Figure 16. The same as Fig. [8] but for the s path close to 
lso Ta m . During the 13 C(a, n) le O neutron irradiation 179 Hf is 
stable and the s path follows the thick red lines; during 22 Ne(a, 
n) 25 Mg neutron burst (thin red lines) the isomeric state of 179 Hf 
becomes unstable (with h alf-life strongly temperature dependent; 
llTakahashi fc Yokoill 1 08?) )). allowing the branch of the s flow to¬ 
wards 179 Ta, and the stable 180 Ta m . The long-lived 182 Hf (fi /2 
= 8.9 Myr) can be considered almost stable during the s process, 
and its radiogenic contribution to 182 Ta occurs after the TP-AGB 
phase (red dashed-thick line). 
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Figure 17. The same as Fig. [9] but for the isotopic abundances 
Of 180 Ta m , l 79 ,181 Tai 180 W and 178,179,180 


5.4 The s-only isotopes mainly sensitive to stellar 
temperature (and electron density) 

The branches at issue are 179 Hf ( responsible for the s 
contribution to 180 Ta m ; see, e.g., iKappeler et al.l l200-ll : 
IWisshak et al.l 120041 '). 164 Dy (which produced about half 
of solar amount of the proton-rich isotope 164 Er; 
Ijaag fc Kappeleilll996lt. a nd 128 I (which regulates the s pro¬ 
duction to 128 Xe, Reifarth et al j|2004lh 

Because of their low sensitivity to the neutron density, the 
branches at 1S0 Ta m and 128 I have provided information on 
the convective mixing time scale. However, the large un¬ 
certainties of the /3-decay rates of the respective branch 
point nuclei have led to problems in the interpretation of 
the branchings towards 164 Er and 180 Ta m . 


5.4-1 The s contribution to 180 Ta m (the branches at 
179 Hf, 179 Ta and 1S0 Hf) 

In the solar system, 180 Ta occurs only in its stable isomer 
(fi /2 > 10 15 yr). 180 Ta m is the rarest stable isomer in Na¬ 
ture and constitutes 0.012% of solar Ta. The ground state 
of 180 Ta is unstable and decays with a half-life of 8.15 li by 
(3~ and /3 + emission to 180 W and 180 Hf, respectively. 

The origin of 180 Ta m has been a challenge for years. Be¬ 
ing proton-rich, the expected contributions from p(y)~ a nd 
t^p-processes have been largely investigated (e .g, iRavet et al.l 
Il995l : IRauscher et alf2002t Heeer et alj|200a ) . As illustrated 
in Fig. [16] the s flow is bypassing 180 Ta m because 179 ’ 180 Hf 
are stable during the 13 C(a, n) le O irradiation. At the higher 
temperatures during TPs, however, higher lying nuclear 
states in 179 Hf and 180 Hf are thermally populated, opening 
two branches towards 180 Ta m : (i) via /3~ decay of the ther- 
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mally populated 179 Hf state at 214 keV, which /3 _ decays 
to 179 Ta and produces 180 Ta m via neutron captures, and (ii) 
via the 179 Hf partial neutron capture cross section towards 
the weakly populated 180 Hf isomer, which /3“ decays quickly 
to 180 Ta m . 


Several studies have been dedicate d to the s- 

proce s s nucleosynthe sis _ of 180 Ta m (e.g., iNemeth^ffiJ 

19921: iKaDDeler et al.l 12004 IWisshak et al.l l20()4 |2006b; 


Briefly, about half of the s contribution to 180 Ta m is pro¬ 
duced via the branch at 179 Hf m . At Tg = 3, the 179 Hf half- 


life becomes sensitive to te mperature and electron d ensity 
(e.g., t 1 / 2 from 85 to 25 vr: ITakahashi fe Yokoi|[l987l ). feed¬ 
ing 1,9 Ta in small amounts. Neutron captures on 179 Ta are 
dominating over the small /3 + decay back to 179 Hf (the ter¬ 
restrial 179 Ta half-life increases by a factor of two at TP 
temperature), producing mainly 180 Ta g , which /3 + - and /3~ 
decays into 1S0 Hf and 180 W, respectively. This channel yields 
about 5% of the solar p-rich 180 W. A minor amount of 179 Ta 
(IR = 0.04) feeds the stable isomer 180 Ta m , contributing 44% 
of its solar abundance. At the end of the TP, the residual 
179 Ta abundance /3 + decays to its stable Hf isobar IFig. llTI) . 
An additional 39% s contribution to 180 Ta m derives from 
the /3~ decay of 180 HP 1 . The production of the 180 HP 1 iso¬ 
meric state is weak: only ~1.3% of the total 179 Hf MACS 
feeds 180 HP 1 (at 30 keV cr[ 179 Hf(n, 7 ) 180 Hf] = 922 ± 8 
mbarn; cr[ 179 Hf(n, 7 ) 180 HP 1 ] = 11.4 ± 0.6 mbarn; KADo- 
NiS). 180 Hf m has short half-life (ti /2 = 5.5 h), and decays 
directly to 180 Ta m . This contribution is independent of tem¬ 
perature and is completely determined by the partial cross 
section to 180 HP 1 . 


The final s abundance of 180 Ta m is larg ely sensitive 
to te mperature and electron density llTakahashi fe Yokoil 
Il987t) . and needs to be properly evaluated by account¬ 
ing of the temperature and densi ty gradie nts over the 
convective He flashes. IBelic et al.l (Il999l . 120021 1 carried 
out a photoactivation experiment of 180 Ta m to study the 
probability for connecting isomer and ground state in 180 Ta 
via thermally induced transitions to higher lying mediating 
states. Whereas the direct internal decay of the isomer 
to the ground state is highly forbidden by selection rules, 
thermal excitations of such mediating states are allowed 
under stellar conditions. They found that the two states 
are fully thermalised at Tg = 3: thus, 180 Ta m should be 
destroyed in the bottom layers of the advanced TPs, where 
such temperatures are reached. Instead, the fast convective 
mixing (of the order of a few hours) occurring during TP 
prevents the destru ction of 180 Ta m (see Fig. flTl) . 


Starting from Wisshak et al.l (1200 ll . l2004l l and 
iKappeler et al.l (120041 ). our AGB models account for 
the thermally induced destruction of 1 80 Ta m obtai ned by 
the photoactivation measurement by IBelic et al.l d 19991 . 
Hi), by following the strong half-life variation with 
temperature and density gradients together with convective 
mixing at each TP. This method firstly pointed out that 
180 Ta m receives a dominant contribution from the main 
component, incr e asing the previous estimate from ^49% 
jAHandrni_et_ak _ 19991 4 to ~80% dKappeler et alJ l2004 
1 Wisshak et alj|2004h . In addition, it provides information 
on the convective turnover time during He-shcll flashes, and 
underlines that temperature gradient and neutron freeze- 


out effects are not sufficient to analyse the abundances of 
such peculiar isotopes correctly. 

The main component reproduces about 83% of solar 
180 Ta m (see filled circle in Fig. EJ. This prediction is af¬ 
fected by a number of uncertainties, however. Therefore, the 
fraction of ~20%, which have to be contributed by other 
sources, could well be significantly larger. 

In order to improve the estimated s contribution to 180 Ta m , 
we have included the same treatment to the 179 Ta /3 + -decay 
and 179 Hf j3~- decay rates over the full convective zone. The 
s contribution to 780 Ta m is reduced from ~83% to ~68% 
(see filled diamond in Fig. El). This value remains strongly 
se nsitive to the uncertain /3-decay rate of 179 Hf estimated 
by IGorielvl (jl999|) who suggested that the 179 Hf /3 -decay 
rate may decrease by a factor of three, resulting in a strong 
impact on the contribution to solar 180 Ta m (—30%). The 
smaller 1,9 Ta /3 + -decay uncertainty estimated by IGorielvl 
(119991 ) does not affect 180 Ta m . 

As a consequence, AGB stars remain to be the major nu¬ 
cleosynthesis sites for 180 Ta m , but leave a larger fraction for 
abundance contributions by other processes. 

In light of this result, we suggest that an investigation of 
the 179 Hf half-life under stellar conditions would be impor¬ 
tant in order to assess the contribution to 180 Ta m from AGB 
stars more reliably. 

We examine additional uncertainties associated to the 
180 Ta m s prediction. Although the 180 Ta m MACS is rather 
well kno wn with 6.8% uncer tainty (1465 ± 100 mbarn at 
30 keV; IWisshak et alJ 120041 ), it is affected by a non neg¬ 
ligible stellar enhancement factor (SEF = 0.87 at 30 keV; 
KADoNiS). Moreover, the 179 Ta MACS is only estimated 
theoretically with ~30% uncertainty (1334 ± 422 at 30 
keV; KADoNiS), producing ~10% variations of the 180 Ta m 
s abundance. 

Finally, the 180 Ta m s-prediction decreases by ~12% with in¬ 
creasing the 22 Ne(a, n) 25 Mg rate by a factor of four (see 
Table Al, Supporting Information). 


5-4-2 Updated s contributions to 182 W and the short-lived 

182 H f 

The s and r process contributions to solar 182 W and 182 Hf 
are mainly regulated by the branch at 181 Hf (see Fig. 1161) . 
Present AGB predictions marginally produce the short-lived 
182 Hf (fi /2 = 8.9 Myr), because the terrestrial half-life of 
lsl Hf is strongly redu ced during TPs (fro m t \/2 = 42.39 
d to 1.26 d at Tg = 3; ITakahashi fe Yokoilll987l ) and the s 
path mainly feeds 182 W. We estimate that 65% of solar 182 W 
is synthesised by the main component (see Table A.l, Ap¬ 
pendix A, Supporting Information), which corresponds to a 
residual 35% r contribution , in agre emen t with previous cal¬ 
culati ons bv I Wisshak et ahl ll2006bl ) and lYockenhuber et alJ 
ll2007l) . This value includes the radiogenic contribution of 
182 Hf at the end of the TP-AGB phase. Indeed, no /3-decay 
occurs during TPs because the decrease of the 182 Hf half-life 
is not sufficient to compete with neutron captures (ti /2 = 
22 yr at Tg = 3; Takahashi fe Yokoilll987l ). 
iLugaro et al.1 (l2014al ) have recently demonstrated that the 
present AGB contributions to 182 Hf and 182 W have so far 
be en underestimated . Start ing from the experimental work 
bv iBondarenko et al.1 (I2002I) . who did not find any evidence 
for the existence of the 181 Hf levels responsible for the strong 
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enhancem ent of the /3-d ecay ra te (at 68 keV, 170 keV, and 
298 keV), iLuearo et al.l (l2014al l argued that the terrestrial 
/3-decay rate of 181 Hf remains rather unchanged in stellar 
condition. This favours the neutron capture channel towards 
182 Hf, which increases the abundance of 182 W after the TP- 
AGB phase. By assuming a constant 181 Hf /3-decay rate, 
we find a negligible r contribution to solar 182 W, which 
would be ful ly produced by AGBs. On the other hand, 
iLugaro et al. (l2014a f estimated a new 181 Hf /3-decay rate 
by removing the 68 keV level (see their Fig. S2; top panel, 
blue line). By including this value in our calculations, the s 
contribution of 182 W increases from 65% to 80%, reducing 
the r component by a factor of two. The lower limit given 
by the authors (obtained by removing the three levels; top 
panel, red line) further increases the solar s abundance of 
182 W to 86%. The 182 W s prediction shows —3% variations 
within the 22 Ne(a, n) 25 Mg rate uncertainties. 

This result is of great re levance to the chronome try of the 
early solar system (see IWasserburg et ahl 120061 for a re- 
view on shor t-lived isotopes), or in presolar stardust grains 
llAvila et al.ll2012h . 

In addition, it provides another key to explain the ap¬ 
parent anomaly in the r-process residuals at A = 182, to¬ 
gether with the uncertainty ass ociated to the neut r on ca p- 
ture cross sections advanced by I Vockenhuber et al] J2007I '). 

5-4-3 Additional branches mainly sensitive to stellar 
temperature (and/or electron density) 

Similarly to 180 Ta m , also the p-rich isotope 164 Er receives a 
dominant contribution from the s process owing to the ac¬ 
tivation of the branch at 163 Dy (which becomes unstable in 
stellar environment). 

About 80% of solar 164 Er is produced during TPs ( filled 
circle in Fig. GO). This value rises to 87% by including an 
improved treatment of the /3-decay rates of nearby unsta¬ 
ble isotopes over the convective TPs (see filled diamond in 
Fig- E}. Unfortunately, the s contribution to 164 Er can not 
be accurately assessed, being largely sensitive to the com¬ 
petition between the /3~ and /3 + -decay rates of 164 Ho, and 
in particular to the uncertainty of the 164 Ho /3~- decay rate 
estimated under stellar conditions (10-20%). We discuss in 
Section B (Supporting Information) the main uncertainties 
affe cting the s cont r ibutio n to 164 Er. 

iReifarth et al.l (120041 ') provided a test for the convec¬ 
tive mixing time scale during TPs, by examining the branch 
at 128 I, which regulates the production of the s-only pair 
i28,i3C>x e The s predictions of 128,130 Xe exhibit a weak de¬ 
pendence on neutron density owing to the short half-lives of 
127 Te (t- 1/2 = 9 h) and 128 I (t\/ 2 = 25 min), and the com¬ 
petition between 128 I /3-decay and electron captures may 
constrain the convective mixing timescale during the He 
shell flash. While the /3 _ -decay channel of 128 I shows only a 
weak dependence on temperature, the electron capture rate 
is strongly temperature-dependent: from t 1 / 2 (i?C') ~ 6 h at 
Ts = 0.5 up to ~8 d at Ts = 3. By assuming a sufficiently 
short turnover mixing timescale, the 128 I produced via neu¬ 
tron captures in the hot bottom layers of the TP is promptly 
brought to the cooler external layers of the convective zone, 
allowing a partial activation of the EC channel with the re¬ 
sult that ~5 to 6% of the s-process flow is bypassing 128 Xe. 
The corresponding s predictions for 128 Xe and 128 Xe are in 


agreement with the 128,130 X e ratio observed in SiC grains 
and confirm the results of IReifarth et all d2004l ) . (see Sec¬ 
tion B, Supporting Information). The 128,130 Xe s-predictions 
are rather well determined, because the branching is com¬ 
pletely regulated by the decay of 128 1, so that the influence 
of uncertain theoretical MACS values of nearby unstable 
isotopes is negligible at stellar temperatures. 


6 SUMMARY AND CONCLUSIONS 


Despite a Galactic chemical evolution model would provide 
a more realistic description of the abundances observed in 
the solar system, the main component is still a useful tool 
to investigate the s-process nucleosynthesis of nuclei with 
atomic mass between 90 and 204. 

We have studied the major uncertainties affecting the nu¬ 
clear network. The analysis has been carried out with the 
most recent neutron capture cross sections and with updated 
solar abundance data. 


We have examined the impact of the present uncertain¬ 
ties of the two neutron sources operating in AGB stars, the 
13 C(a, n) 16 0 and 22 Ne(a, n) 25 Mg reactions, focusing on 
the s-only isotopes sensitive to the most important branch 
points of the main component. 

The overall s distribution of isotopes heavier than A 
~ 90 shows negligible variations (up to ~1%) by changing 
the 13 C(a, n) le O rate by about factor of two. Only the two 
neutron-magic nuclei 86 Kr and 87 Rb are influenced by the 
13 C(q, n) le O neutron irradiation, because of the marginal 
activation of the branch at 85 Kr. Note, however, that in the 
AGB models we adopt to reproduce the solar main com¬ 
ponent (M = 1.5 and 3 M@ at half solar metallicity), the 
13 C abundance in the pocket is exhausted radiatively dur¬ 
ing the interpulse period (only a negligible amount of 13 C 
is engulfed in the subsequent TP). Thus, our AGB mod¬ 
els do not experience a partial convective burning of 13 C 
during the first TPs, which could affect the production of 
a few neutr on-rich isotopes (e.g. , 86 Kr and 87 Rb, or 96 Zr 
as well, see iCristallo et ahl 120061) . The new measurement 
bv lLa Cognata et al.l ( 20131 ') suggests that the 13 C(a, n) 16 0 
rate adopted in our models should be increased by 40%, 
rather than decreased. In this case, the convective 13 C burn¬ 
ing during TPs is significantly reduced. 

However, beside the un deniable progress made b y means of 
indirect measurements (lLa Cognata et ahl l2013h , the exis¬ 
tence of a sub-threshold state makes the evaluation of the 
rate at astrophysical energies still uncertain. A further in¬ 
crease of this rate would likely have marginal consequences 
on the main component, a substantial reduction may in¬ 
crease the amount of 13 C engulfed in the convective pulse 
and burned at relatively high temperature. 

The prese nt unc ertainty of 2 2 Ne(a , n) 25 Mg (e.g.. 


Kappeler et al.l 1 1994 I Jaeger et al.1 1200 ll : iLongland et all 


2012!) mainly influences the isotopes close to and within the 


branchings of the s path. 

We have provided new evaluated values of the 22 Ne(a, 
n) 25 Mg and 22 Ne(a, 7 ) 26 Mg rates that account for all known 
and potential resonances as well as all nuclear data avail¬ 
able. In the temperature range of AGB stars (Ts = 2.5-3), 
the recommended 22 Ne(a, n) 25 Mg value is about a factor of 
two lower than that adopted in our models so far, while the 
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recommended 22 Ne(a, 7 ) 26 Mg rate is essentially unchanged. 
The recommended rates are in agreement with in the errors 
with the recent values pres ented in literature ll.Iaeger et al .1 
l200ll : lLongland et al.ll 2012 h . 

However, the 22 Ne(a, n) 25 Mg and 22 Ne(a, 7 ) 26 Mg rates are 
mainly based on the knowledge of the 832 keV resonance. As 
discussed for the 13 C(a, n) 16 0 neutron source, the presence 
of unknown low-energy states, which have been identified in 
several indirect experiments, makes the 22 Ne(a, n) 25 Mg and 
22 Ne(a, 7 ) 26 Mg recommended values uncertain at stellar en¬ 
ergy. A direct determination of these reaction rates at tem¬ 
peratures below the present experimental limits will shed 
light on the actual efficiency of the AGB neutron sources. 
Because the available low-energy measurements were mainly 
limited by a significant neutron background, deep under¬ 
ground laboratories are the most promising places to plan 
future experimental investigations on the two AGB neutron 
sources (LUNA, DIANA). 

For this reason we have analysed the solar main component 
within a more cautious range of uncertainty, corresponding 
to variations by a factor of four starting from our recom¬ 
mended 22 Ne(a, n) 25 Mg rate. 

Major variations are shown by the s-only nuclides close to 
the branchings, which are most sensitive to neutron density 
as 96 Mo, 142 Nd, and 170 Yb (due to the branch points at 95 Zr, 
141 Ce, 170 Tm), 134 Ba and 152 Gd (which also depend on the 
134 Cs and 151 Sm half-lives, respectively, both strongly re¬ 
duced in stellar environments), and 176 Lu (affe cted b y the 
contribution of isomeric mediating states, see [Mohr ct al. 

l200dl . 

In low mass AGB stars (which are the major contributors to 
the s olar s abundances: iTravaglio e t al. 2004; Bis terzo et al] 
l2014l l. the 22 Ne(ct, 7 ) 26 Mg reaction is not efficiently acti¬ 
vated during TPs owing to the rather small temperature 
reached at the bottom of the advanced convective He flashes. 
Thus, the competition between 22 Ne(a, 7 ) 26 Mg and 22 Ne(a, 
n) 25 Mg reactions is marginal and does not influence the so¬ 
lar s distribution. 

A larger impact of the 22 Ne(a, 7 ) 26 Mg and 22 Ne(a, n) 25 Mg 
rates is expected from intermediate AGB stars, where the 
22 Ne(a, n) 25 Mg is the major neutron source, or from low- 
metallicity AGB models. The study of these models pro¬ 
vides information about the s -process nucleosynthesis in, 
e.g., globular clusters, dwarf galaxies, intrinsic or extrinsic 
peculiar stars showing s enhancement. To this purpose, we 
have analysed the rate uncertainties in a 3 M@ model at 
[Fe/H] = —1 and a half-solar metallicity 5 M@ model cho¬ 
sen as representative of an extended range of AGB stars. 


The status of the (n, 7 ) stellar cross sections has been 
significantly improved in the last decade, with accuracies 
of less than 5% for a number of isotopes in the mass re¬ 
gion 120 < A < 180. On the other hand, the accuracy of 
neutron capture and /3-decay rates of isotopes that act as 
important branching points of the s path plays a crucial 
role in the production of a few s-only nuclides. First, be¬ 
cause the MACS of unstable isotopes are barely accessible to 
direct measurements, only theoretical estimates (with large 
uncertainty) are available in most cases. Second, because the 
stellar /3-decay rates are poorly known for branch point nu¬ 
clei, in particular if they are extremel y sensitive to tempera - 
ture and electron density. Moreover, Rauscher et all ([2011) 


and iRauscherl (l2012ll highlighted that the (experimentally 
measured) ground state cross section may constitute only a 
minor fraction of the MACS. Thus, the theoretical uncer¬ 
tainties associated to th e MACS may be in a few cases un¬ 
derestimated (note that fRauscher et al.ll201ll provide upper 
limits of MACS uncertainties; specific theoretical investiga¬ 
tions carried out individually for each nucleus are strongly 
needed for a few isotope^). 

In this study, we have discussed the present major nuclear 
uncertainties that affect the s-only isotopes. We have dis¬ 
tinguished the s-only nuclei in different classes, according 
to the type of information that can be deduced from their 
abundances with respect to the physical conditions in AGB 
stars: unbranched nuclides (useful to constrain the s distri¬ 
bution), isotopes mainly affected by neutron density, and 
isotopes strongly sensitive to temperature and electron den¬ 
sity (which help to address the characteristics of stellar mod¬ 
els). For each class, the specific problems and suggestions 
for a possible improvement are given. We suggest that an 
investigation of the /3-decay rates as a function of the stellar 
environment would be important for 134 Cs, lsl Sm, 179 Hf, 
and 164 Ho in order to improve the s contributions of 134 Ba, 
152 Gd, 180 Ta m , and 164 Er. 

I 11 conclusion, we find that the solar main component 
may still reproduce the s-only isotopes within the present 
uncertainties. 
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SUPPORTING INFORMATION 

Additional Supporting Information may be found in the on¬ 
line version of this article: 

Appendix A. The results for the updated main com¬ 
ponent obtained in the analysis described in Section 2 are 
listed in Table Al for the isotopes with A > 70. We also 
report the effects of the tests discussed in Section 4.1. 

Appendix B. This Appendix completes the informa¬ 
tion about the most important branches of the main com¬ 
ponent discussed in the paper. Following the classification 
given in Section 5, we analyse here the additional branch¬ 
ings of each category, which are only briefly outlined in the 
text. 
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Appendix C. In this Appendix we have analysed the 
uncertainties of the neutron sources in a 3 M© model at 
[Fe/H] = — 1 and a half-solar metallicity 5 M© model chosen 
as representative of an extended range of AGB stars. 
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APPENDIX A: THE MAIN COMPONENT OF THE SOLAR SYSTEM 

The results for the updated main component shown in Fig. 5 are listed in Table A1 for the isotopes with A > 70. The 
recommended 22 Ne(a, n) 25 Mg and 22 Ne(a, 7 ) 26 Mg rates described in Section 4 have been adopted. 

We also report the effects of a conservative range of uncertainty associated to the 22 Ne(a, n) 25 Mg and 22 Ne(a, 7 ) 26 Mg rates 
(see discussion in Section 4.1). 

The main component is listed in column 2. All values are normalised to 150 Sm (see Section 5.1.1). In columns 3 to 6 we list 
the results of two tests on the 22 Ne(a, n) 25 Mg rate discussed in Section 4.1: Test A corresponds to column 3 (recommended 
22 Ne(a, n)x4, top panel of Fig. 6), Test B corresponds to column 5 (lower limit estimated in Section 4, bottom panel of 
Fig. 6). The s-only isotopes are highlighted in bold, while s-only nuclides ( 80 Kr, 86 Sr, 128 Xe, 152 Gd) and s isotopes ( 90 Zr, 96 Zr, 
164 Er, 180 Ta m ) with significant p-process contributions are written in italic. Isotopes with negligible s-process contribution 
(<5%) are excluded from this Table. Uncertainties (in %) listed for s-only isotopes in column 7 refer to solar abundances by 
Lodders et al. (2009). 

A complete version of Figs. 3, 6 and 7 is displayed in Figs. Al, A2, and A3, respectively. 

Table Al. 


Isotope 

Main component (in %) 
Recomm. 22 Ne(a, n) 

Test A (in %) 

Recomm. 22 Ne(a, n)x4 

Ratio 

(3)/(2) 

Test B (in %) 

Recomm. 22 Ne(a, n)x0.9 

Ratio 

(5)/(2) 

Solar 

uncertainty (in %) 

70 Ge 

6.0 

4.9 

.81 

6.1 

1.02 

10 

72 Ge 

7.5 

6.2 

.82 

7.6 

1.01 


73 Ge 

7.3 

6.2 

.84 

7.4 

1.00 


74 Ge 

8.4 

7.4 

.89 

8.4 

1.00 


75 As 

5.9 

5.3 

.90 

5.9 

1.00 


70 Se 

13.2 

11.7 

.88 

13.3 

1.01 

7 

77 Se 

6.5 

5.6 

.87 

6.5 

1.01 


78 Se 

14.3 

12.7 

.88 

14.4 

1.00 


80 Se 

8.0 

7.6 

.95 

8.0 

1.00 


79 Br * 

7.9 

7.5 

.96 

7.8 

.99 


81 Br* 

10.3 

9.1 

.88 

10.5 

1.01 


SO Kr 

9.9 

5.7 

.58 

10.5 

1.06 

— 

82 Kr 

23.3 

19.5 

.84 

23.5 

1.01 

— 

83 Kr 

8.2 

7.0 

.85 

8.3 

1.00 


84 Kr 

10.5 

10.4 

.99 

10.4 

.99 


86 Kr 

16.4 

25.8 

1.57 

16.1 

.98 


85 Rb 

12.0 

15.2 

1.27 

11.7 

.97 


87 Rb 

18.3 

41.1 

2.25 

17.2 

.94 


86 Sr 

58.5 

39.1 

.67 

60.0 

1.03 

7 

87 Sr 

58.5 

37.7 

.64 

59.7 

1.02 

7 

88 Sr 

82.6 

72.3 

.87 

83.0 

1.00 


89 Y 

85.8 

85.6 

1.00 

85.5 

1.00 


90 Zr 

76.0 

69.1 

.91 

76.3 

1.00 


91 Zr 

85.4 

86.0 

1.01 

84.9 

.99 


92 Zr 

83.3 

89.5 

1.07 

82.8 

.99 


94 Zr 

107.4 

103.5 

.96 

107.8 

1.00 


96 Zr 

14.3 

86.8 

6.06 

12.0 

.84 


93 Nb* 

69.2 

7 4.3 

1.07 

69.0 

1.00 


95 Mo 

60.4 

61.7 

1.02 

59.7 

.99 


90 Mo 

99.2 

120.3 

1.21 

98.2 

.99 

10 

97 Mo* 

53.4 

65.1 

1.22 

53.4 

1.00 


98 Mo 

74.1 

70.5 

.95 

74.8 

1.01 


99 Ru* 

27.2 

25.0 

.92 

27.5 

1.01 


100 Ru 

105.5 

88.6 

.84 

106.7 

1.01 

6 

101 Ru 

16.8 

14.1 

.84 

17.0 

1.01 


102 Ru 

57.2 

48.1 

.84 

57.6 

1.01 


103 Rh 

15.6 

13.0 

.84 

15.7 

1.01 


104 Pd 

110.7 

92.7 

.84 

111.5 

1.01 

5 

105 Pd 

14.2 

12.0 

.84 

14.3 

1.01 


i06p d 

52.4 

45.0 

.86 

52.7 

1.00 


108 pd 

66.4 

58.4 

.88 

66.7 

1.00 


107 Ag* 

15.1 

13.0 

.86 

15.2 

1.01 


109 Ag 

29.7 

26.2 

.88 

29.8 

1.00 


110 Cd 

101.7 

91.3 

.90 

102.0 

1.00 

7 

m Cd 

32.7 

29.5 

.90 

32.7 

1.00 


112 Cd 

65.1 

60.0 

.92 

65.1 

1.00 


113 Cd 

37.1 

34.3 

.92 

37.2 

1.00 


114 Cd 

76.6 

72.2 

.94 

76.5 

1.00 



Continued on next page 
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Table Al- Continued from previous page 


Isotope 

Main component (in %) 
Recomm. 22 Ne(a, n) 

Test A (in %) 

Recomm. 22 Ne(o-, n)x4 

Ratio 

(3)/(2) 

Test B (in %) 

Recomm. 22 Ne(a, n)x0.9 

Ratio 

(5)/(2) 

Solar 

uncertainty (in %) 

11B Cd 

8.6 

20.0 

2.33 

7.7 

.89 


H5in 

38.6 

36.4 

.94 

38.6 

1.00 


116 Sn 

87.7 

82.8 

.94 

87.7 

1.00 

15 

117 Sn 

48.9 

47.2 

.97 

48.9 

1.00 


118 Sn 

69.3 

70.0 

1.01 

69.1 

1.00 


119 Sn 

58.5 

58.7 

1.00 

58.4 

1.00 


120 Sn 

76.2 

76.3 

1.00 

76.1 

1.00 


122 Sn 

32.6 

64.3 

1.97 

31.2 

.96 


121 Sb 

37.6 

37.2 

.99 

37.7 

1.00 


122 Te 

89.1 

88.1 

.99 

89.1 

1.00 

7 

123 Te 

89.5 

88.4 

.99 

89.6 

1.00 

7 

124 Te 

92.7 

93.7 

1.01 

92.8 

1.00 

7 

125 Te 

20.6 

20.6 

1.00 

20.7 

1.00 


126 Te 

42.7 

41.4 

.97 

42.7 

1.00 


128 Xe 

894 

84.8 

.95 

89.6 

1.00 

- 

130 Xe 

98.2 

92.0 

.94 

98.5 

1.00 

— 

131 Xe 

7.4 

7.0 

.94 

7.4 

1.00 


132 Xe 

27.5 

27.9 

1.01 

27.4 

.99 


133 Cs 

13.9 

14.0 

1.01 

13.8 

1.00 


134 Ba 

115.0 

103.8 

.90 

116.6 

1.01 

6 

135 Ba* 

26.3 

29.6 

1.13 

26.0 

.99 


136 Ba 

108.9 

101.5 

.93 

109.2 

1.00 

6 

137 Ba 

60.3 

70.0 

1.16 

60.1 

1.00 


138 Ba 

93.5 

89.9 

.96 

93.5 

1.00 


139 La 

76.2 

73.3 

.96 

76.3 

1.00 


140 Ce 

95.3 

85.2 

.89 

95.8 

1.00 


142 Ce 

9.5 

37.5 

3.97 

8.1 

.86 


141 Pr 

50.4 

49.6 

.98 

50.0 

.99 


142 Nd** 

104.5 

88.0 

.84 

105.4 

1.01 

5 

143 Nd 

33.4 

32.1 

.96 

33.4 

1.00 


144 Nd 

52.8 

51.1 

.97 

52.8 

1.00 


145 Nd 

26.4 

25.8 

.98 

26.4 

1.00 


146 Nd 

67.3 

63.5 

.94 

67.6 

1.00 


148 Nd 

15.6 

20.1 

1.29 

14.7 

.94 


147 Sm 

27.1 

25.3 

.93 

27.2 

1.00 


148 Sm 

105.3 

100.9 

.96 

105.9 

1.01 

5 

149 Sm 

12.9 

12.9 

1.01 

12.8 

1.00 


150 Sm*** 

100.0 

100.0 

1.00 

100.0 

1.00 

5 

152 Sm 

22.5 

23.1 

1.02 

22.4 

1.00 


151 Eu 

6.0 

6.0 

1.00 

6.0 

1.00 


153 Eu 

6.1 

6.2 

1.02 

6.0 

1.00 


152 Gd 

854 

69.0 

.81 

87.6 

1.03 

5 

154 Gd 

92.0 

95.4 

1.04 

91.6 

1.00 

5 

155 Gd 

5.8 

6.1 

1.04 

5.8 

1.00 


156 Gd 

17.6 

17.9 

1.02 

17.6 

1.00 


157 Gd 

10.9 

11.0 

1.02 

10.9 

1.00 


158 Gd 

27.1 

27.4 

1.01 

27.2 

1.00 


159 Tb 

7.9 

7.9 

1.01 

7.9 

1.00 


160 Dy 

89.0 

89.7 

1.01 

89.4 

1.00 

5 

161 Dy 

5.2 

5.3 

1.01 

5.2 

1.00 


1B2 Dy 

15.9 

15.9 

1.00 

16.0 

1.01 


164 Dy 

24.0 

20.9 

.87 

24.5 

1.02 


165 Ho 

8.5 

8.2 

.96 

8.6 

1.01 


iSiEr 

87.1 

90.8 

1.04 

87.0 

1.00 

5 

166 Er 

15.9 

15.8 

.99 

16.0 

1.00 


l6 7 Er 

9.7 

9.6 

1.00 

9.7 

1.00 


168 Er 

32.7 

30.2 

.92 

33.2 

1.01 


170 Er 

12.6 

12.1 

.96 

12.2 

.97 


169 Tm 

14.8 

13.3 

.90 

15.0 

1.02 


170 Yb 

99.3 

89.5 

.90 

100.6 

1.01 

5 

171 Yb 

23.9 

20.3 

.85 

24.2 

1.01 


172 Yb 

44.7 

40.7 

.91 

44.7 

1.00 


173 Yb 

28.0 

26.3 

.94 

28.0 

1.00 


174 Yb 

61.9 

56.9 

.92 

62.2 

1.00 



Continued on next page 
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Tabic A1- Continued from previous page 


Isotope 

Main component (in %) 
Recomm. 22 Ne(o;, n) 

Test A (in %) 

Recomm. 22 Ne(n, n)x4 

Ratio 

(3)/(2) 

Test B (in %) 

Recomm. 22 Ne(o, n)x0.9 

Ratio 

(5)/(2) 

Solar 

uncertainty (in %) 

17B Yb 

6.0 

8.4 

1.40 

5.6 

.93 


175 Lu 

18.4 

16.8 

.91 

18.5 

1.01 


176 Lu 

112.4 

92.0 

.82 

115.0 

1.02 

5 

176 Hf 

101.5 

105.4 

1.04 

100.6 

.99 

5 

177 Hf 

17.4 

16.4 

.95 

17.5 

1.01 


178 Hf 

58.7 

53.3 

.91 

59.1 

1.01 


179 Hf 

42.1 

38.2 

.91 

42.3 

1.01 


1S0 Hf 

91.6 

83.6 

.91 

92.0 

1.00 


180 Ta m 

67.6 

59.2 

.88 

67.6 

1.00 

10 

181 Ta 

47.7 

43.6 

.91 

47.9 

1.00 


i82 w * 

64.9 

61.0 

.94 

65.2 

1.00 


18 3yy 

65.0 

56.2 

.87 

65.9 

1.01 


1 S4'W 

82.4 

71.2 

.86 

83.0 

1.01 


i86 W 

42.3 

42.5 

1.00 

41.2 

.98 


185 Re 

28.3 

24.2 

.85 

28.5 

1.01 


187 Re 

9.3 

9.6 

1.04 

9.0 

.97 


186 Os 

105.8 

98.0 

.93 

106.2 

1.00 

8 

187 Ost 

38.2 

36.2 

.95 

38.5 

1.01 

8 

188 Os 

28.2 

28.4 

1.00 

28.1 

.99 


190 Os 

14.7 

14.1 

.96 

14.8 

1.00 


192 Pt 

80.6 

78.1 

.97 

81.1 

1.01 

8 

IMpt 

6.6 

6.2 

.94 

6.6 

1.00 


wept 

14.0 

13.2 

.94 

14.1 

1.00 


197 Au 

6.1 

5.8 

.94 

6.2 

1.00 


198 Hg 

84.3 

80.2 

.95 

84.4 

1.00 

20 

i"Hg 

22.1 

21.0 

.95 

22.1 

1.00 


200 Hg 

53.1 

49.6 

.94 

53.1 

1.00 


201 Hg 

40.5 

38.4 

.95 

40.5 

1.00 


202 Hg 

71.8 

67.1 

.94 

72.0 

1.00 


204 Hg 

8.5 

5.8 

.69 

00 

bo 

1.03 


203rpj 

80.9 

79.0 

.98 

80.7 

1.00 


205rpp 

70.8 

50.9 

.72 

72.4 

1.02 


204p b 

91.3 

89.9 

.99 

91.1 

1.00 

7 

206PI-J** 

70.2 

64.2 

.92 

70.4 

1.00 


207 Pb 

61.3 

58.5 

.96 

61.3 

1.00 


208 pk 

47.4 

45.9 

.97 

47.5 

1.00 


209 Bi 

6.3 

9.5 

1.52 

6.1 

.97 



5 T1 include the decay of their short-lived isobars 79 Se. 


81 Kr, 93 Mo and 93 Zr, 97 Tc, "Tc, 107 Pd, 135 Cs, 182 Hf, 205 Pb. 

(**) 142 Nd and 206 Pb account for the a-decay of 146 Sm and 210 Bi. 

(***) Values are normalised to 150 Sm: we obtain X( 130 3m) jX rru ( 130 Sin) = 1133.9 
and 22 Ne(a, 7 ) 26 Mg rates (column 2). Test A provides X( ls 0 Sm)/Xi n j( 150 Sm) = 

X( 159 Sm)/Xi 7 l i( 150 Sm) = 1128.4 (—1%). See Table 3 of the manuscript. 

(t) 187 Os is a nominal s-only isotope, despite it receives a cosmogenic contribution by the decay of the very long-lived 187 Re. 


with the recommended 22 Ne(a, n) 25 Mg 
1275.5 (+12.5%), and Test B yields 


Noteworthy isotopes are 96 Mo and 176 Lu: with TEST A the s contribution to 96 Mo is about 20% higher than solar, while 
176 Lu is clearly underestimated (about 9% of solar 176 Lu is missing). This may suggest that the adopted conservative upper 
limit of the recommended 22 Ne(a, n) 25 Mg rate (x4) is overestimated. However, a large uncertainty is associated to these 
isotopes (e.g., the branches at 95 Zr and the unknown contribution of isomeric 176 Lu mediating states). 

Apart from some s isotopes, several nuclides show variations larger than 10%. Among them, we highlight 96 Zr (whose s 
prediction increases from 14% to 87%), 13 'Ba and 205 T1, which receive a dominant s-process contribution, and 86 Kr, 87 Rb, 
116 Cd, 122 Sn, 135 Ba, 142 Ce, 148 Nd, 185 Re. We remind that 86 Kr, 87 Rb, 96 Zr (and 137 Ba) are also produced in IMS stars (see 
Section C). Note that the neutron-rich isotopes 176 Yb and 208 Hg, as well as 209 Bi are marginally produced by the main 
component (<10%). In Fig |A4l we display the results of both tests for isotopes that receive an s-process contribution higher 
than 50%. 
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le+10 


le+09 


le+08 



Figure Al. Temporal evolution of the neutron density ( N n ) and 
the maximum temperature at the bottom of the convective zone 
(T) during the 15th He shell flash in AGB stars with 1.5 and 3 
Mq (thick and thin lines, respectively) and half-solar metallicity. 
The time-scale starts when the temperature at the bottom of the 
convective TP reaches Tg = 2.5, which corresponds to the onset 
of the 22 Ne(a, n) 25 Mg reaction. 


All isotopes with 90 < A < 205: 
main s-process component 



Atomic mass 

Figure A2. Ratios between the main component obtained with a 
two times higher and a three times lower 13 C(o!, n) ie O rate than 
our adopted rate shown in Fig 2 (Test 1 and Test 2 correspond¬ 
ing to the top and bottom panels). This is a complete version of 
Fig. 3 showing all isotopes from 90 ^ A ^ 210. 



80 90 100110120130140150160170180190200210 
Atomic mass 


Figure A3. Ratios between the main component obtained within 
a conservative range of uncertainty associated to the 22 Ne(a, 
n) 25 Mg and 22 Ne(a, 7 ) 26 Mg rates, as discussed in Section 4.1. 
This is a complete version of Fig. 6 showing isotopes from 90 ^ A 
^210 that receive an s-process contribution higher than 50%. As 
in Fig. 6 , an improved treatment for the branches close to 134 Ba, 
152 ’ 1 54 Gd? 164 Er? 180 Xa m and 204p b ig i nc l uded . 



Atomic mass 

Figure A4. Ratios between the main component obtained with 
the recommended 22 Ne(a, n) 25 Mg and 22 Ne(a, 7) 26 Mg rates 
compared to Test C in which the upper limit for the 22 Ne(o;, 
7 ) 26 Mg rate is adopted, while the 22 Ne(o:, n) 25 Mg reaction is 
unchanged. This is a complete version of Fig. 7 showing all iso¬ 
topes from 90 ^ A ^ 210. 
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Figure Bl. Schematic representation of the chart of nuclei in 
the atomic mass region close to the s-only isotope 1 60 1)y (red 
square). Thick lines represent the s-process nucleosynthesis dur¬ 
ing the 13 C(a, n) le O neutron irradiation, while thin lines corre¬ 
spond to the neutron capture channels. The terrestrial half-life of 
160Tb (ti /2 = 72.3 d) is reduced to 1.45 d at Tg = 2 and to 11 h 
at Ts = 3 (Takahashi & Yokoi 1987). Thus, neutron captures are 
marginally activated during the 22 Ne(«, n) 25 Mg neutron burst, 
and 160 Dy is almost unbranched. ( This and the following Figures 
are adapted from http://www.kadonis.org/. See the electronic edi¬ 
tion of the Journal for a colour version of the Figures.) 

APPENDIX B: ADDITIONAL BRANCHES 

This Appendix completes the information about the most 
important branches of the main component discussed in 
the paper. Following the classification given in Section 5, 
we analyse here the additional branchings of each category, 
which are only briefly outlined in the text. 

Bl Additional unbranched isotopes 

Bl.l The s-only isotope 160 Dy 

The s-only isotope 160 Dy (2.34% of solar Dy) is mainly syn¬ 
thesised during the 13 C-pocket phase and remains almost 
unbranched during TPs (Fig. IB1I) . The related branching 
points 159 Gd and 160 Tb are both short-lived at s-process 
temperatures. The decay of 159 Gd (ti /2 = 19 h) is constant 
under AGB conditions, while the terrestrial half-life of 160 Tb 
is strongly reduced at Tg = 3 (from tj / 2 = 72.3 d to 11 h; 
Takahashi & Yokoi 1987). Thus, less than 2% of the s flow 
is bypassing 160 Dy. 

The s contribution to solar 160 Dy is 90%. 

The solar abundances of Dy and Sm are both uncertain by 
5% (Lodders et al. 2009). The 160 Dy MACS is very well 
determined (cr[ 160 Dy(n, 7 )] = 890 ± 12 mbarn at 30 keV, 
1.4%). However, its large stellar enhancement factor (SEF 
= 1.12 at 30 keV) suggests that the neutron capture rate of 
160 Dy may be affected by an additional theoretical uncer¬ 
tainty. 

Given the negligible activation of the branches at 159 Gd and 
160 Tb, the s contribution to 160 Dy is only marginally af¬ 
fected by the MACS uncertainties of these branching points 
(e.g., ~2% variation for a 16% change of the 160 Tb MACS, 
3240 ± 510; KADoNiS). Similarly, the large 160 Tb /3-decay 
uncertainty (up to a factor of three at Ts = 3; Goriely 1999) 
has negligible effects on 160 Dy. Accordingly, an increase of 
the recommended 22 Ne(a, n) 25 Mg rate by factor of four af¬ 
fects the s contribution of 160 Dy by less than 2%. 
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Figure B2. The same as Fig. ED but for the s path close to the 
s-only isotope 198 Hg (red square). 

In conclusion, the prediction of the solar abundance of 
160 Dy is mainly affected by the theoretical evaluation of 
MACS of 160 Dy (see also Rauscher et al. 2011). 

B1.2 The s-only isotope 196 Hg 

The s-only isotope 198 Hg is almost unbranched (Fia, IB2l) . be¬ 
cause the short half-lives of 19 'Pt and 198 Au strongly favour 
the /3-decay channels compared to neutron captures. The 
half life of 197 Pt (ti /2 = 19.89 h) is rather constant under 
stellar conditions, while the terrestrial half-life of 198 Au (H /2 
= 2.7 d) decreases by a factor of four during TP (Takahashi 
& Yokoi 1987). 

Accordingly, the s prediction for 198 Hg varies by less than 
2% with increasing the recommended 22 Ne(a, n) 25 Mg rate 
by a factor of four. 

The major nuclear uncertainty derives from the 198 Hg 
MACS (8.7%, 173 ± 15 mbarn; KADoNiS). 

About 84% of solar 198 Hg is produced by the main com¬ 
ponent. We remind that the solar Hg abundance is largely 
uncertain for this very volatile element. Originally, it was 
estimated using the s-process systematics via the relation 
crN s = constant (N(Hg) = 0.34 ± 0.04, Si = 10 6 , Beer and 
Macklin, 1985). While this value was adopted in the abun¬ 
dance compilation of Anders and Grevesse (1989), Lodders 
(2003) determined the Hg abundance as a weighted average 
of measurements from two meteorites (Orgueil and Ivuna) 
with 50% of uncertainty due to the difficulty of finding clean 
samples. Recently, Lodders et al. (2009) improved the mea¬ 
sured Hg abundance, claiming an uncertainty of 20% (N(Hg) 
= 0.458 ± 0.092, Si = 10 6 ). This higher solar Hg abundance 
explains the lower s contribution of 198 Hg (84%) compared 
to previous evaluations (~100% using the solar Hg abun¬ 
dance of Anders and Grevesse 1989; see e.g., Kappeler et al. 
2011, their Fig. 15). The 20% difference is, however, com¬ 
patible with the respective solar and nuclear uncertainties. 

B2 Additional branches sensitive to the neutron 
density 

B2.1 The s-only isotopes 86,87 Sr (the branch points at 
S5 Kr and 86 Rb) 

The abundances of the s-only isotopes are sensitive to the 
branching at 85 Kr and to a lesser extent to that at 86 Rb (see 











32 S. Bisterzo, et al. 



Figure B3. The same as Fig. lBlI but for the s path close to the 
s-only isotopes 86,87 Sr. The branch at 85 Kr is represented by a 
red circle. The neutron-magic nuclei at N = 50 are indicated by 
a blue box. 


f^. nu. 

In fact, the branching starts already with neutron captures 
on 84 Kr, which are feeding either the short-lived isomeric 
state of 85 Kr m ( t 1 / 2 = 4.48 h) or the long-lived ground state 
of 85 Kr s (fj /2 = 10.73 yr). The two 85 Kr states are not ther- 
malised and have to be treated independently (Ward 1977). 
At least ~80% of the short-lived 85 Kr m /3~ decay to 85 Rb, 
while the remaining 20 % decay to 85 Kr E by internal transi¬ 
tions. Therefore, about 40% of the s flow coming from 84 Kr 
feed 85 Rb by /3~ decay of 85 Kr m , and about 60% produce 
85 Kr g . 

Due to the short half-life of 86 Rb (ti /2 = 18.6 d), this addi¬ 
tional branching is only partially open at high neutron den¬ 
sities. The combined effect of both branchings is completely 
determined by the neutron density, because the half-lives of 
the branch point isotopes are insensitive to the temperatures 
during the AGB phase (Takahashi and Yokoi 1987). 

B2.1.1 During the 13 C pocket. The neutron densities 
typical of the 13 C pocket are too low to activate the neutron 
captures branch at 85 I<r efficiently. The isomer with a half- 
life of the order of a few hours immediately B~ decays to 
85 Rb, and the ground state mainly /3 - decays as well in spite 
of its longer half-life (f n < 0.1 even at the maximum neutron 
density of ~2xl0 7 cm -3 ). This means that less than ~10% 
of the s-process flow proceed towards 86 Kr, which receives a 
correspondingly small s-process contribution during the 13 C 
pocket (~10% of solar 86 Kr). 

The second branch point at 86 Rb is not active during the 
13 C pocket and the s-flow directly feeds 86 - 87 ’ 88 Sr. Only ~7% 
of solar 87 Rb is produced during the 13 C(a, n) 16 0 neutron 
irradiation. 

The present contributions of the main component to 
86 Kr and 87 Rb are 16% and 18%, respectively (Table Al). 
By including the recent 13 C(a, n) le O rate by La Cognata 
et al. (2013), these values increases to 19% and 20%, re¬ 
spectively, (see Table IbTI) . and are given with less than 3% 
uncertainty. By assuming an extended range of uncertainty 
based on analyses given in literature, variations up to 20 - 
35% are obtained by changing the 13 C(a, n) 16 0 rate by 
factors of two and three. 

All stable Br, Kr and Rb isotopes have well known (n, 
7 ) cross sections (Mutti et al. 2005, Heil et al. 2008a). Conse- 


Table Bl. Updated solar s contributions to 86 Kr and 87 Rb us¬ 
ing our 13 C(ct, n) le O rate by Denker et al. (1993; “D93”, see 
Table Al) and the 13 C(ct, n) le O rate by La Cognata et al. (2013; 
“L13”). The present s contributions are compared with those eval¬ 
uated by changing the 13 C(a, n) le O rate by factors of two and 
three. 


D93 

L13 

D93x2 

D93/3 

86 Kr 16.4 

18.6 

19.6 (+19.5%) 

10.7 (-34.7%) 

S7 Rb 18.3 

19.8 

20.5 (+12.0%) 

13.7 (-25.1%) 


quently, the major uncertainties affecting the s-process pre¬ 
dictions in this mass region are due to the branch point nu¬ 
clei 85 Kr and 86 Rb, for which only theoretical estimations 
are available. The ~80% uncertainty of the 85 Kr(n, 7 ) 86 I<r 
MACS recommended by KADoNiS (55 ± 45 mbarn at 30 
keV; Bao et al. 2000) has been recently reduced: Raut et al. 
(2013) have measured photo-disintegration and ( 7 , 7 ’) cross 
sections on 86 Kr for improving the prediction of the 85 Kr(n, 
7 ) cross section. Their recommended MACS of 83 mb 
is about 50% higher than in Bao et al. (2000) with a two 
times smaller uncertainty. Nevertheless, the 85 Kr MACS by 
Bao et al. (2000) and Raut et al. (2013) are in agreement 
within the respective uncertainties. 

The MACS uncertainty of 85 Kr of Raut et al. (2013) re¬ 
sults in variations of 16-30% and 10-20% in the contribu¬ 
tion of the main component to 86 Kr and 87 Rb during the 
13 C-pocket phase, respectively. 

B2.1.2 During thermal pulses. The neutron density 
reached during TPs by the 22 Ne(a, n) 25 Mg reaction (up to 
a few 10 10 cm -3 ) allows the activation of both branches at 
85 Kr and 86 Rb, largely overcoming the uncertainties affect¬ 
ing 86 Kr and 87 Rb discussed above. 

About 40% of the s flux coming from 84 Kr feeds 85 Rb by 
/3 - decay of 85 Kr m (ti / 2 = 4.48 h), and 60% produces 85 Kr g 
(ti /2 = 10.73 yr). The longer-lived 85 Kr g easily captures 
neutrons, and the s path proceeds towards 86 Kr, 87 Rb (a 
very long-lived isotope with t 1 / 2 = 4.8xl0 10 yr), and 88 Sr, 
bypassing 86 , 87 Sr. At the end of the convective instability, all 
the 85 Kr g stored during TP decays into 85 Rb (see Fig. IB4I) . 
At 86 Rb (ti /2 = 18.7 d), neutron capture and /3~ decay com¬ 
pete for neutron densities larger than ~ 10 10 cm -3 , further 
increasing the 87 Rb abundance. The neutron-magic isotopes 
86 Kr, 87 Rb and 88 Sr (at N = 50) operate as bottlenecks of 
the s path, owing to their very small neutron capture cross 
sections, thus enhancing the effect of the branches at 85 Kr g 
and 86 Rb. 

86 Kr and 87 Rb are largely affected by the 22 Ne(a, n) 25 Mg 
rate (see Table Al): the s-contributions to 86 Kr and 87 Rb 
increase by a factor of 1.6 and 2.3 by increasing the rec¬ 
ommended 22 Ne(a, n) 25 Mg rate by a factor of four. The 
86,87 Sr s-predictions are less affected than 86 Kr and 87 Rb, 
with decreases of ~30%. The overabundances of 88 Sr, which 
has been accumulated during the 13 C(a, n) 16 0 neutron ir¬ 
radiation (MACS at 30 keV is 6.13 ± 0.11 mb; KADoNiS), 
is almost unchanged. Owing to the normalisation effect to 
150 Sm (Section 5.1.1), the s contribution to 88 Sr shows vari¬ 
ations of 12.5%. 

Additional uncertainties affecting 86 Kr and 85 ' 87 Rb de- 
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Figure B4. Temporal evolution of the neutron density, and the 
isotopic abundances of 86 ’ 87 Sr, 85 Kr^, 86 Kr and 85,86,87 Rb dur _ 
ing the 15 th He shell flash in an AGB star with 1.5 Mq and 
half-solar metallicity. The abundance values are given as number 
fractions normalised to 150 Sm at the end of the He shell flash. 
The time-scale starts when the temperature at the bottom of the 
convective TP reaches Tg = 2.5, which corresponds to the onset 
of the 22 Ne(a, n) 25 Mg reaction. 


Table B2. The abundance contributions by the main component 
(in %) to the isotopes from 84 Kr to 88 Sr compared with the corre¬ 
sponding values obtained with the upper and lower MACS limits 
of 85 Kr and 86 Rb (Raut et al. 2013, ”R13”; KADoNiS, ”KAD”). 
Variations higher than 10% are listed between brackets. 



Tests of the 85 Kr(n, 7 ) 

86 Kr MACS 

Isotope 

[R 13 ] 

[R13]LL 

[R13]UL 

84 Kr 

10.5 

10.5 

10.6 

86 Kr 

16.4 

9.5 (0.58) 

20.4 (1.24) 

85 Rb 

12.0 

13.4 (1.12) 

11.4 (0.95) 

87 Rb 

18.3 

13.8 (0.75) 

21.0 (1.15) 

86 Sr 

58.4 

59.2 

58.6 

87 Sr 

58.5 

59.1 

58.7 

88 Sr 

82.6 

83.3 

83.1 


Tests of the s6 Rb(n, 7 ) 87 Rb MACS 


Isotope 

KAD 

KAD LL 

KAD UL 

84 Kr 

10.5 

10.5 

10.5 

86 Kr 

16.4 

16.5 

16.6 

85 Rb 

12.0 

12.0 

12.0 

87 Rb 

18.3 

15.7 (0.86) 

21.9 (1.20) 

86 Sr 

58.4 

60.3 

56.8 

s7 Sr 

58.5 

59.8 

57.7 

88 Sr 

82.6 

83.2 

83.2 


rive from the theoretical MACS values of 85 Kr and 86 Rb 
(Raut et al. 2013; KADoNiS). In Table lB2l we evaluate the 
effects of the MACS uncertainties discussed above on the 
s predictions of the isotopes from 84 I<r to 88 Sr. While the 
MACS uncertainty of 86 Rb influences only 87 Rb up to 20%, 
the one of 85 Kr affects 86 I<r as well as 87 Rb (variations up 
to 40 and 25%, respectively). Also 85 Rb changes by 5-10%. 
The s-only isotopes 86,87 Sr show smaller variations (($2%). 

Note that the solar s abundances of isotopes with A < 
90 receive additional contributions. 
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Figure B5. Schematic representation of the chart of nuclei in 
the atomic mass region close to the s-only isotope 142 Nd (red 
square). The neutron-magic nuclei at N = 82 are indicated by a 
blue box. 


Firstly, the weak s process in massive stars synthesises iso¬ 
topes from Cu up to Kr-Rb-Sr (Rauscher et al. 2002; Pig- 
natari et al. 2013; Chieffi & Limongi 2013). In this case, the 
s predictions are largely affected by the 12 C + 12 C rate (see, 
e.g., Pignatari et al. 2013). 

Second, intermediate-mass AGB stars contribute to solar 
86 Kr and 87 Rb. The s contribution is difficult to be assessed 
because of the large uncertainty affecting the treatment of 
mass loss in IMS AGB stars. Bisterzo et al. (2014) estimate 
that the predicted contribution to the solar 86 Kr abundance 
may increase from 15% up to ~+35% and a still larger effect 
is found for 87 Rb (from ~25% up to ~70%). In Section C 
(Supplementary Information), we discuss the result of a M 
= 5 Mq model at half-solar metallicity, in the light of the 
neutron source uncertainties. 

B2.2 The s-only isotope li2 Nd (the branches at 141 Ce and 
142 Prj 

142 Nd belongs to the magic nuclei at N = 82, and has a fairly 
small MACS (35.0 ± 0.7 mbarn; KADoNiS). The reproduc¬ 
tion of the solar 142 Nd abundance was the first evidence that 
the classical model was not adequate to describe the main 
component close to the neutron-magic numbers (N ~ 50, 
82, and 126), located at the deep steps of the <cr>N s curve 
(Seeger et al. 1965, Clayton 1968, Kappeler et al. 1989). The 
success of the stellar model in reproducing the s-process pat¬ 
tern of several branchings (e.g., 142 Nd, and 116 Sn, 134,136 Ba 
as well) was demonstrated by Arlandini et al. (1999). 

The s-only isotope 142 Nd is shielded against the re¬ 
process by 142 Ce. Most of solar 142 Nd is synthesised during 
the 13 C-pocket phase where the s path is restricted to the 
neutron-magic isotopes 140 Ce, 141 Pr, and 142 Nd (Fig. IB5ll 
because of the low neutron density. The branchings at 141 Ce 
and 142 Pr, which open at higher neutron densities, are af¬ 
fecting the s abundance of 142 Nd only during the TPs. 

In spite of the relatively long half-life of 141 Ce (ti /2 = 
32.51 d, constant at s-process temperatures), this branching 
opens only at high neutron densities (/„ = 0.2 for N„ = 
4xl0 9 cm -3 ) because of the small MACS of 76 ± 33 mb at 
kT = 30 keV. About 8% of the reaction flow proceed via 
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Figure B6. Schematic representation of the chart of nuclei in the 
atomic mass region close to the s-only isotope 186 Os (red square). 
The branch at 185 W, mainly activated during TPs, is represented 
with a red circle. 187 Os is shielded from the r -process by its long- 
lived isobar 187 Re (£ 1/2 = 4.12xl0 10 yr), which contributes to 
its cosmoradiogenic enrichment. 
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Figure B7. The same as Fig. IB-41 but for the isotopic abundances 
G f 186,187,189 q s 185,186,187,188 J^ e an( j 185,186,187^ 


141 Ce, contributing ~20% to the solar abundance of 142 Ce. 
The second branch at 142 Pr is characterized by a shorter 
half-life (19.25 h, increasing by a factor of four at Tg = 3; 
Takahashi and Yokoi 1987), only partially compensated by a 
larger MACS of 415 ± 178 mb (KADoNiS). Therefore, this 
branching remains weak even during TPs, reaching f n > 0.1 
only for neutron densities N n > 3xl0 9 cm -3 . 

The present contribution of the main component to the 
142 Nd abundance is 104.5%. Note that this value includes 
the small contribution from the a decay of the long-lived 
146 Sm (tj /2 = 103 Myr). 

The uncertainties of the theoretical MACS values of 
141 Ce and 142 Pr (43% each) are causing variations of about 
4% in the 142 Nd abundance. A similar effect of 5% is found 
if the 22 Ne(a, n) 25 Mg rate is increased by a factor of four. 
The s contribution to solar 142 Nd shows up to 15% varia¬ 
tions (owing to the additional effect of the 22 Ne(a, n) 25 Mg 
tests of the 150 Sm s-prediction; Section 5.1.1). 

B2.3 The s-only isotope 186 Os (the branches at 185 W and 
1S6 Re) 

The abundance of 186 Os is mainly regulated by the branch 
at 185 W, which is efficiently activated during TPs CFig. IBGI) . 
Actually, the main component overestimates the solar abun¬ 
dance of 186 0s by 30%, essentially due to the recommended 
MACS for the branch point nucleus 185 W (KADoNiS). This 
recommendation represents a weighted average between two 
measurements of the inverse 186 W(y, n) 185 W reaction by 
Sonnabend et al (2003; cr( 185 W) = 687 ± 110 mbarn at 30 
keV) and Mohr et al. (2004; <r( 185 W) = 553 ± 60 mbarn at 
30 keV). 

This problem was firstly emphasized by Sonnabend et al. 
(2003), who found a 20% excess of 186 Os. By adopting the 
upper limit of this experimental value, the s abundance of 
186 Os is still overestimated by 14%, incompatible with the 
8% uncertainty of the solar Os abundance (Lodders et al 
2009). 

The half-life of 185 W decreases slightly during TPs, from 
ii /2 = 75.1 d to ~53 d at Tg = 3 (Takahashi & Yokoi 1987, 
depending on electron density). The 10% uncertainty of the 


stellar half-life estimated by Goriely (1999) is also not suffi¬ 
cient to explain the predicted overabundance of 186 Os. 

The neutron cross sections of the Os isotopes are well 
determined by Mosconi et al. (2010): the 186 Os MACS is 
known at 4.1% (Mosconi et al. 2010, 430.6 ± 17 mbarn), 
and produces marginal variations of the 186 Os s-contribution 
(2%). Possible uncertainties related to the SEF estimation 
have been substantially reduced by Fujii et al. (2010). 

On the other hand, the second branch point at 186 Re 
is only marginally activated and can also be excluded as 
an explanation of the 186 0s problem. Indeed, 186 Re decays 
quickly into 186 Os (ti /2 = 4.09 d) and neutron captures are 
favoured only at very high neutron density (/„ > 0.5 for N n 
> 4xl0 9 cm -3 ; see Fig. IB7I) . The 186 Re /3 + -decay channel 
is negligible. Moreover, 186 Re does not build up to a signifi¬ 
cant abundance because of its large calculated MACS (1550 
± 250 mb, 16% estimated uncertainty) and has, therefore, 
a negligible effect on the 186 Os s-abundance (~1% ). 

Also the long-lived 186 Re isomer does not provide any ap¬ 
preciable contribution at s-process temperatures, and the 
branch at 186 Re remains small (Mohr et al. 2008, Hayakawa 
et al. 2005). Finally, the experimental photodisintegration 
measurement of the inverse 187 Re( 7 , n) 186 Re reaction car¬ 
ried out by Muller et al. (2006) provides a 186 Re MACS 
similar to that given by Bao et al. (2000). 

In order to reproduce the solar 186 Os abundance, we 
suggest that the MACS of 185 W should be 50% higher than 
the first value recommended by Bao et al. (2000). This guess 
is based on the fact that the uncertainty of the 185 W MACS 
may be underestimated because only the inverse ( 7 , n) reac¬ 
tion has been measured. Rauscher (2014) estimates that our 
guess may be plausible. However, specific theoretical anal¬ 
yses on the 185 W(n, 7 ) reaction should be carried out to 
provide a more realistic uncertainty. 

The 186 Os s-contribution, computed by including the above 
guess on the 185 W MACS, is overestimated by 6 %, in agree¬ 
ment within the uncertainties. By increasing the recom¬ 
mended 22 Ne(a, n) 25 Mg rate by a factor of four, the solar s 
186 0 s is reduced to ~98%. 

187 Os is shielded from the r-process by its stable iso- 
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Figure B8. Schematic representation of the chart of nuclei in the 
atomic mass region close to the s-only isotope 192 Pt (red square). 
The branch at 192 Ir, mainly activated during TPs, is represented 
with a red circle. 
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Figure B9. The same as Fig. [B4l but for the isotopic abundances 
Q f 192,193,194p t 190,191,192 Q g 191,192,193 


bar ls 'Re. Apart from neutron captures on 186 0s, 18 'Os is 
also produced by the cosmoradiogenic /3~ decay of 187 Re. 
The half-life of the ground state (fi /2 = 4.35 x 10 10 yr) is 
reduced in stellar interiors (Takahashi & Yokoi 1987; Bosch 
et al. 1996), but remains quasi stable within the time scale 
of TPs. The stable daughter 18 'Os becomes unstable under 
stellar conditions (Takahashi and Yokoi 1987), but behaves 
also like a stable isotope during the 6 yr duration of a TP. 
Therefore, its temporal evolution follows essentially that of 
186 Os, which is depleted during the peak neutron density 
and re-established towards the end of the TP (see Fig. IB7I) . 
A comprehensive investigation of the 187 Os/ 18 ' Re ratio with 
chemical evolution models will be the topic of a future work. 


B2-4 The s-only isotope 192 Pt (the branch at 192 Ir) 

The abundance of 192 Pt is mainly determined by the branch 
at 192 Ir ('Fig. IB8I) . At neutron densities higher than 2x10 s 
cm~ 3 , 192 Pt is mainly bypassed by neutron captures on 192 Ir 
(the terrestrial half-life of 77.54 d is reduced to 55 d at Tg 
= 3), to be partially restored (78% of the 192 Pt abundance 


produced during the previous 13 C pocket; see Fig. IB9I) at 
the end of TP. 

The branch at 191 Os has a smaller influence on 192 Pt, be¬ 
cause of the shorter half-life (the terrestrial value of 15.04 
d is reduced by a factor of two at Tg = 3, Takahashi and 
Yokoi 1987). Therefore, neutron captures on 191 Os (MACS 
= 1290 ± 280 mbarn at 30 keV; KADoNiS) can compete 
only at peak neutron densities (N n > 3.6xl0 9 cm -3 ; f n 
> 0.5). The s-process prediction underestimates the solar 
192 Pt by about 20%. 

Previously, the missing s contribution was compatible with 
the 20% uncertainties that were estimated for all theoretical 
192 Pt, 191 Os and 192 Ir MACS (cr( 192 Pt) = 590 ± 120 mbarn, 
cr( 191 Os) = 1290 ± 280 mbarn, and cr( 192 Ir) = 2080 ± 450 
mbarn at 30 keV; KADoNiS). The 192 Ir /1-decay rate has a 
marginal impact on the 192 Pt abundance. 

Recently, Koehler et al. (2013) measured the neutron cap¬ 
ture cross sections of Pt isotopes with much improved ac¬ 
curacy (e.g., er( 192 Pt) = 483 ± 20 mbarn at 30 keV; 4% 
uncertainty), and used their experimental results to provide 
a new theoretical estimation of the 192 Ir MACS (3220 ± 720 
mbarn at 30 keV). With the new MACS by Koehler et al. 
(2013) the uncertainty was significantly reduced but the s 
contribution of 81% remains too small. Allowing for a 20- 
uncertainty of the theoretical MACS of 192 Ir, the s contri¬ 
bution of 192 Pt can be increased to 95%, consistent with the 
8% uncertainty of the solar Pt abundance (Lodders et al. 
2009). 

More detailed analyses concerning the MACS of 192 Ir would 
help to improve the understanding of this branching point 
and to reproduce the solar abundance of 192 Pt more accu¬ 
rately. 

The impact of the 22 Ne(a, n) 25 Mg neutron source is 
marginal in this case. An increase of the recommended 
22 Ne(a, n) 25 Mg rate by a factor of four is affecting the s 
abundance of 192 Pt by less than 4%. 


B3 Additional branches strongly sensitive to 

stellar temperature (and/or electron density) 
and neutron density 

B3.1 The s-only isotopes 152,154 Gd (the branches at 
151 Sm and 154 Eu) 

The production of the s-only isotopes 152,154 Cd is affected 
by two main branches at 151 Sm and 154 Eu, respectively 

(Fig. mo}. 

During the 13 C-pocket phase most of the s path feeds 
the Sm isotopes from A = 150 to 153, mainly bypassing 
1S2 Gd (~0.2% of solar Gd). Although the terrestrial half-life 
of 151 Sm (ti /2 = 90 yr) may decrease to t \/2 = 55 yr (at Tg = 
1 and an electron density n e = 30xl0 26 cm~ 3 , Takahashi & 
Yokoi 1987), the neutron density is still sufficient to bypass 
152 Gd (N n > 1x10® cm -3 ; /„ > 0.5). 

At the peak neutron densities of the TPs ( N n > 10 9 cm' 3 , 
fn > 0.2) when the half-life of 151 Sm is reduced by a factor 
of 30 (ti /2 ~ 3 yr at Tg = 3), the branching at 151 Sm has 
a strong effect on 152 Gd. Accordingly, lsl Eu and 152 Gd are 
depleted near the maximum of N„, before they are partially 
restored during the freeze-out phase. As shown in Fig. IB 111 
the 151 Eu abundance is even increasing (~20%) above the 
value after the 13 C-pocket phase, because of the additional 
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Figure BIO. Schematic representation of the chart of nuclei in 
the atomic mass region close to the s-only isotopes 152 ’ 154 Gd (red 
squares). The terrestrial half-life of 151 Sm (£ 1/2 = 90 yr), which 
depends strongly on temperature and neutron density (Takahashi 
and Yokoi 1987), is reduced to a few years at Tg = 3. From 152 Eu, 
the s-process flow proceeds directly to 152 Gd, because the half-life 
of 152 Eu is reduced to about 5 hours at stellar temperatures. 
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Figure Bll. The same as Fig. lB4l but for the isotopic abun¬ 
dances Of 152,153,154,155 Gd! 151,152,153 Sm and 151,152,153,154 Eu 


feeding by the decay of 151 Sm after the end of the TPs. In 
contrast, 152 Gd, which feels the strong variations of neutron 
density and temperature in the He flash, reaches only about 
20% of the abundance produced in the 13 C pocket. 

Partial neutron captures on 151 Sm are allowed at the 
bottom of the advanced TPs when peak neutron densities 
are reached ( N n > 10 9 cm~ 3 , f n > 0.2), and 151 Eu and 
1S2 Gd are mostly destroyed (see Fig. IBlll . At the decreas¬ 
ing tail of the neutron density, most of 151 Sm present in the 
hot TP layers /3“ decays into 151 Eu, feeding again 152 Gd. 
The 151 Sm abundance stored at the end of the TP decays 
into 151 Eu during the interpulse phase, increasing the lsl Eu 
amount present in the previous 13 C pocket (~+20%). As 
shown in Fig. IBlll 152 Gd feels the strong variations of neu¬ 
tron density and temperature in the He flash, being de¬ 
stroyed at peak neutron density and partially rebuilt at the 
end of the TP: only ~20% of 152 Gd yielded during the pre¬ 
vious 13 C pocket is restored. 

The branch at 154 Eu has important effects on the s- 


prediction of 154 Gd. 

The terrestrial 154 Eu half-life (ti /2 = 8.8 yr) strongly de¬ 
creases during convective instabilities: from a few years 
down to a few days at the bottom layers of the advanced 
TPs (ii /2 = 10.67 d at Is = 3, Takahashi & Yokoi 1987). 
Given the large probability for 154 Eu to capture neutrons 
(MACS( 154 Eu) = 4420 ± 670 mb at 30 keV; KADoNiS), 
the s path mainly proceeds towards 155 Eu starting from 
N n ~9xl0 8 cm -3 (f n > 0.5), and bypasses 154 Gd. When 
the neutron density decreases, the abundance of 154 Gd pro¬ 
duced during the previous 13 C pocket is almost completely 
restored (~94%, Fig. IBlll) . 

The additional branchings at 153 Sm, 152 Eu, and 153 Gd 
are relatively weak. Almost the entire s-process flow runs 
from 152 Sm to 153 Eu due to the short half-life of 153 Sm (fi /2 
= 1.1 d at 7*8 = 3). While the rather slow terrestrial decay 
of 152 Eu (half-lives of 13.51 and 47.84 yr for the /3 + and EC 
channels, respectively), the /3~ decay dominates under stel¬ 
lar conditions (tf/ 2 = 20 and 5 h at 7s = 2 and 3, Takahashi 
and Yokoi 1987), directing the s-process flow from 151 Eu di¬ 
rectly to 152 Gd. 

The next unstable isotope of the s path is 153 Gd, which 
slowly /3 + decays with ti /2 = 240 d (rather constant un¬ 
der stellar conditions, Takahashi & Yokoi 1987). Due to the 
large MACS recommended for 153 Gd (4550 ± 700 mb at 
30 keV; KADoNiS) the reaction flow through 152 Gd contin¬ 
ues to 154 Gd already at low neutron densities (N n > 4x10' 
cm -3 ; f n > 0.5). 

Previous AGB models (see filled circles in Fig. 5) pre¬ 
dict that 81% of solar 152 Gd and 87% of solar 154 Gd are 
produced by the main component. These results were com¬ 
puted by adopting constant /1-decay rates for 151 Sm and 
154 Eu. These values increase to 85% and 92%, respectively, 
by including an improved treatment of the /3-decay rates of 
unstable isotopes close to 152 ’ 154 Gd ( lsl Sm, 154 Eu, 153 Sm, 
152 Eu, 153 Gd) over the convective TPs (see filled diamond 
in Fig. 5). 

However, the s predictions of 152,154 Gd are largely influ¬ 
enced by the uncertainties of the decay rates of 151 Sm and 
154 Eu. According to Goriely (1999) the rate of 151 Sm is 
increasing by factors of ~2 and 3 at 7s = 2 and 3, respec¬ 
tively, and a similar uncertainty is estimated for the 154 Eu 
/3-decay rate. This may produce up to 10-15% variations 
of the 154 Gd s-contribution, and an extreme impact on the 
152 Gd s-prediction (up to a factor of 2). 

Consequently, an appropriate treatment of the /3-decay rates 
close to 152,154 Gd improves the s predictions, but the uncer¬ 
tainty remains large. Investigation of the 151 Sm and 154 Eu 
half-lives at stellar temperature is suggested to assess the s 
contribution to 152,154 Gd. 

Marginal effects on 152,154 Gd are produced by the 
MACS uncertainties. Wisshak et al. (2006) and Marrone et 
al. (2006) have provided important constraints on the branch 
at 151 Sm (3031 ± 69 mbarn and 3080 ± 150 mbarn, respec¬ 
tively). In our calculations we adopt a weighted average of 
the two measurements. The MACS of 152,154 Gd are well de¬ 
termined with less than 3% of uncertainty (KADoNiS; Bao 
et al. 2000). The 15% uncertainties of the measured MACS 
of the stable isotopes 152,154 Eu and the theoretical MACS 
of lo3 Gd are somewhat larger, but have no significant effect 
on the predictions for 152,154 Gd (<1-2%, respectively). 
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Figure B12. Schematic representation of the chart of nuclei in 
the atomic mass region close to the s-only isotope 176 Lu (red 
rectangle). 176 Lu has a short-lived isomeric state (m) with t -±/ 2 = 
3.66 h and a long-lived ground state (g) with ti /2 = 38 Gyr. Dur¬ 
ing the 13 C(a, n) ie O neutron irradiation these states are treated 
independently; during 22 Ne(o!, n) 25 Mg neutron burst thermally 
induced transitions between isomer and ground states are allowed, 
largely affecting the 176 Lu/ 176 Hf ratio (see text). 
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Figure B13. The same as Fig. lB4l but for the isotopic abun- 
dances of 176 ’ 177 H f, 175,176,177 Lu and 175,176 Yb . 


An additional uncertainty comes from the 22 Ne(a, 
n) 25 Mg reaction: increasing the recommended rate by a fac¬ 
tor of four, the s contribution to solar 152 Gd decreases to 
69%, whereas 154 Gd is affected by only ~4%. 


B3.2 The s-only isotopes ll6 Lu and 176 Hf 

The 176 Lu/ 176 Hf ratio is largely affected by the branch at 
176 Lu, which is strongly sensitive to the temperature during 
TPs (Takahashi & Yokoi 1987; Klay et al. 1991; Doll et al. 
1999). 

176 Lu has a short-lived isomer (£™ 2 = 3.66 h) and an 
extremely long-lived ground state (£ 9 , 2 = 38 Gyr; see 
Fig. IB12I) . In this case, internal transitions are highly for¬ 
bidden by nuclear selection rules and any coupling can only 


be provided by states with intermediate quantum numbers 
and higher excitation energies. Because the low tempera¬ 
tures during the 13 C-pocket phase are not sufficient for that 
coupling, the s-process flow via both states has been treated 
independently. Accordingly, the production of 176 Lu 9 is de¬ 
termined by the partial cross section of 175 Lu (1 — IR = 1 
— 0.84; Wisshak et al. 2006a). 

At temperatures above Ts = 1.5, thermally induced transi¬ 
tions by the hot stellar photon bath start to reach the first 
known mediating states at an excitation energy of 838.6 keV. 
As this mediating state can decay to the long-lived ground 
state as well, this fraction contributes to the production of 
176 Lu 9 and to an enhancement of the (n, 7 ) branch bypass¬ 
ing 1 , 6 Hf. Above Tg = 3, ground and isomeric states are in 
thermal equilibrium, and internal transitions are much faster 
than the timescales for f}~ decay and neutron capture. 

176 Lu 9 is actually produced via thermal transitions in the 
bottom layers of the advanced TPs, where temperatures are 
Ts ~ 3. Once produced, the long-lived 176 Lu 9 survives in the 
cooler external layers of the convective flashes, outside of 
the burning zone. Accordingly, the detailed neutron density 
and temperature profiles in the TPs have to be considered 
for the s-process calculations. For an extensive discussion of 
the treatment of the branch at 176 Lu in our AGB models 
see Heil et al. (2008c). I 11 the advanced TPs, the tempera¬ 
ture reached at the bottom layers is high enough to provide 
additional feeding of the 176 Lu long-lived ground state via 
thermal population of the mediating state. Consequently, at 
the beginning of the TP 176 Hf is destroyed, while 176 Lu is 
produced (+23%; see Fig. IB 131) . The abundance of 1,c Hf 
produced in the previous 13 C pocket is restored to 93% at 
the end of the TP (Fig. IB 131) . 

Note that the branch at 175 Yb (£ 1/2 = 4.18 d, reduced 
up to a factor of three at stellar temperature) has marginal 
effects on the production of 176 Hf, because almost all the s 
path coming from 174 Yb /3“decays into 175 Lu (99%). 

Both 176 Lu and its daughter 1,6 Hf are of pure s-process 
origin, because they are shielded against the r process by 
their stable isobar 176 Yb. About 10% of the long-lived 176 Lu 
ground state decays in the interstellar medium prior to 
the formation of the solar system, with a negligible increase 
of 1 , 6 Hf, which is about ten times more abundant. There¬ 
fore, s contributions between 100 % and 110 % were found 
acceptable for 176 Lu, and between 95% and 105% for 176 Hf 
(Heil et al. 2008c). 

With respect to the results of Heil et al. (2008c), the 
MACS of 1,e Hf has been reduced by 5% due to a new mea¬ 
surement (CTi76 H f) = 594 ± 16 mbarn at 30 keV, Wisshak 
et al. 2006b) and to a revision of the SEF value (KADo- 
NiS). Also the adopted solar abundance ratio by Lodders 
et al. (2009) is ~4% higher than previous assumed (Anders 
and Grevesse 1989). Moreover, the recommended 22 Ne(a, 
n) 25 Mg rate (about a factor of two smaller than previously 
adopted) increases the s abundance of 1,e Lu from 101% to 
112%, and reduces the s prediction of 176 Hf from 104% to 
101%. The increase found for 176 Lu is clearly exceeding the 
5% uncertainty of the solar abundance given by Lodders et 
al. (2009). Otherwise, 176 Lu is particularly sensitive to the 
22 Ne(a, n) 25 Mg rate during TP: an increase of the recom¬ 
mended 22 Ne(a, n) 25 Mg rate by a factor of four is underpro¬ 
ducing the solar 176 Lu abundance by about 8 %, while the s 
abundance of 176 Hf remains almost unchanged. This high- 
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Figure B14. Schematic representation of the chart of nuclei in 
the atomic mass region close to the s-only isotope 204 Pb (red 
rectangle). The branch at 204 T1, mainly activated during TPs, is 
represented with a red circle. The 205 T1 abundance is influenced 
by the radiogenic /3~*~ decay of the long-lived 205 Pb (ti /2 = 17.3 
Myr; blue arrow). At Ts = 3, the stable 205 T1 becomes unstable, 
while the half-life of the long-lived 205 Pb is strongly reduced (see 
text). 
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Figure B15. The same as Fig. lB4l but for the isotopic abun¬ 
dances Of 204,205 p b; 202,203,204 Hg and 203,204,205 TL 


lights that large uncertainties affect the s contribution to 
176 Lu, despite the MACS of the Lu and Hf isotopes are very 
well known, with uncertainties of 1% to 3% (KADoNiS). 

Note that Mohr et al. (2009) reported a significantly 
stronger coupling between isomer and ground state than 
adopted here. Therefore, the longer-lived ground state is 
largely populated, thus strongly increasing the overproduc¬ 
tion of 176 Lu. As evidenced by Mohr et al. (2009), it is not 
possible to find a consistent solution with our AGB models 
within the given experimental errors of the neutron capture 
cross sections of the Lu and Hf isotopes, and the uncertainty 
of the thermal coupling. Mohr et al. (2009) suggested that a 
possible solution of the problem may originate from rather 
small modifications of the temperature profile or convective 
mixing during the helium shell flashes. Cristallo et al. (2010) 
show that convective velocities have to be reduced by a fac¬ 
tor of 1000 with respect to the value estimated by means 
of the mixing length theory, to lead to a lower Lu produc¬ 
tion and to a lower Hf destruction. This strong discrepancy 
seems incompatible with 3D hydrodynamical calculations of 
shell He flashes, which justify uncertainties of a factor 2-3 
in the determination of convective velocities. The solution 
may be found in the nuclear coupling scheme between the 
populated levels in 1,6 Lu (Gintautas et al. 2009; Dracoulis 
et al. 2010; Gosselin et al. 2010). At present, the discrepancy 
between experimental data and astrophysical determination 
remains to be solved. 

B3.3 The s-only isotope 204 Ph (the branch at 204 Tl) 

The most important branch point affecting the s-only nu¬ 
clide 204 Pb is the unstable isotope 204 T1 with a temperature- 
sensitive decay rate (Takahashi & Yokoi 1987; Fig. IB14I) . 
Although 204 Pb is partly produced during the 13 C(a, n) 16 0 
irradiation, most of the AGB s contribution derives from 
the 22 Ne(a, n) 25 Mg neutron source (Ratzel et al. 2004; 
Domingo-Pardo et al. 2007). 

During the 13 C-pocket phase at Tg ~ 0.9, only about 25% of 
204 Pb are produced because of the branching at 204 T1 (ti /2 
= 3.78 yr; N n > 10 7 cm" 3 , /„ > 0.1). 

With increasing stellar temperature, the 204 T1 half-life is 


strongly reduced (fi/ 2 ( 204 Tl) = 7 d at Tg = 3; Takahashi & 
Yokoi 1987), and most of the s path feeds 204 Pb. At the peak 
neutron density 204 Pb is slightly depleted by neutron cap¬ 
ture (N n ^ 2xl0 9 cm -3 , f n ^ 0.2), but is recovered again 
as soon as the neutron density decreases (Fig. IB15I) . At the 
end of the TP, the abundance of 204 Pb has increased by 9%. 

The branch at 203 Hg has no significant impact on the 
204 Pb abundance, because the half-life is strongly reduced 
to less than one day during TPs (Takahashi & Yokoi 1987). 

Previous AGB results (see filled circles in Fig. 5) pre¬ 
dict that about 87% of solar 204 Pb is produced by the main 
component. This value increases to 91% by including an im¬ 
proved treatment of the /3-decay rates of nearby unstable iso¬ 
topes ( 203 Hg and 204 T1) over the convective TPs (see filled 
diamonds in Fig. 5). Accordingly, it is well compatible with 
the solar Pb abundance, which is given with a 10% uncer¬ 
tainty (Lodders et al. 2009). 

The MACS of 204 Pb is well known (81.0 ± 2.3 mbarn with 
~3% uncertainty at 30 keV; KADoNiS). The major un¬ 
certainties affecting the 204 Pb s-prediction derive from the 
branch at 204 T1: the 204 T1 MACS is evaluated theoretically 
with ~18% uncertainty (215 ± 38 mbarn); moreover, Goriely 
(1999) estimated that the 204 Tl (3~ -decay rate may increase 
up to a factor of two, increasing the s contribution to 204 Pb 
by a few percentages. 

Recently, Gonzalez (2014) suggests that the current solar 
log(Pb/Sm) ratio may be overestimated by +0.1 dex, owing 
to a systematic error (related to a trend with condensation 
temperature) that affects the evaluation of the Pb/Sm solar 
abundance. This would reconcile the present s prediction of 
204 Pb. 

In spite of its long terrestrial half-life of 17.3 Myr, 
205 Pb acts as an additional branching point because of 
the strong dependence on temperature and electron density 
(tjy 2 ( 205 Pb) = 5 yr at Tg = 3; Takahashi & Yokoi 1987). 
Similarly, the stable daughter isotope 205 T1 becomes unsta¬ 
ble during TPs (tr /2 ( 205 Tl) = 0.55 yr at Tg = 3; Takahashi 
& Yokoi 1987) and its /3~ decay is competing with the (3 + 
decay of 205 Pb. As the peak neutron density decreases, the 
205 Pb decay starts to prevail, however, and at the end of the 
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Figure B16. Schematic representation of the chart of nuclei in 
the atomic mass region close to 164 Er. The stable isotope 163 Dy 
becomes unstable during TPs and acts, therefore, as the branch 
point isotope responsible for the s production of 164 Er (see text). 
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Figure BIT. The same as Fig. lB4l but for the isotopic abun¬ 
dances Of 164,165,166 Eri 163,165,166 Ho and 163,164 Dy 


TP-AGB phase the radiogenic decay of 205 Pb is responsible 
for the production of 205 T1 (Fig. IB15I see also Mowlavi et 
al. 1998, Busso et al. 1999). 


B4 Additional branches mainly sensitive to stellar 
temperature (and/or electron density) 

B^.l The s contribution to 164 Er (the branches at 163 Dy 
and 163,164 Ho/ 

The origin of 164 Er was first understood by Takahashi and 
Yokoi (1987), who showed that the stable isotope 163 Dy be¬ 
comes unstable under stellar conditions. 

164 Er is completely bypassed during the 13 C-pocket phase 
when the 13 C(a, n) 16 0 operates radiatively at Tg ~ 0.9, be¬ 
cause 163 Dy is still stable at that temperature. During TPs, 
however, 163 Dy becomes unstable and decays to 163 Ho with 
a rapidly decreasing half-life, down to 40 - 120 d (at Tg 
= 3), partially feeding the unstable 163 Ho (Fig. IB161) . Neu¬ 
tron captures on 163 Ho are efficiently activated during TPs. 
Although the half-life of 163 Ho decreases from 4.57 kyr to 



Figure B18. Schematic representation of the chart of nuclei in 
the atomic mass region close to the s-only isotopes 12 8,i30x e . 
The s abundance of 128 Xe is determined only by the competi¬ 
tion between /3~ decays and electron captures of 128 I (red circle), 
completely independent of the s-process neutron density. 


values of 7 to 20 yr at Ts = 3 (Takahashi and Yokoi 1987) it 
remains larger than the time scale for neutron captures dur¬ 
ing the 22 Ne(a, n) 25 Mg irradiation so that the s path runs 
almost completely from 163 Ho to 164 Ho, which decays via 
f}~ to 164 Er and via /3+ to 164 Dy (Fig. IbTFI) . At the end of 
the neutron irradiation, the 163 Ho abundance stored during 
a TP (3 + decays into 163 Dy during the following interpulse 
period, when the s flow bypasses 164 Er. 

The improved treatment of the half-life of 163 Dy and 
163,164 jj 0 over the convective TP provides an increase of so¬ 
lar 164 Er from 81% to 87% (see filled circle and diamond in 
Fig. 5). 

This contribution is however uncertain. Although the 164 Er 
MACS is rather well determined, with 4.7% uncertainty 
(1084 ± 51 mbarn at 30 keV; KADoNiS), it is affected by 
a non-negligible theoretical stellar enhancement factor, SEF 
= 1.08 at 30 keV (KADoNiS; see also Rauscher et al. 2011). 
Moreover, the competition between the 164 Ho /3~ and /3 + - 
decay rates is influenced by the large uncertainty of the /3~- 
decay rate under stellar conditions. Goriely (1999) estimated 
that the 164 Ho /3 ~-decay rate may change up to ~30-40% at 
Ts = 3, producing changes from 75% up to 105% of the 164 Er 
s-prediction. The /3 + decay of 164 Ho has only a marginal ef¬ 
fect, because there are no unknown transition probabilities 
in this case, and the /3~ decay of 163 Dy is also of minor im¬ 
portance (a 30% uncertainty of the decay rate produces up 
to 4% variations of 164 Er). 

This holds for the 22 Ne(a, n) 25 Mg rate as well. An in¬ 
crease of the recommended 22 Ne(a, n) 25 Mg rate by a factor 
of four produces up to 4% variations in the s abundance of 
164 Er. 

B4-2 The s-only pair 128,130 Ae (the branch at 128 1) 

The s predictions of 128,130 Xe depend mainly on the branch 
at 128 I, whereas the short 127 Te half-life (ti /2 = 9.34 h) pre¬ 
vents neutron captures towards 128 Te (Fig. IB181) . 

128 I partly /3~ decays into 128 Xe and partly feeds 128 Te via 
electron captures (tf, 2 = 26.55 m; tff% ~ 7 h). While the /3~ 
half-life shows negligible variations under stellar conditions 
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(ij^ 2 = 42 m at Tg = 3), electron captures on 128 Xe are 
strongly temperature and electron density dependent (e.g., 
at Tg = 3, increases to 10 d for density of p = 1000 
g/cm 3 ; Takahashi & Yokoi 1987). Note that the decay of 
128 I is dominated by the /3~ channel, both during the 13 C 
pocket at Tg ~ 0.9 and during the convective TP at Tg ~ 3, 
so that only 6% of the s path bypasses 128 Xe. 

The branching at 128 1 is unique for the fact that it is com¬ 
pletely independent of the neutron density and that the s 
abundance of 128 Xe is only determined by the competition 
between /3~ decay and electron captures. 

As outlined by Reifarth et al. (2004), the branch at 128 1 
is a good indicator of the convective mixing timescale during 
He shell flashes. Although decays dominate over electron 
captures, the EC half-life of 128 1 decreases to ~6 h in the 
cooler layers of the TP (at Tg = 0.5), directing about 5 to 6% 
of the s path to 128 Te, thus bypassing 128 Xe. By assuming 
a sufficiently fast time scale for turnover mixing in the TP, 
freshly produced 128 I is rapidly transported from the hot 
bottom layers to the cooler external layers of the convective 
zone, allowing a partial activation of the EC channel. The 
resulting ratio of the s abundances of 128 Xe and 130 Xe is 
in agreement with the 128 Xe/ 130 Xe ratio observed in SiC 
grains. 

According to the present s predictions 89% of solar 
128 Xe and 98% of solar 130 Xe are produced by the main 
component. 

The MACS of 128,130 Xe are well determined, with less than 
2% uncertainty (cr( 128 Xe) = 262.5 ± 3.7 mbarn and <r( 130 Xe) 
= 132.0 ± 2.1 mbarn; KADoNiS). Both values are consistent 
with the uncertainty of the solar Xe abundance, which was 
estimated theoretically by interpolation of measured cross- 
sections and abundances of neighbouring elements, through 
the relation crN = constant (see Reifarth et al. 2002; Lod- 
ders et al. 2009). The above s predictions are practically in¬ 
dependent of the 22 Ne(a, n) 25 Mg rate because of the small 
influence of the stellar neutron flux: by increasing the rec¬ 
ommended 22 Ne(a, n) 25 Mg reaction by a factor of four, 
i28,i30-^- e s _p rec hctions show marginal variations (<6%). Ac¬ 
cordingly, the uncertainties of the theoretical MACS values 
of 12 'Te and 128 I of up to a factor of two are affecting the 
128 Xe/ 130 Xe ratio by less than 2%. 
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APPENDIX C: 

The aim of this paper is to provide an overview of the ma¬ 
jor uncertainties that affect the main component: we have 
focused the analysis on two AGB models (M = 1.5 and 3 
Mq at [Fe/H] = —0.3), selected in the mass and metal- 
licity ranges useful to reproduce the solar s distribution 
of isotopes between 90 ^ A ^ 204. However, the study of 
intermediate-mass or low-metallicity AGB models is essen¬ 
tial to understand the s-process nucleosynthesis in differ¬ 
ent environments (e.g., globular clusters, dwarf spheroidal 
galaxies, intrinsic or extrinsic peculiar stars showing s en¬ 
hancement; see Section 1). 

To this purpose, we analyse in this Section two AGB mod¬ 
els chosen as illustrative of extended mass and metallicity 
ranges: a half-solar metallicity 5 Mq model and a 3 Mq 
model at [Fe/H] = — 1. 

The post-process 3 and 5 Mq models are based on stel¬ 
lar input data (e.g., temperature during TPs, TDU, and 
He-intershell masses) of full evolutionary FRANEC models 
(Straniero et al. 2003, 2000). The maximum temperature 
at the bottom of the intershell grows slightly as the mass of 
the H-exhausted core increases (Iben & Renzini 1983): larger 
core masses are attained by more massive stars and, for a 
given mass, by more metal-poor stars (see, e.g., Straniero et 
al. 2006). This allows an efficient activation of the 22 Ne(a, 
n) 2S Mg neutron source in these stars. 

In IMS stars, the mass of the He intershcll is about E 

a factor of ten lower than low-mass models, the efficiency 
of the TDU decreases by about one order of magnitude or 
more, and the interpulse phase is much shorter (~6.5xl0 3 " 

yr). The 13 C(a, n) 16 0 neutron source is expected to have $ 

a small or even negligible effect in these stars. Thus, IMS fo 

AGB stars play a minor role in the Galactic enrichment of § 

s isotopes, but they are crucial for globular clusters (e.g., <= 

D’Orazi et al. 2013). The major uncertainty of IMS stars 
is the treatment of mass loss. In our IMS models we have CC 

adopted an efficient mass loss, but we can not exclude dif¬ 
ferent prescriptions (see Ventura & Marigo 2010, Karakas et 
al. 2012, Straniero et al. 2014). 

Our 5 Mq model at [Fe/H] = 0.3 is characterised by an 
He-intershell mass of ~0.003 Mq, a TDU mass of ~0.001 
Mq and experiences 24 TDU episodes. The maximum tem¬ 
perature at the bottom of the TPs reaches Ts ~ 3.6. The 
dominant neutron source is the 22 Ne(a, n) 25 Mg reaction, 
which is efficiently activated and produces high peak neu¬ 
tron density (N„ ~ 10 11 cm -3 ). 

In low-mass models of close-to-solar metallicity, the 
maximum temperature at the bottom of TPs, although grad¬ 
ually increasing with the pulse number, barely reaches Ts 
= 3. The 3 Mq models at [Fe/H] = — 1 is instead charac¬ 
terised by a rapidly increasing temperature at the bottom 
of the convective zone, whose maximum reaches Ts = 3.5 
in the advanced TPs inducing a considerable production of 
the nuclei up to zirconium (and of the few neutron-rich iso¬ 
topes involved in those s-process branchings) via the 22 Ne(a, 
n) 2S Mg neutron burst. As for the half-solar metallicity mod¬ 
els, the mass loss is considered through the Reimers for¬ 
mula (rj = 1), which allows the model to experience 35 TDU 
episodes. 

We refer to Bisterzo et al. (2010, 2014) for more details on 
these two AGB models. 


M = 5 M 0 : [Fe/H] = -0.3 



Figure Cl. The s distribution of a M = 5 Mq model at half-solar 
metallicity ( top panel) obtained with the recommended 22 Ne(a, 
n) 25 Mg and 22 Ne(a, y) 26 Mg rates (black symbols), compared 
to the results computed with the lower and upper limits of the 
22 Ne(ct, n) 25 Mg rate (red and blue symbols) and with the recom¬ 
mended 22 Ne(a, n) 25 Mg x4 (green symbols), while the 22 Ne(a, 
7) 26 Mg reaction is unchanged. 


All isotopes with 90 < A < 205: 
M = 5 M 0 ; [Fe/H] = -0.3 



Atomic mass 

Figure C2. Ratios between the s distribution of a M = 5 
Mq model obtained with the recommended 22 Ne(ct, n) 25 Mg and 
22 Ne(«, 7 ) 26 Mg rates (black symbols) compared to the results 
computed with the 22 Ne(o, 7 ) 26 Mg rate upper limit (blue sym¬ 
bols). 

In Fig. ED the s distribution of a 5 Mq model at [Fe/H] 
= —0.3 is shown. The results obtained with the recom¬ 
mended 22 Ne(a, n) 25 Mg and 22 Ne(a, 7 ) 26 Mg rates (black 
symbols) are compared with those computed with the lower 
and upper limits of the 22 Ne(a, n) 25 Mg rate (red and blue 
symbols) and with the test 22 Ne(a, n) 25 Mg x4 (which corre¬ 
sponds to Test A in Section 4.1; green symbols). No changes 
are applied to the 22 Ne(a, 7 ) 26 Mg reaction. All values are 
normalised to 150 Sm to overcome the major uncertainties of 
AGB models. 

Larger variations are seen for isotopes with A < 100. How¬ 
ever, the whole s distribution is widely affected: the s- 
production factors increase by up to a factor of 1.7 with 
the upper limit of the recommended 22 Ne(a, n) 25 Mg rate, 
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and by one order of magnitude with Test A. The 22 Ne(a, 
n) 25 Mg neutron burst mostly feeds isotopes between 80 < A 
< 100 (with abundances about a factor of ten higher than 
that of lso Sm), and with a minor extent up to the neutron 
magic isotopes with N = 82. Note that the 150 Sm abun¬ 
dance (in mass fraction) produced by the 5 Mq model (with 
the recommended 22 Ne(a, n) 25 Mg rate) is about 20 times 
higher than its initial (solar-scaled) value. In general, the s- 
predictions obtained by IMS stars are almost negligible with 
respect to the s contribution predicted by low-mass AGB 
models, which produce values of one thousand or more. Ex¬ 
ceptions are a few neutron-rich isotopes, 86 Kr, 87 Rb, 96 Zr: 
they are efficiently produced by the branches at 85 Kr and 
95 Zr, which are essentially open during TPs. 

The competition between the 22 Ne(a, n) 2o Mg and 
22 Ne(a, 7 ) 26 Mg rates affects mostly the Mg isotopes, while 
the s isotopes heavier than A = 90 show variations smaller 
than ~ 6 % (Tig. IC2I) . 


In Fig. IC3I the s predictions (normalised to 150 Sm) of a 
3 Mq model are shown. We compare the half-solar metallic- 
ity model adopted for the main component with a 1/10 so¬ 
lar metallicity model (top and bottom panels, respectively). 
The results obtained by adopting the recommended 22 Ne(a, 
n) 25 Mg and 22 Ne(a, 7 ) 26 Mg rates are represented by black 
symbols. The s predictions computed with the lower and 
upper limits of the 22 Ne(a, n) 25 Mg rate (red and blue sym¬ 
bols) and with the recommended 22 Ne(a, n) 25 Mg x4 (Test 
A in Section 4.1; green symbols) are also displayed for com¬ 
parison. 

Focusing on black symbols, the large dependence on metal¬ 
licity of the s distribution is evident: because the neutron 
exposure increases with decreasing the iron seeds, the abun¬ 
dances of isotopes with neutron magic numbers N = 50 and 
82 are overcome (thus reducing the whole distribution be¬ 
tween 90 < A < 130 and 140 < A < 204, respectively) and 
208 Pb is progressively produced. 

As discussed for the main component, the 22 Ne(a, n) 25 Mg 
neutron source does not substantially modify the s abun¬ 
dances in the half-solar metallicity 3 Mq model, except in 
a few branchings. In the low-metallicity model, not only 
branches are affected, but the whole s distribution is widely 
modified (see green symbols). 

The competition between the 22 Ne(a, n) 25 Mg and 
22 Ne(a, 7 ) 26 Mg rates affects the Mg isotopes. Marginal vari¬ 
ations are found for s isotopes heavier than A = 90 (<2%) 
by adopting the upper limit of the 22 Ne(a, 7 ) 26 Mg rate 

(Fig. El- 

The 13 C(q, n) 16 0 rate operates efficiently in the M 
= 3 Mq model at [Fe/H] = —1, providing a substantial 
contribution to 208 Pb. The present uncertainty associated 
to the 13 C(a, n) le O reaction has marginal effects on s-only 
nuclei (see Fig. IC5l) . 
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Figure C3. The s distribution of a M = 3 Mq model at [Fe/H] 
= —0.3 and [Fe/H] = —1 ( top and bottom panels respectively) 
obtained with the recommended 22 Ne(a:, n) 25 Mg and 22 Ne(a, 
7 ) 26 Mg rates (black symbols), compared to the results computed 
with the lower and upper limits of the 22 Ne(a, n) 25 Mg rate (red 
and blue symbols) and with the recommended 22 Ne(a, n) 25 Mg 
x4 (green symbols), while the 22 Ne(a, 7 ) 26 Mg reaction is un¬ 
changed. 
















Main s-process Branchings and Neutron Source Uncertainties 


All isotopes with 90 < A < 205: 
M = 3 M 0 ; [Fe/H] = -1 


E 

w 

O 

LO 

O 

"O 

CD 

U) 

"cc 

E 

o 


o 

CO 

cc 


1.08 

1.04 

1.00 

0.98 

0.94 

0.90 


iiiiiiiii|iiiiiiiii|iiiiiiiii|iiiiiiiii|iiiiiiiii|iiiiiiiii|iiiiiiiii|iiiiiiiii|iiiiiiiii|iiiiiiiii|iiiiiiiii|iiiiiiiii 

- New (a, n) recommeded - 

- + New (a, y) upper limit - 



.I.I.I.I.I.I.I.I.I.I.I. 


90 100 110 120 130 140 150 160 170 180 190 200 210 

Atomic mass 


Figure C4. Ratios between the s distribution of a M = 3 
Mq model obtained with the recommended 22 Ne(a, n) 25 Mg and 
22 Ne(a, 7 ) 26 Mg rates (black symbols) compared to the results 
computed with the 22 Ne(a, y) 26 Mg rate upper limit (blue sym¬ 
bols). 


M = 3 M 0 : [Fe/H] = -1 



Figure C5. The s distribution of a M = 3 Mq model at [Fe/H] 
= —1 computed with our adopted 13 C(ct, n) 16 0 rate by Denker et 
al. (1995; black symbols), with the test [ 13 C(o, n) 16 0 ]x2 (blue 
symbols) and [ 13 C(o, n) le O ]/3 (red symbols). 
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